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Abstract 

In this review article, we discuss the current status of particle dark 
matter, including experimental evidence and theoretical motivations. We 
discuss a wide array of candidates for particle dark matter, but focus on 
neutralinos in models of supersymmetry and Kaluza-Klein dark matter in 
models of universal extra dimensions. We devote much of our attention to 
direct and indirect detection techniques, the constraints placed by these 
experiments and the reach of future experimental efforts. 
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1 Introduction 



1.1 Overview 

A great deal of effort has been made since 1687, the year of pubUcation of Issac 
Newton's classic work "Philosophiae Naturalis Principia Mathematica", towards 
explaining the motion of astrophysical objects in terms of the laws of gravitation. 
Since then, the deviations of observed motions from expected trajectories have 
proved very effective in deepening our understanding of the Universe. When- 
ever anomalies were observed in the motion of planets in the Solar system, the 
question arose: should such anomalies be regarded as a refutation of the laws of 
gravitation or as an indication of the existence of unseen (today we would say 
"dark") objects? 

The second approach proved to be correct in the case of the anomalous 
motion of Uranus, which led the French astronomer U. Le Verrier and the En- 
glish astronomer John Couch Adams to conjecture the existence of Neptune, 
eventually discovered in 1846 by J.G. Galle. Conversely, the attempt to ex- 
plain the anomalies in the motion of Mercury as due to the existence of a new 
planet, called Vulcan, failed, and the final solution had to wait for the advent 
of Einstein's theory of general relativity, i. e. the introduction of a more refined 
description of the laws of gravitation. 

The modern problem of dark matter is conceptually very similar to the old 
problem of unseen planets. We observe in large astrophysical systems, with sizes 
ranging from galactic to cosmological scales, some "anomalies" that can only be 
explained either by assuming the existence of a large amount of unseen, dark, 
matter, or by assuming a deviation from the known laws of gravitation and the 
theory of general relativity. 

About ten years ago, Jungman, Kamionkowski and Griest wrote a review of 
supersymmetric dark matter for Physics Reports |H19j . This article, although 
incredibly useful, complete and popular, has gradually become outdated over 
the last decade. With this in mind, we have endeavored to write a new review of 
particle dark matter. As with the Jungman et al. article, our review is intended 
to be suitable for a wide range of readers. It could be used as an introduction for 
graduate students interested in this subject or for more experienced scientists 
whose research focuses in other areas. It is also intended to be a useful reference 
in day-to-day research for particle physicists and astrophysicists actively working 
on the problem of dark matter. Unlike the review by Jungman et at, we do not 
limit our discussion to supersymmetric dark matter. 

The article is organized as follows: we first present, in this chapter, a brief 
review of the Standard Model of particle physics and cosmology, and review our 
present understanding of the history of the Universe. We focus in particular on 
the freeze-out of dark matter particles and on the calculation of their relic abun- 
dance, and discuss the possible relationship between dark matter and physics 
beyond the Standard Model of particle physics. 

Chapter 2 is devoted to the compelling evidence for dark matter at all astro- 
physical length scales. We review the key observations and discuss the theoret- 
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ical predictions (from N-body simulations) for the distribution of dark matter, 
focusing in particular on the innermost regions of galaxies, and discuss how they 
compare with observations. Particular attention is devoted to the galactic cen- 
ter, where the presence of a supermassive black hole could significantly modify 
the dark matter distribution. 

Dark matter candidates are presented in chapter 3. We start with an in- 
troduction to the "dark matter zoo" , i.e. a description of the many candidates 
that have been proposed in the literature. We then focus on two particularly 
interesting dark matter candidates: the supersymmetric neutralino and Kaluza- 
Klein dark matter. For each of these candidates, we give a brief introduction to 
the physical motivations and underlying theories. We conclude chapter 3 with 
a review of the constraints put on dark matter from collider experiments, and 
discuss the prospects for future experiments. 

The second part of this review is devoted to the astrophysical constraints 
on particle dark matter. We begin in chapter 4 with a review of existing and 
next-generation experiments that will probe the nature of dark matter. This 
chapter is propedeutical to chapter 5, which discusses the many possible direct 
and indirect searches of dark matter and which constitutes the heart of this 
review. We give our conclusions in chapter 6. Some useful particle physics 
details are given in the appendices. 

1.2 Standard Cosmology 

Although the exact definition of the Standard cosmological model evolves with 
time, following the progress of experiments in measuring the cosmological pa- 
rameters, most cosmologists agree on a fundamental picture, the so-called Big 
Bang scenario, which describes the Universe as a system evolving from a highly 
compressed state existing around 10^" years ago. 

This picture has its roots in the discovery of Bubble's law early in the past 
century, and has survived all sorts of cosmological observations, unlike alter- 
native theories such as the "steady state cosmology" , with continuous creation 
of baryons, which, among other problems, failed to explain the existence and 
features of the cosmic microwave background. 

We now have at our disposal an extremely sophisticated model, allowing us 
to explain in a satisfactory way the thermal history, relic background radiation, 
abundance of elements, large scale structure and many other properties of the 
Universe. Nevertheless, we are aware that our understanding is still only partial. 
It is quite clear that new physics is necessary to investigate the first instants of 
our Universe's history (see section n~H|l . 

To "build" a cosmological model, in a modern sense, three fundamental 
ingredients are needed: 

• Einstein equations, relating the geometry of the Universe with its matter 
and energy content 

• metrics, describing the symmetries of the problem 
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• Equation of state, specifying the physical properties of the matter and 
energy content 

The Einstein field equation can be derived almost from first principles, as- 
suming that: 1) the equation is invariant under general coordinate transforma- 
tions, 2) the equation tends to Newton's law in the limit of weak fields, and 3) 
the equation is of second differential order and linear in second derivatives |4U(J) . 
The resulting equation reads 



where i?^^ and R are, respectively, the Ricci tensor and scalar (obtained by 
contraction of the Riemann curvature tensor). Qf^i, is the metric tensor. Gat is 
Newton's constant, r^,y is the energy-momentum tensor, and A is the so-called 
cosmological constant. 

Ignoring for a moment the term involving the cosmological constant, this 
equation is easily understood. We learn that the geometry of the Universe, 
described by the terms on the left-hand-side, is determined by its energy con- 
tent, described by the energy-momentum tensor on the right-hand-side. This 
is the well known relationship between the matter content and geometry of the 
Universe, which is the key concept of general relativity. 

The addition of the cosmological constant term, initially introduced by Ein- 
stein to obtain a stationary solution for the Universe and subsequently aban- 
doned when the expansion of the Universe was discovered, represents a "vacuum 
energy" associated with space-time itself, rather than its matter content, and is 
a source of gravitational field even in the absence of matter. The contribution 
of such "vacuum energy" to the total energy of the Universe can be important, 
if one believes recent analyses of type la supernovae and parameter estimates 
from the cosmic microwave background (for further discussion see section 

To solve the Einstein equations one has to specify the symmetries of the 
problem. Usually one assumes the properties of statistical homogeneity and 
isotropy of the Universe, which greatly simplifies the mathematical analysis. 
Such properties, made for mathematical convenience, are confirmed by many 
observations. In particular, observations of the Cosmic Microwave Background 
(CMB) have shown remarkable isotropy (once the dipole component, interpreted 
as due to the Earth motion with respect to the CMB frame, and the contribu- 
tion from the galactic plane were subtracted) . Isotropy alone, if combined with 
the Copernican principle, or "mediocrity" principle, would imply homogene- 
ity. Nevertheless, direct evidence of homogeneity comes from galaxy surveys, 
suggesting a homogeneous distribution at scales in excess of ~ 100 Mpc. More 
specifically, spheres with diameters larger than ^ 100 Mpc centered in any place 
of the Universe should contain, roughly, the same amount of matter. 

The properties of isotropy and homogeneity imply a specific form of the 
metric: the line element can in fact be expressed as 



(1) 
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where a{t) is the so-called scale factor and the constant k, describing the spatial 
curvature, can take the values k = —1, 0, +1. For the simplest case, fc = 0, the 
spatial part of Eq. |21 reduces to the metric of ordinary (flat) Euclidean space. 

The Einstein equations can be solved with this metric, one of its components 
leading to the Friedmann equation 

, k SttGn ,„x 
^ — Ptot, (3) 



,2 



where ptot is the total average energy density of the universe. It is common to 
introduce the Hubble parameter 

A recent estimate 1 159j of the present value of the Hubble parameter, Hq, 
(also referred to as the Hubble constant) is Hq = 73 ± 3 km s~^ Mpc^^. We see 
from Eq. 13 that the universe is flat (fc = 0) when the energy density equals the 
critical density, pc- 

Pc ^ (5) 

In what follows we will frequently express the abundance of a substance in the 
Universe (matter, radiation or vacuum energy), in units of pc- We thus define 
the quantity fi^ of a substance of species i and density pi as 

n.^^ . (6) 

Pc 

It is also customary to define 

Pc 

in terms of which the Friedmann equation (Eq. ^ can be written 

"-'-ift' (8) 

The sign of k is therefore determined by whether is greater than, equal to, or 
less than one (see table 1). 

Following Ref. [HE] , we note that the various fii evolve with time differently, 
depending on the equation of state of the component. A general expression for 
the expansion rate is 



Hi 



+nM (1 + zf + flail + zf (9) 
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P< Pc 


n < 1 


k = -i 


open 


Pc 


n = 1 




flat 


P> Pc 


n > 1 


k = 1 


closed 



Table 1: Classification of cosmological models based on the value of the average density, 
p, in terms of the critical density, p^. 



where M and R are labels for matter and radiation, VLk = -i !^-2 and X refers 
to a generic substance with equation of state px — cexPx (in particular, for the 
cosmological constant, = —1)- z is the redshift. 

We discuss in Sec. 12.31 recent estimates of cosmological parameters using 
CMB measurements, combined with various astrophysical observations. 



1.3 The Standard Model of Particle Physics 

The Standard Model (SM) of particle physics has, for many years, accounted 
for all observed particles and interactions ^. Despite this success, it is by now 
clear that a more fundamental theory must exist, whose low-energy realization 
should coincide with the SM. 

In the SM, the fundamental constituents of matter are fermions: quarks and 
leptons. Their interactions are mediated by integer spin particles called gauge 
bosons. Strong interactions are mediated by gluons Ga, electroweak interaction 
by W^, Zq, 7 and the Higgs boson H*^. The left-handed leptons and quarks are 
arranged into three generations of SU (2) l doublets 



d'). ( s' ). { b' 



(10) 



(11) 

L \ " / L \ " / L 

with the corresponding right-handed fields transforming as singlets under SU (2)l. 
Each generation contains two flavors of quarks with baryon number 5 = 1/3 
and lepton number L = and two leptons with B — and L ^ I. Each particle 
also has a corresponding antiparticle with the same mass and opposite quantum 
numbers. 

The quarks which arc primed are weak eigenstates related to mass eigenstates 
by the Cabibbo-Kobayashi-Maskawa (CKM) matrix 




= VbKM \ s \. (12) 



^ It is a matter of definition whether one considers neutrino masses as part of the SM or as 
physics beyond the SM. 
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Gauge symmetries play a fundamental role in particle physics. It is in fact in 
terms of symmetries and using the formalism of gauge theories that we describe 
electroweak and strong interactions. The SM is based on the SU{3)c(^SU (2)i(g) 
U{1)y gauge theory, which undergoes the spontaneous breakdown: 

SU{3)c SU{2)l U{1)y SU{3)c ® U{1)q (13) 

where Y and Q denote the weak hypercharge and the electric charge genera- 
tors, respectively, and 5[/(3)c describes the strong (color) interaction, known 
as Quantum ChromoDynamics (QCD). This spontaneous symmetry breaking 
results in the generation of the massive and Z gauge bosons as well as a 
massive scalar Higgs field. 

1.4 A very brief history of the Universe 

Our description of the early Universe is based on an extrapolation of known 
physics back to the Planck epoch, when the Universe was only t = 10~'''^seconds 
old, or equivalently up to energies at which the gravitational interaction becomes 
strong (of the order of the Planck mass, Mpi = 10^^ GeV). Starting at this epoch 
we take now a brief tour through the evolution of the Universe: 

• T ^ IQi^ GcV. It is thought that at this scale, some (unknown) grand 
unified group, G, breaks down into the Standard Model gauge group, 
SU {3)c SU (2) L (E)U {1)y ■ Little is known about this transition, however. 

• T 10^ GeV. The Standard Model gauge symmetry breaks into SU{3)c'Si 
U{1)q (see Eg. I13f) . This transition, called electroweak symmetry break- 
ing, could be the origin of baryogenesis (see e.g. Ref. ^|) and possibly of 
primordial magnetic fields {e.g. Ref. |317p . 

• T ^ 10^ — 10'^ GeV. Weakly interacting dark matter candidates with GeV- 
TeV scale masses freeze- out^ as discussed in next section. This is true in 
particular for the neutralino and the B^^^ Kaluza-Klein excitation that we 
discuss in Chapter O 

• T ~ 0.3 GeV. The QCD phase transition occurs, which drives the confine- 
ment of quarks and gluons into hadrons. 

• T ^ 1 MeV. Neutron freeze-out occurs. 

• T ^ 100 keV. Nucleosynthesis: protons and neutrons fuse into light ele- 
ments (D, ■^He, ^He, Li). The standard Big Bang nucleosynthesis (BBN) 
provides by far the most stringent constraints to the Big Bang theory, and 
predictions remarkably agree with observations (see Fig.^. 

• T ~ 1 eV. The matter density becomes equal to that of the radiation, 
allowing for the formation of structure to begin. 

• T ^ 0.4 eV. Photon decoupling produces the cosmic background radiation 
(CMB), discussed in Sec.O 
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• T = 2.7K ~ 10-'' eV. Today. 



1.5 Relic Density 

We briefly recall here the basics of the calculation of the density of a thermal 
relic. The discussion is based on Refs. |34UI 12601 1447) and we refer to them for 
further comments and details. 

A particle species in the early Universe has to interact sufhciently or it 
will fall out of local thermodynamic equilibrium. Roughly speaking, when its 
interaction rate drops below the expansion rate of the Universe, the equilibrium 
can no longer be maintained and the particle is said to be decoupled. 



10 



1.5.1 The standard calculation 

The evolution of the phase space distribution function, /(p,x), is governed by 
the Bohzmann equation 

m = c[fi (14) 

where L is the LiouviUe operator, and C is the cohision operator, describing the 
interactions of the particle species considered. 

After some manipulation, the Boltzmann equation can be written as an 
equation for the particle number density n: 

J + 3Hn - -{av) {n^ - (n^'^ f) , (15) 

where av is the total annihilation cross section multiplied by velocity, brackets 
denote thermal average, H is Hubble constant, and n'^'^ is the number density 
at thermal equilibrium. For massive particles, i.e. in the non-relativistic limit, 
and in the Maxwell-Boltzmann approximation, one has 

n^^ = ff(^)'^'e-"/^, (16) 

where m is the particle mass and T is the temperature. We next introduce the 
variables 

Y=- , y^« = — (17) 

s s 

where s is the entropy density s = 2Tr^g^T^ /A5 and counts the number of 
relativistic degrees of freedom. Using the conservation of entropy per co-moving 
volume {sa^ =constant), it follows that n + 3Hn — sY and Ea. [TBI reads 

If we further introduce the variable x = ^, Eq. ^1 can be expressed as 

For heavy states, we can approximate (av) with the non-relativistic expansion 
in powers of v'^ 

(av) = a + b{v^) + ©((w")) ^a + 6b/x , (20) 
which leads to our final version of Eq.^lin terms of the variable A = Y — Y^'^: 
A' = -y-"' - /(a;)A(2r"« + A), (21) 
where prime denotes d/dx and 



f{x) = J^m Mpi (a -f 6 h/x) x-\ (22) 
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Following Ref. |34(J| we introduce the quantity xp = m/TF, where Tp is the 
freeze-out temperature of the relic particle, and we notice that Eq. 1211 can be 
solved analytically in the two extreme regions x <ti xp and x ^ xp 

^= ~ 2/fa;7y-^ forx<t:xp (23) 
A' = -,fix)A^ for x > xp. (24) 

These regions correspond to long before freeze-out and long after freeze-out, 
respectively. Integrating the last equation between x p and oo and using A^p 2> 
Aoo, we can derive the value of A^o and arrive at 



Y^^^^^Mpimxp\a + 3b/xp). (25) 

The present density of a generic relic, X, is simply given by px — mxnx = 
'TijfSoFoo, where sq = 2889.2 cm~'^ is the present entropy density (assuming 
three Dirac neutrino species). The relic density can finally be expressed in 
terms of the critical density (see Eq. O 

^ ,2 1.07 X lO^GeV"^ xp 1 

^xh^ ~ —-, ; — r, (26) 

Mpi ^{a + 3b/xp)' ^ ^ 

where a and b are expressed in GeV"'^ and is evaluated at the freeze-out 
temperature. It is conventional to write the relic density in terms of the Hubble 
parameter, h — iJo/lOOkms^^ Mpc~^. 

To estimate the relic density, one is thus left with the calculation of the 
annihilation cross sections (in all of the possible channels) and the extraction 
of the parameters a and b, which depend on the particle mass. The freeze-out 
temperature x p can be estimated through the iterative solution of the equation 



Xp = hi 




m Mpi{a + 6b/ Xp) 

1/2 1/2 
5* p 



(27) 



where c is a constant of order one determined by matching the late-time and 
early-time solutions. 

It is sometimes useful to perform an order-of-magnitude estimate using an 
approximate version of Eq. (319^ : 

^xh^~ ■ 28 

{(Tv) 

We note that the approximation introduced in Eq. [^Olis not always justified 
(see e.g. Ref. |319p . For example, Ref. |437| suggests a scenario where the 
presence of a scalar field in the early Universe could significantly affect the 
value of the relic density. Furthermore, a dramatic change in the relic density 
can be induced by resonance enhancements or so-called coannihilations. We 
discuss the effects of coannihilations in the next section. 
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1.5.2 Including coannihilations 

Following earlier works (see Ref. llOSp . Griest and Seckel |279j noticed that 
if one or more particles have a mass similar to the relic particle and share a 
quantum number with it, the standard calculation of relic density fails. 

Let us consider N particles Xi {i — 1, . . . , N) with masses and internal 
degrees of freedom (statistical weights) gi. Also assume that toi < TO2 < • • • < 
rriN-i < rriN, and that the lightest particle is protected against decay thanks 
to some symmetry (i.e. R-parity or KK-parity, for neutralinos or Kaluza-Klein 
particles, respectively. See section OJ. We will also denote the lightest particle 
byXi. 

In this case, Eg. llSlbecomes 

^ = -3Hn ~ ^ (cTy Uy ) {n.^j - n^n]"^) , (29) 

where n is the number density of the relic particle and n = J2^=i "iue to the 
fact that the decay rate of particles, Xi, other than the lightest is much faster 
than the age of the Universe. Here, 

= ;^a(X,X, ^Xsm) (30) 

X 

is the total annihilation rate for XiXj annihilations into a Standard Model 
particle. Finally, 



v^i = ^ (31) 

is the relative particle velocity, with pi and Ei being the four-momentum and 
energy of particle i. 

The thermal average {uijVij) is defined with equilibrium distributions and 
is given by 

where fi are distribution functions in the Maxwell-Boltzmann approximation. 

The scattering rate of supcrsymmetric particles off particles in the thermal 
background is much faster than their annihilation rate. We then obtain 



— = -3Hn - ((Toffw) ("^ - "cq) 7 (33) 



dn 

where 

^3 



oq „cq 



Tcsv) = }^ (ay ) . (34) 
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Edsjo and Gondolo |2U2| reformulated the thermal average into the more 
convenient expression 

(^cfff) = : (35) 



where Ki are the modified Bessel functions of the second kind and of order i. 
The quantity Wcs is defined as 



(36) 

where Wij = AEiEjaijVij and is the momentum of the particle Xi (or Xj) 
in the center-of-mass frame of the pair XiXj, and s = mf + m| + 2EiEj — 
2\pi\\pj\ cos9, with the usual meaning of the symbols. 

The details of coannihilations in the framework of supersymmetric models 
are well established (see e.g. the recent work of Edsjo et al. 2Q63). and numerical 
codes now exist including coannihilations with all supersymmetric particles, 
e.g. MicrOMEGAs [HHI and the new version of DarkSusy gMIEMI, pubhcly 
released in 2004. The case of coannhilations with a light top squark, such as 
the one required for the realization of the electroweak baryogenesis mechanism, 
has been discussed in R,ef.|HH]. 



1.6 Links with Physics Beyond the Standard Model 

The concepts of dark energy and dark matter do not find an explanation in the 
framework of the Standard Model of particle physics. Nor are they understood in 
any quantitative sense in terms of astrophysics. It is interesting that also in the 
realm of particle physics, evidence is accumulating for the existence of physics 
beyond the Standard Model, based on theoretical and perhaps experimental 
arguments. 

On the experimental side, there is strong evidence for oscillations of atmo- 
spheric neutrinos (originating from electromagnetic cascades initiated by cosmic 
rays in the upper atmosphere) and solar neutrinos. The oscillation mechanism 
can be explained under the hypothesis that neutrinos do have mass, in contrast 
to the zero mass neutrinos of the Standard Model (see Ref. |3f)9j for a recent 
review) . 

On the theoretical side, many issues make the Standard Model unsatisfac- 
tory, for example the hierarchy problem, i.e. the enormous difference between 
the weak and Planck scales in the presence of the Higgs field (this will be dis- 
cussed in some detail in Sec. 13.2. or the problem of unification addressing 
the question of whether there exists a unified description of all known forces, 
possibly including gravity. 

The list of problems could be much longer, and it is natural to conjecture 
that our Standard Model is the low-energy limit of a more fundamental theory. 
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Two examples of popular extensions of the Standard Model include: 

• Super symmetry. As a complete symmetry between fermions and bosons, 
supersymmetry's theoretical appeal is very great | 498| . So great, in fact, 
is this appeal, that it appears to many as a necessary ingredient of future 
extensions of the Standard Model. Many interesting features make it at- 
tractive, including its role in understanding the fundamental distinction 
between bosons and fermions, and the problems of hierarchy and unifica- 
tion discussed above. Last, but not least, it provides an excellent dark 
matter candidate in terms of its lightest stable particle, the neutralino. 
We will present the basics of supersymmetry and the properties of the 
neutralino in Sec. 13.21 

• Extra dimensions. In the search of a fundamental theory with a unified 
description of all interactions, physicists developed theories with extra 
spatial dimensions, following an early idea of Kaluza |322| . who extended 
to four the number of space dimensions to include electromagnetism into 
a "geometric" theory of gravitation. In theories with unified extra dimen- 
sions, in which all particles and fields of the Standard Model can propagate 
in the extra dimensions, the lightest Kaluza-Klein particle, i.e. the lightest 
of all the states corresponding to the first excitations of the particles of 
the Standard Model, is a viable dark matter candidate, as we discuss in 
Sec.rOl 

Despite the fact that neutrinos are thought to be massive, they are essen- 
tially ruled out as dark matter candidates (see Sec. 13.1(1 . Consequently, the 
Standard Model does not provide a viable dark matter candidate. This is fur- 
ther supported by the fact that most of the dark matter is non-baryonic (see 
section lT^ . Dark matter is therefore a motivation to search for physics beyond 
the Standard Model (others might say that this is evidence for physics beyond 
the Standard Model). 

This is a typical example of the strong interplay between particle physics, 
theoretical physics, cosmology and astrophysics. From one side, theoretical par- 
ticle physics stimulates the formulation of new theories predicting new particles 
that turn out to be excellent dark matter candidates. On the other side, cosmo- 
logical and astrophysical observations constrain the properties of such particles 
and consequently the parameters of the new theories. 

2 Evidence and Distribution 

2.1 The Galactic Scale 

The most convincing and direct evidence for dark matter on galactic scales 
comes from the observations of the rotation curves of galaxies, namely the graph 
of circular velocities of stars and gas as a function of their distance from the 
galactic center. 
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Figure 2: Rotation curve of NGC 6503. The dotted, dashed and dash-dotted lines are 
the contributions of gas, disk and dark matter, respectively. From Ref. (5U| . 

Rotation curves are usually obtained by combining observations of the 21cni 
line with optical surface photometry. Observed rotation curves usually exhibit 
a characteristic flat behavior at large distances, i.e. out towards, and even far 
beyond, the edge of the visible disks (see a typical example in Fig. |2J|. 

In Newtonian dynamics the circular velocity is expected to be 



where, as usual, AI{r) = An J p(r)r^dr, and p{r) is the mass density profile, 
and should be falling cx l/\/r beyond the optical disc. The fact that v{r) is 
approximately constant implies the existence of an halo with M(r) oc r and 
p oc 

Among the most interesting objects, from the point of view of the observa- 
tion of rotation curves, are the so-called Low Surface Brightness (LSB) galaxies, 
which are probably everywhere dark matter-dominated, with the observed stel- 
lar populations making only a small contribution to rotation curves. Such a 
property is extremely important because it allows one to avoid the difficulties 
associated with the deprojection and disentanglement of the dark and visible 
contributions to the rotation curves. 

Although there is a consensus about the shape of dark matter halos at large 
distances, it is unclear whether galaxies present cuspy or shallow profiles in their 
innermost regions, which is an issue of crucial importance for the effects we will 
be discussing in the following chapters. 

Using high-resolution data of 13 LSB galaxies, de Blok et al. |179j recently 
showed, that the distribution of inner slopes, i.e. the power-law indices of the 
density profile in the innermost part of the galaxies, suggests the presence of 




(37) 



16 



-2 -1.5 -1 -0.5 0.5 
inner' slope a 



-2-10 1 2 

log(r,„/kpc) 



Figure 3: Left panel: The distribution of inner slopes, a, of dark matter density profiles 
in LSB galaxies. The hatched (blank) histogram represents well-resolved (unresolved) 
galaxies. Right panel: The value of a as a function of the radius of the innermost 
point. From Ref. [H^. 

shallow, or even flat, cores (see Fig. El . Furthermore, the highest values of the 
power-law index are obtained in correspondence to galaxies with the poorest 
resolution, as can be seen from the right panel of the same figure. 

Following Salucci and Borriello |439 | . rotation curves of both low and high 
surface luminosity galaxies appear to suggest a universal density profile, which 
can be expressed as the sum of an exponential thin stellar disk, and a spher- 
ical dark matter halo with a flat core of radius rg and density po = 4.5 x 
lO~^(ro/kpc)~^/^M0pc~^ (here, M© denotes a solar mass, 2 x 10^° kg). In a 
similar way the analysis of Reed et al. 425 leads to the conclusion that sim- 
ulated halos have significantly steeper density profiles than are inferred from 
observations. 

Nevertheless, claims have been made in the literature about the possibil- 
ity of reconciling these results with the steep profiles predicted by numerical 
simulations (see section 12.41 for a discussion on the state of art of N-body sim- 
ulations and for further discussions, see Refs. |179l 14831 1427| ). In particular, 
Hayashi et al. j291| have claimed consistency between most observations and 
their simulated profiles and have argued that the remaining discrepancies could 
be explained by taking into account the difference between the circular velocity 
and gas rotation speed, likely to arise in gaseous disks embedded within realistic, 
triaxial cold dark matter halos. 

Another area of contention is that of the dark matter content in the inner 
halos of massive disk galaxies. It has been argued that barred galaxies cannot 
contain substantial amounts of dark matter out to the outermost extent of the 
observed bars, otherwise the rapidly rotating bars would have slowed down due 
to dynamical friction on the dark matter |177[ 117^ . One counterargument is 
the contention that bars may be dynamically young systems that formed by 
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secular evolution of unstable cold disks and hence poor dynamical probes |158| . 
Another is that the slowing down of bars, perhaps in an earlier phase of the 
forming galaxy, actually heated the dark matter and generated a core. 

Despite the uncertainties of the slope in the innermost regions of galaxies, 
rotation curves of disk galaxies provide strong evidence for the existence of a 
spherical dark matter halo. The total amount of dark matter present is difficult 
to quantify, however, as we do not know to what distances halos extend. Addi- 
tional evidence for dark matter at galactic scales comes from mass modelling of 
the detailed rotation curves, including spiral arm features. Submaximal disks 
are often, although not always, required |455| . 

Some elliptical galaxies show evidence for dark matter via strong gravita- 
tional lensing |341| . X-ray evidence reveals the presence of extended atmo- 
spheres of hot gas that fill the dark halos of isolated ellipticals and whose hy- 
drostatic support provides evidence for dark matter. In at least one case, an 
elliptical contains a cold gas disk whose HI rotation curve is flat out to about 
5 half light radii. In contrast, however, planetary nebula studies to a similar 
distance for other ellipticals can be explained only with a constant mass-to-light 
ratio. There may be some dark matter in these cases, but its relative dominance 
does not appear to increase with increasing galactocentric distance. Rather, it 
is asociated with the stellar distribution. 

Other arguments for dark matter, both on subgalactic and inter-galactic 
scales, also comes from a great variety of data. Without attempting to be 
complete, we cite among them: 

• Weak modulation of strong lensing around individual massive elliptical 
galaxies. This provides evidence for substructure on scales of ~ lO^M© 

• The so-called Oort discrepancy in the disk of the Milky Way (see e.g. Ref. \51\). 
The argument follows an early suggestion of Oort, inferring the existence 

of unobserved matter from the inconsistency between the amount of stars, 
or other tracers in the solar neighborhood, and the gravitational potential 
implied by their distribution. 

• Weak gravitational lensing of distant galaxies by foreground structure (see, 
e.g. Ref. gMj). 

• The velocity dispersions of dwarf spheroidal galaxies which imply mass- 
to-light ratios larger than those observed in our "local" neighborhood. 
While the profiles of individual dwarfs show scatter, there is no doubt 
about the overall dark matter content (see Refs. |486[ 1373] ). 

• The velocity dispersions of spiral galaxy satellites which suggest the exis- 
tence of dark halos around spiral galaxies, similar to our own, extending 
at galactocentric radii ^ 200 kpc, i.e. well behind the optical disc. This 
applies in particular to the Milky Way, where both dwarf galaxy satellites 
and globular clusters probe the outer rotation curve (see Refs. |507[ I46| '). 
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Figure 4: Chandra X-ray (left) and Hubble Space Telescope Wide Field Planetary 
Camera 2 optical (right) images of AbeU 2390 {z = 0.230) and MS2137.3-2353 {z = 
0.313). Note the clear gravitational arcs in the Hubble images. From Ref. |225| . 



2.2 The Scale of Galaxy Clusters 

A cluster of galaxies gave the first hints of dark matter (in the modern sense) . In 
1933, F.Zwicky |51(J| inferred, from measurements of the velocity dispersion of 
galaxies in the Coma cluster, a mass-to-light ratio of around 400 solar masses 
per solar luminosity, thus exceeding the ratio in the solar neighborhood by 
two orders of magnitude. Today, most dynamical estimates 1331L I139| are 
consistent with a value fl^ ~ 0.2—0.3 on cluster scales. A convenient calibration 
is riM = (M/L)/1000. 

The mass of a cluster can be determined via several methods, including 
application of the virial theorem to the observed distribution of radial velocities, 
by weak gravitational lensing, and by studying the profile of X-ray emission that 
traces the distribution of hot emitting gas in rich clusters. 

Consider the equation of hydrostatic equilibrium for a system with spherical 
symmetry 

(38) 
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where P, p, and a are, respectively, the pressure, density, and gravitational 
acceleration of the gas, at radius r. For an ideal gas, this can be rewritten in 
terms of the temperature, T, and the average molecular weight, ^ ~ 0.6, 



where uip is the proton mass. The temperature of clusters is roughly constant 
outside of their cores and the density profile of the observed gas at large radii 
roughly follows a power-law with an index between —2 and —1.5. We then find 
that the temperature should obey the relation 



for the baryonic mass of a typical cluster, where Mr is the mass enclosed within 
the radius r. The disparity between the temperature obtained using Ea. l4Ul and 
the corresponding observed temperature, T s:! 10 keV, when Mr is identified 
with the baryonic mass, suggests the existence of a substantial amount of dark 
matter in clusters. 

These conclusions can be checked against estimates from gravitational lens- 
ing data (see Fig. 0J. Following Einstein's theory of general relativity, light 
propagates along geodesies which deviate from straight lines when passing near 
intense gravitational fields. The distortion of the images of background objects 
due to the gravitational mass of a cluster can be used to infer the shape of 
the potential well and thus the mass of the cluster (see e.g. Ref. |477| for a 
spectacular demonstration of gravitational lensing in clusters) . 

The fraction of baryons inside a cluster, crucial to disentangle the contri- 
butions of ordinary (visible) and dark matter, can also be inferred through the 
so-called Sunyaev-Zel'dovich effect by which the cosmic microwave background 
(see section gets spectrally distorted through Compton scattering on hot 
electrons. 

Despite general agreement between dark matter density profiles at large 
radii and numerical simulations (see section 12.4(1 , it is unclear whether there 
is agreement with the predicted profiles in the cores of clusters. Gravitational 
lensing measurements appear to be in confiict with cuspy profiles, excluding 
at the 99% confidence level cusps with power-law indices of about —1 (see 



This argument is strengthened by use of radial arcs which probe the mass 
gradient, but is weakened if the cluster is not spherically symmetric. Indeed 
an asymmetry of a few percent alows the cluster profiles to be consistent with 
NFW. Moreover, recent Chandra observations of X-ray emission from Abell 2029 
suggest a full compatibility of dark matter distributions with cuspy profiles (see 
Ref. |358j ). For a critique of gravitational lensing constraints on dark matter 
halo profiles, see Ref. I171j . 




(39) 
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e.g. Ref. PD)'!. 
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2.3 Cosmological Scales 



We have seen in the previous sections that, on distance scales of the size of 
galaxies and clusters of galaxies, evidence of dark matter appears to be com- 
pelhng. Despite this, the observations discussed do not allow us to determine 
the total amount of dark matter in the Universe. We discuss in this section how 
such information can be extracted from the analysis of the Cosmic Microwave 
Background (CMB). 

Excellent introductions to CMB theory exist in the literature [3121 \'61'6\ . 
Here, we limit ourselves to a brief review of the implications of recent CMB data 
on the determination of cosmological parameters. In particular, we discuss the 
stringent constraints on the abundances of baryons and matter in the Universe 
placed by the Wilkinson Microwave Anisotropy Probe (WMAP) data. 

The existence of background radiation originating from the propagation of 
photons in the early Universe (once they decoupled from matter) was predicted 
by George Gamow and his collaborators in 1948 and inadvertently discovered 
by Arno Penzias and Robert Wilson in 1965. After many decades of experi- 
mental effort, the CMB is known to be isotropic at the 10~^ level and to follow 
with extraordinary precision the spectrum of a black body corresponding to a 
temperature T = 2.726 K. 

Today, the analysis of CMB anisotropics enables accurate testing of cosmo- 
logical models and puts stringent constraints on cosmological parameters. 

The observed temperature anisotropics in the sky are usually expanded as 



where Yim{0, (f) are spherical harmonics. The variance of a^n is given by 



If the temperature fluctuations are assumed to be Gaussian, as appears to be 
the case, all of the information contained in CMB maps can be compressed into 
the power spectrum, essentially giving the behavior of Ci as a function of 
Usually plotted is t{£+l)Ci/2iT (see Fig.p. 

The methodology, for extracting information from CMB anisotropy maps, 
is simple, at least in principle. Starting from a cosmological model with a fixed 
number of parameters (usually 6 or 7), the best-fit parameters are determined 
from the peak of the N-dimensional likelihood surface. 

From the analysis of the WMAP data alone, the following values are found 
for the abundance of baryons and matter in the Universe 



Taking into account data from CMB experiments studying smaller scales (with 
respect to WMAP), such as ACBAR ^3iS and CBl j411j . and astronomical 




(41) 




(42) 



VLbh^ = 0.024 ± 0.001 VtMh^ = 0.14 ± 0.02. 



(43) 
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Figure 5: CMB Temperature fluctuations: A comparison between COBE and WMAP. 
Image from http://map.gsfc.nasa.gov/. 
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Figure 6: The observed power spectrum of CMB anisotropies. From Ref. |470| . 



measurements of tlie power spectrum from large scale structure (2dFGRS, see 
Ref. |414j ') and the Lyman a forest (see e.g. Ref. |167j ). the constraints become 

17b/i^ = 0.0224 ±0.0009 and f^wh^ = 0.135to^o^. (44) 

The value of flbh^ thus obtained is consistent with predictions from Big Bang 
nucleosynthesis {e.g. \4:0'd\ ) 

0.018 < ribh^ < 0.023. (45) 

Besides those provided by CMB studies, the most reliable cosmological mea- 
surements are probably those obtained by Sloan Digital Sky Survey (SDSS) 
team, which has recently measured the three-dimensional power spectrum, P{k), 
using over 200,000 galaxies. An estimate of the cosmological parameters com- 
bining the SDSS and WMAP measurements can be found in Ref. 469 . 

2.4 N-Body Simulations 

Our understanding of large scale structure is still far from a satisfactory level. 
The description of the evolution of structures from seed inhomogeneities, i.e. pri- 
mordial density fluctuations, is complicated by the action of many physical pro- 
cesses like gas dynamics, radiative cooling, photoionization, recombination and 
radiative transfer. Furthermore, any theoretical prediction has to be compared 
with the observed luminous Universe, i.e. with regions where dissipative effects 
are of crucial importance. 
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a P J R (kpc) 



Kra 
NFW 
Moore 

Iso 



2.0 3.0 0.4 10.0 

1.0 3.0 1.0 20.0 

1.5 3.0 1.5 28.0 

2.0 2.0 3.5 



Table 2: Parameters of some widely used profile models for the dark matter density 
in galaxies (See Eg. I46|l . Values of R can vary from system to system. 

The most widely adopted approach to the problem of large-scale structure 
formation involves the use of N-body simulations. The first simulation of inter- 
acting galaxies was performed by means of an analog optical computer (Holm- 
berg 1941 j3()l| ) using the flux from 37 light-bulbs, with photo-cells and gal- 
vanometers to measure and display the inverse square law of gravitational force. 
Modern, high resolution simulations make full use of the tremendous increase 
in computational power over the last few decades. 

The evolution of structure is often approximated with non-linear gravita- 
tional clustering from specified initial conditions of dark matter particles and 
can be refined by introducing the effects of gas dynamics, chemistry, radiative 
transfer and other astrophysical processes. The reliability of an N-body sim- 
ulation is measured by its mass and length resolution. The mass resolution is 
specified by the mass of the smallest ("elementary") particle considered, be- 
ing the scale below which fluctuations become negligible. Length resolution 
is limited by the so-called softening scale, introduced to avoid infinities in the 
gravitational force when elementary particles collide. 

Recent N-body simulations suggest the existence of a universal dark matter 
profile, with the same shape for all masses, epochs and input power spectra 
| 393| . The usual parametrisation for a dark matter halo density is 



Various groups have ended up with different results for the spectral shape in the 
innermost regions of galaxies and galaxy clusters. In particular, several groups 
have failed to reproduce the initial results of Navarro, Frcnk and White ' 393| . 
which find a value for the power-law index in the innermost part of galactic 
halos of 7 = 1. In table 2, we give the values of the parameters {a,f3,j) for 
some of the most widely used profile models, namely the Kravtsov et al. (Kra, 
Navarro, Frenk and White (NFW, ^^), Moore et al. (Moore, [HMj ) 
and modified isothermal (Iso, e.g. Ref. jHOl) profiles. 



Although it is definitely clear that the slope of the density profile should 
increase as one moves from the center of a galaxy to the outer regions, the precise 
value of the power-law index in the innermost galactic regions is still under 



p{r) 



(r/i?)-r[l + (r/i?)"](/'-'r)/. 



a 



(46) 
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debate. Attention should be paid when comparing the results of different groups, 
as they are often based on a single simulation, sometimes at very different length 
scales. 

Taylor and Navarro [¥681 ISM] studied the behaviour of the phase-space den- 
sity (defined as the ratio of spatial density to velocity dispersion cubed, p/cr^) 
as a function of the radius, finding excellent agreement with a power-law ex- 
tending over several decades in radius, and also with the the self-similar solution 
derived by Bertschinger |96| for secondary infall onto a spherical perturbation. 
The final result of their analysis is a "critical" profile, following a NFW profile in 
the outer regions, but with a central slope converging to the value ^^tn — 0.75, 
instead of ^nfw = 1- 

The most recent numerical simulations (see Navarro et al. |395j . Reed et 
al. '425? and Fukushige et al. |242| ') appear to agree on a new paradigm, 
suggesting that density profiles do not converge to any specific power-law at 
small radii. The logarithmic slope of the profile continuously flattens when 
moving toward the galactic center. The slope at the innermost resolved radius 
varies between 1 and 1.5, i.e. between the predictions of the NFW and Moore 
profiles. It is important to keep in mind that predictions made adopting such 
profiles probably overestimate the density near the Galactic center and should 
be used cautiously. 

Recently, Prada et al. j421| have suggested that the effects of adiabatic 
compression on the dark matter profile near the Galactic center could play 
an important role, possibly enhancing the dark matter density by an order of 
magnitude in the inner parsecs of the Milky Way. 

The extrapolations of cuspy profiles at small radii have appeared in the past 
(and still appear to some) to be in disagreement with the flat cores observed in 
astrophysical systems, such as low surface brightness galaxies mentioned earlier. 
Such discrepancies prompted proposals to modify the properties of dark matter 
particles, to make them self-interacting, warm etc. Most of such proposals 
appear to create more problems than they solve and will not be discussed here. 

Today, the situation appear less problematic, in particular after the analysis 
of Hayashi et al. |291j . Our approach, given the uncertainties regarding ob- 
served and simulated halo profiles, will be to consider the central slope of the 
galactic density profile as a free parameter and discuss the prospects of indirect 
detection of dark matter for the different models proposed in literature. 

2.5 The Case of the Milky Way 

Since the Milky Way is prototypical of the galaxies that contribute most to the 
cosmic luminosity density, it is natural to ask how the results discussed in the 
previous section compare with the wide range of observational data available 
for our galaxy. 

One way to probe the nature of matter in our neighborhood is to study 
microlensing events in the direction of the galactic center. In fact, such events 
can only be due to compact objects, acting as lenses of background sources, and 
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it is commonly believed that dark matter is simply too weakly interacting to 
clump on small scales^. 

Binney and Evans (BE) |l(J4j recently showed that the number of observed 
microlensing events implies an amount of baryonic matter within the Solar circle 
greater than about 3.9 x 10^°Mq. Coupling this result with estimates of the 
local dark matter density, they exclude cuspy profiles with power-law index 
7^0.3. 

Nevertheless, Klypin, Zhao and Somerville (KZS) |334| find a good agree- 
ment between NEW profiles (7 = 1) and observational data for our galaxy 
and M31. The main difference between these analyses is the value of the mi- 
crolensing optical depth towards the Galactic center used. Observations of this 
quantity disagree by a factor of ~ 3 and a low value within this range permits 
the presence of a dark matter cusp. Another difference arises from the modeling 
of the galaxy: KZS claim to have taken into account dynamical effects neglected 
by BE and to have a "more realistic" description of the galactic bar. 

An important addition is adiabatic compression of the dark matter by bary- 
onic dissipation. This results in a dark matter density that is enhanced in the 
core by an order of magnitude. This result can be reconciled with modelling of 
the rotation curve if the lower value of the microlensing optical depth found by 
the EROS collaboration is used rather than that of the MACHO collaboration. 
In the latter case, little dark matter is allowed in the central few kpc. The mi- 
crolensing result constrains the stellar contribution to the inner rotation curve, 
and hence to the total allowed density. 

2.5.1 The Galactic center 

The dark matter profile in the inner region of the Milky Way is even more 
uncertain. Observations of the velocity dispersion of high proper motion stars 
suggest the existence of a Super Massive Black Hole (SMBH) lying at the center 
of our galaxy, with a mass, Afgj^gjj w 2.6 x IO^Mq |252p . 

Recently, near-infrared high-resolution imaging and spectroscopic observa- 
tions of individual stars, as close as a few light days from the galactic center, 
were carried out at Keck |251| and ESO/VLT telescopes (see Ref. |445| . for an 
excellent and updated discussion of the stellar dynamics in the galactic center, 
based on the most recent observations at ESO/VLT). The analysis of the orbital 
parameters of such stars suggest that the mass of the SMBH could possibly be 
a factor of two larger with respect to the above cited estimate from the velocity 
dispersion. In Fig. [3 we show a plot of the enclosed mass as a function of the 
galactocentric distance, along with a best-fit curve, which corresponds to a dark 
object with a mass of 2.87 ± 0.15 x IO^Mq. 

^It was noticed by Berezinsky et al. I7til that if microlensing was due to neutralino stars (see 
the definition of "neutrahno" in the chapter on dark matter candidates), i.e. self-gravitating 
systems of dark matter particles, then the gamma-ray radiation originated by annihilations 
in these object would exceed the observed emission. 

''The existence of a SMBH at the center of the galaxy is not surprising. There is, in fact, 
mounting evidence for the existence of 10^—10* Mq black holes in the centers of most galaxies 
with mass amounting to approximately 0.1% of the stellar spheroid (see, e.g. Ref. 13421 '). 
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Figure 7: The mass distribution in the galactic center, as derived by different obser- 
vations, down to a 10~'*pc scale. Lines represent fits under different assumptions, as 
specified by the text in the figure. In particular, the solid line is the overall best fit 
model: a 2.87 ± 0.15 x lO®Af0 central object, plus a stellar cluster distributed with a 
power-law of index 1.8. For more details see Ref. |445| . 



It has long been argued (see e.g. Peebles, Ref. 14121 1 that if a SMBH exists 
at the galactic center, the process of adiabatic accretion of dark matter on 
it would produce a "spike" in the dark matter density profile. Gondolo and 
Silk 261 have recently applied such a process to study the enhancement of the 
the annihilation signal from the galactic center. 

If we consider an initial power-law type profile of index 7, similar to those 
discussed in Sec. 12.41 the corresponding dark matter profile, p'(r'), after this 
accretion process is, following Ref. iMLi ■ 
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3-7' 



PD9[r)[-Y\ (47) 
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where = (9 — 27) /(4 — 7), D ~ 8kpc is the solar distance from the Galac- 
tic center and po — 0.3GeV/c^/cm^ is the density in the solar neighborhood. 
The factors and g~j{r) cannot be determined analytically (for approximate 
expressions and numerical values see Ref. |2til| ). Eq. 07|is only valid in a cen- 
tral region of size Rsp — ajD{M/ pDD^y^^^~'^\ where the central black hole 
dominates the gravitational potential. 

It is easy to understand the basics of adiabatic accretion under the assump- 
tions of circular orbits. Assuming an initial power-law distribution, p cx r"'', 
and a final distribution, p cx r~''=p, the equations of conservation of mass and 
angular momentum can be expressed, respectively, as 

Pirfdri = Pfrjdrf (48) 
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and 

nM,{r) = rfMf(r) « rjMBH, (49) 
which imply, respectively, 

(3-7sp)/(3-7) /KfW 

ri cx ^ (50) 



and 



ri (X r 



1/(4-7) 
/ 



(51) 



The final distribution will thus have a power-law index 



7.P = (52) 

which assumes values in the range of 2.25 to 2.5 as 7 varies in the interval of 
to 2. 

If we take into account the annihilation of dark matter particles, the density 
cannot grow to arbitrarily high values, the maximal density being fixed by the 
value ^ 

Pcore — T j ("O) 
av tBH 

where Ibh ~ 10^° yr is the age of the central black hole. The final profile, 
resulting from the adiabatic accretion of annihilating dark matter on a massive 
black hole is 

Pdm{r) = -TT— (54) 

P [r) + Pcore 

following a power-law for large values of r, and with a flat core of density, Pcorei 
and dimension, 

piRs,)^''^''^ 



Rcore — Rsp I 1 ■ (55) 

\ Pcore / 

We will use these equations when discussing the prospects for indirect detec- 
tion of dark matter in the presence of a spike. We recall, nevertheless, that they 
have been derived under the simplifying assumption that the SMBH formed at 
a position coinciding exactly with the center of the galactic potential well, and 
neglecting all dynamical effects. 

It was shown by Ullio, Zhao and Kamionkowski |48H that if the black hole 
forms from a low-mass seed, then its spiral-in to reach the exact center of 
the galaxy could take a length of time longer than the age of the galaxy. If, 
conversely, the seed black hole is massive, the back-reaction to the spiral-in of 
the black hole leads to the formation of a weak-density cusp, with p oc r^"-^. 
Fig. IHI shows the modification of an NFW dark matter profile due to the off- 
center formation of the seed black hole. The solution found by Gondolo and 
Silk would be indistinguishable from the case of an initial light seed of roughly 
10~^ -^SMBH' starting its growth very near to the galactic center. 

The spike could also be destroyed by hierarchical mergers, as discussed by 
Merritt et al. |381| , but such mergers are unlikely to have occurred in the recent 
history of the Milky Way. What can be stated with considerable confidence is 
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Figure 8: Modification of an NFW profile due to tfie off-center formation of a black- 
hole seed of mass Mgjj, its spiral-in the center of the dark matter system and its 
adiabatic growth to the present-day mass, -^SMBH (called in the figure) . From 

Ref. BEI]. 
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that the Milky Way galaxy underwent one significant merger about 12 billion 
years ago. This resulted in the formation of the bulge, and therefore presumably 
of the SMBH, and of the thick disk. The chemical evidence for a unique merger 
origin in the case of our Milky Way's thick disk is compelling |254[ I5U(J) . as the 
continuity between thin disk, thick disk, and bulge would have been destroyed 
had anything significant happened more recently in the way of a merger (see 
also the discussion of Bertone, Sigl and Silk j92j ) . 

Furthermore, the scattering of dark matter particles by stars in the dense 
stellar cusp observed around the SMBH could substantially lower the dark mat- 
ter density near the Galactic center over 10"'^'' years, due both to kinetic heating, 
and to capture of dark matter particles by the SMBH |38U| . 

The existence of such spikes would produce a dramatic enhancement of the 
annihilation radiation from the galactic center. The implications for indirect 
detection of dark matter particles have been discussed in Refs. |261l 12621 19111921 . 

2.5.2 The local density 

Very important to the prospects for direct and indirect detection is the density 
of dark matter in the region of our solar system. Although this quantity is 
considerably more well known than the density near the galactic center, there 
are still uncertainties associated with the local density, which we will discuss 
here. 

The local density of dark matter is determined by observing the rotation 
curves of the Milky Way. This is somewhat difficult to do from our location 
within the galaxy. Furthermore, rotation curves measure the total mass within 
an orbit, thus the density distributions of the galactic bulge and disk are needed 
to accurately calculate the dark matter profile. 

In addition to the local density, the velocity distribution of dark matter in 
the local region is needed to accurately calculate direct and indirect detection 
rates. This is also best inferred from observed rotation curves. 

Different groups have come to somewhat different conclusions regarding the 
local density and velocity distribution of dark matter. For example, Bahcall et 
al. finds a best-fit value of po = 0.34GeV/cm'^ [5^1, Caldwell and Ostriker find 
po = 0.23GeV/cm3 while Turner calculates po = 0.3 - 0.6 GeV/cm^ |T7^ . 
In figure we show the range of local dark matter densities found to be accept- 
able by Bergstrom, UUio and Buckley [50] for various choices of halo profile and 
galactocentric distance. They find local dark matter densities acceptable in the 
range of about 0.2 - 0.8GeV/cm^. 

The velocity distribution of dark matter is typically described only by its 
average velocity, v =< >i/^~ 270km/s. 

For more discussion on the local dark matter distribution, see section 2.4 of 
Ref. [Si- 
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Figure 9: The range of local dark matter densities acceptable with observations of 
rotation curves for a variety of halo profiles and galactocentric distances. Densities in 
the range of 0.2 — 0.8 GeV/cm^ are shown to be acceptable. From Ref. |80|. 

3 Candidates 

As we have seen in the previous section, the evidence for non-baryonic dark 
matter is compelling at all observed astrophysical scales. It is therefore natural 
to ask what is the dark matter made of? In this section, we present some of the 
candidates discussed in the literature, and focus our attention especially on two 
popular candidates: the supersymmetric neutralino, probably the most widely 
studied candidate, and the B'-^^ particle, the first Kaluza-Klein excitation of 
the B boson in theories with universal extra dimensions. We will also briefly 
discuss "superheavy" candidates, also referred to as wimpzillas. 

3.1 The Non-Baryonic Candidate Zoo 

There is no shortage of candidates for non-baryonic dark matter. In this section 
we briefly describe some of these candidates. 

• Standard Model neutrinos 

Neutrinos have been considered, until recently, excellent dark matter can- 
didates for their "undisputed virtue of being known to exist" j86| . However, 
a simple calculation shows that, if we call the mass of the i-th neutrino, 
their total relic density is predicted to be 
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The best laboratory constraint on neutrino masses comes from tritium 
/3-decay experiments at Troitsk and Mainz 14961 , pointing to the foUowing 
upper Hmit on the neutrino mass 

< 2.05 eV (95% C.L.) , (57) 

while next-generation experiments are expected to reach a sensitivity of 
approximately 0.2 eV (see Ref. |496| and references therein). The above 
upper limit applies to all three mass eigenvalues |73| . since the mass 
differences among them must be very small to explain solar (Am^ « 
7 • 10~^ eV^) and atmospheric (Am^ « 3 • lO"*^ eV^) neutrino anomalies 
(see e.g. Ref. |255p . This implies an upper bound on the total neutrino 
relic density of 

n^h^ <, 0.07 , (58) 

which means that neutrinos are simply not abundant enough to be the 
dominant component of dark matter. A more stringent constraint on 
the neutrino relic density comes from the analysis of CMB anisotropics, 
combined with large-scale structure data, suggesting Q^h'^ < 0.0067 (95% 
confidence limit). For three degenerate neutrino species this implies nii, < 
0.23 eV. If extra neutrino interactions are allowed, e.g., the coupling of 
neutrinos to a light boson, the neutrino mass limits arising from large scale 
structure can be evaded |74|. 

Being relativistic coUisionless particles, neutrinos erase (moving from high 
to low density regions) fluctuations below a scale of ~ 40 Mpc eV), 
called the free- streaming length |122j . This would imply a top-down forma- 
tion history of structure in the Universe, where big structures form first. 
The fact that our galaxy appears to be older than the Local Group |413) . 
and the discrepancy between the predicted late formation of galaxies, at 
redshift z ^ 1, against observations of galaxies around z > 4 |123) . is a 
further argument against neutrinos as a viable dark matter candidate. 

• Sterile neutrinos 

These hypothetical particles are similar to Standard Model neutrinos, but 
without Standard Model weak interactions, apart from mixing. They were 
proposed as dark matter candidates in 1993 by Dodelson and Widrow |191| . 
Stringent cosmological and astrophysical constraints on sterile neutrinos 
come from the analysis of their cosmological abundance and the study of 
their decay products (see Ref. T and references therein). 

Light neutrinos, with masses below a few keV, would be ruled out as 
dark matter candidates. In fact, if the WMAP result for the reionization 
optical depth is correct, then dark matter structures were in place to form 
massive stars prior to redshift z > 20, which is is simply not possible 
if the dark matter particle mass is smaller than ~ 10 keV |5()4| . An 
alternative explanation for the WMAP optical depth is reionization by 
decaying particles, such as sterile neutrinos (see Ref. |287j and references 
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therein). Sterile neutrinos could also be cold dark matter, if there is a 
very small lepton asymmetry, in which case they are produced resonantly 
with a non-thermal spectrum |448| . 

• Axions 

Introduced in an attempt to solve the problem of CP violation in particle 
physics, axions have also often been discussed as a dark matter candidate. 

Laboratory searches, stellar cooling and the dynamics of supernova 1987A 
constrain axions to be very light 0.01 eV). Furthermore, they are ex- 
pected to be extremely weakly interacting with ordinary particles, which 
implies that they were not in thermal equilibrium in the early universe. 

The calculation of the axion relic density is uncertain, and depends on 
the assumptions made regarding the production mechanism. Neverthe- 
less, it is possible to find an acceptable range where axions satisfy all 
present-day constraints and represent a possible dark matter candidate 
(see e.g. Ref. fl^ l. 

• Supersymmetric candidates 

— Neutralinos 

Neutralinos in models of R-parity conserving supersymmetry are by 
far the most widely studied dark matter candidates. We devote 
Sec. l3.2l to their presentation. 

— Sneutrinos 

The superpartners of the Standard Model neutrinos in supersymmet- 
ric models have long been considered as dark matter candidates. It 
has been shown that sneutrinos will have a cosmologically interest- 
ing relic density if their mass is in the range of 550 to 2300 GeV. 
However, the scattering cross section of a sneutrino with nucleons is 
easily calculated and is much larger than the limits found by direct 
dark matter detection experiments | 226| . 

— Gravitinos 

Gravitinos are the superpartners of the graviton in supersymmetric 
models. In some supersymmetric scenarios, gauge mediated super- 
symmetry for example, gravitinos can be the lightest supersymmetric 
particle and be stable. Gravitinos are thus very strongly theoretically 
motivated. With only gravitational interactions, however, gravitinos 
are very difhcult to observe |234| . 

It has been known for some time that long lived gravitinos can pose 
problems for cosmology pHl IHIH 1^ 1^ 1^ . In particular, 
their presence can destroy the abundances of primordial light ele- 
ments in some scenarios '385J 11701 12191 1233| . Gravitinos may also be 
overproduced in the early universe if the temperature of the reheat- 
ing epoch is not sufficiently low |385| . In some scenarios, however, 
these problems can be circumvented |132l 11521 1249j . 
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— Axinos 

Axinos, the superpartner of the axion, were beheved until recently 
to only be capable of acting as a warm, or hot, dark matter can- 
didate |124l 126^ . It has been shown, however, that for quite low 
reheating temperatures, cold axino dark matter may be possible 
|164 , 165. 166 366 . In many ways, axinos and gravitinos share sim- 
ilar phenomenological properties. 

• Light scalar dark matter 

Considering fermionic dark matter candidates with standard Fermi inter- 
actions, Lee and Weinberg concluded that relic density arguments pre- 
clude such a WIMP with a mass less than a few GeV )355| (see also 
Hut 1977 |314j ). If the dark matter is made up other types of particles, 
however, this limit could be evaded. For example, a 1-100 MeV scalar 
candidate has been proposed [TTTI I118j . 

Such a candidate, although somewhat ad hoc from a particle physics per- 
spective, has recently become experimentally motivated. In Ref. |120j . it 
has been suggested that the 511 keV gamma-ray line emission observed by 
the INTEGRAL satellite from the galactic bulge could be the product of 
light dark matter particles annihilating into positrons which then annihi- 
late producing the observed gamma-ray line. To confirm this hypothesis, 
more tests are needed. In particular, a similar signature could be expected 
from dwarf spheroidal galaxies |304| . 

Very recently, light decaying dark matter particles such as axinos with 
R-parity violation ,311] or sterile neutrinos {416; have been suggested as 
the source of the observed 511 keV emission. 

• Dark matter from Little Higgs models 

As an alternative mechanism (to supersymmetry) to stabilize the weak 
scale, the so-called "little Higgs" models have been proposed and de- 
veloped [SHlEnilinilll]- In these models, the Standard Model Higgs is 
a pseudo-Goldstone boson with its mass protected by approximate non- 
linear global symmetries. The divergences to the Higgs mass which remain 
are present only at the two-loop level and, therefore, the weak scale can be 
stabilized in an effective field theory which is valid up to ^ 10 TeV. Re- 
call that in supersymmetry, the divergences to the Higgs mass are exactly 
cancelled at all orders. 

At least two varieties of little Higgs models have been shown to contain 
possible dark matter candidates. One of these classes of models, called 
"theory space" little Higgs models, provide a possibly stable, scalar par- 
ticle which can provide the measured density of dark matter |ll()j . In 
Ref. |105| . the detection prospects for such a candidate were found to be 
not dissimilar to WIMPs predicted in models of supersymmetry or uni- 
versal extra dimensions. 
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Cheng and Low |147| have developed another variety of httle Higgs model, 
motivated by the problem of the hierarchy between the electroweak scale 
and the masses of new particles constrained by electroweak precision mea- 
surements. They solve this problem by introducing a new symmetry at 
the TeV scale which results in the existence of a stable WIMP candidate 
with a ~TeV mass. 

For a potential dark matter candidate from a little Higgs model to be 
stable, we must assume that the discrete symmetry which protects it from 
decay is fundamental and is not broken by the operators in the UV com- 
pletion. 

• Kaluza-Klein states 

Kaluza-Klein excitations of Standard Model fields which appear in models 
of universal extra dimensions have also been discussed a great deal recently 
as a candidate for dark matter. They are discussed in Sec. 13.31 Addi- 
tionally, a dark matter candidate has been proposed in the framework of 
"warped" universal extra-dimensions: an exotic particle with gauge quan- 
tum numbers of a right-handed neutrino, but carrying fractional baryon- 
number |15j. 

• Superheavy dark matter 

Superheavy dark matter particles, sometimes called Wimpzillas, have in- 
teresting phenomenological consequences, including a possible solution to 
the problem of cosmic rays observed above the GZK cutoff. These are 
discussed in Sec. 13.41 

• Q-balls 12221 \M3l, mirror particles ^EM EZl \HM 125^'; CHArged 
Massive Particles (CHAMPs) 118^ . self interacting dark matter '^4-58^ 
\17^ , D-matter 144^ > cryptons ]21^ superweakly interacting dark 
matter ^2'84^h brane world dark matter fT^llj , heavy fourth generation neu- 
trinos ,4^1, ,3211, etc. 

Although some of these candidates or classifications present some intrigu- 
ing features we will not discuss them here. We refer the interested reader 
to the wide literature on the subject, e.g. the reviews of non-baryonic 
candidates by Ellis |2(J8| and Bergstrom |86| . 

We stress that it is by no means assured that the dark matter is made of a 
single particle species. On the contrary, we already know that Standard Model 
neutrinos contribute to dark matter, but cannot account for all of it. Even 
in supersymmetry models for dark matter, N — 2 supersymmetry allows the 
possibility of two stable dark matter relics (see, for example, Ref. '119'). 

In what follows, we will assume that the abundance of our candidates satisfy 
the limits provided by the analysis of the CMB discussed in Sec l2.3l but we stress 
that, although the upper bound is a strict limit, the lower bound can be relaxed, 
assuming that our candidate is a sub-dominant component of dark matter. The 
interested reader will find in Ref. |201| a detailed discussion on the detection 
prospects of a subdominant density component of dark matter. 
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3.2 Supersymmetry 



It would be impossible to review in only a few pages the history and theory of 
Supersymmetry (SUSY). Instead, we prefer here to review the motivations that 
led to its introduction and to briefly present the concepts and the notations that 
we will use in the following chapters. Furthermore, we present a few of the su- 
persymmetric models discussed in the literature (we reserve the word "scenario" 
for a specific choice of parameters in the framework of a given model) and dis- 
cuss the consequences of various assumptions, involved in the process of model- 
building, on SUSY phenomenology. For further discussions of supersymmetry, 
we refer the interested reader to Refs. [TOTl ML JM 1103 HHBl • 

3.2.1 Basics of supersymmetry 

As we saw in Sec ll.3l in the Standard Model of particle physics there is a 
fundamental distinction between bosons and fermions: while bosons are the 
mediators of interactions, fermions are the constituents of matter. It is therefore 
natural to ask whether a symmetry exists which relates them, thus providing a 
sort of "unified" picture of matter and interactions. 

Another way to state the problem is to ask whether a Lie group exists mixing 
internal (Isospin, etc.) and space-time (Lorentz) symmetries |282| . Although 
apparently uncorrelated to the differing behaviour of bosons and fermions, this 
problem led to the study of the same algebraic structures. Early attempts to 
find a broad Lie group including the Poincare and internal symmetry groups 
had to face the limitations imposed by the so-called no-go theorem of Coleman 
and Mandula. Such limitations were finally circumvented with the introduction 
of graded Lie algebras, i.e. algebras involving fcrmionic generators satisfying 
anticommutation relations (see below). 

For those who are not convinced by these symmetry arguments, there are 
other major reasons for interest in supersymmetry. One reason is its role in 
understanding the hierarchy problem. The hierarchy problem is linked to the 
enormous difference between the electroweak and Planck energy scales. This 
problem arises in the radiative corrections to the mass of the Higgs boson. 

All particles get radiative corrections to their mass, but while fermion masses 
increase only logarithmically, scalar masses increase quadratically with energy, 
giving corrections at 1-loop of 



where A is a high-energy cut-off where new physics is expected to play an 
important role. The radiative corrections to the Higgs mass (which is expected 
to be of the order of the electroweak scale Mw ^ 100 GeV) will destroy the 
stability of the electroweak scale if A is higher than ~ TeV, e.g. if A is near the 
Planck mass. 

An appealing, though not the only, solution to this problem is to postulate 
the existence of new particles with similar masses but with spin different by one 




(59) 
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Figure 10: The measurements of the gauge couphng strengths at LEP do not (left) 
evolve to a unified value if there is no supersymmetry but do (right) if supersymmetry 
is included 

half. Then, since the contribution of fcrmion loops to Sm^ have opposite sign 
to the corresponding bosonic loops, at the 1-loop level. Eg. 1591 becomes 

^ (^) (A^ + ml) - (^) (A^ + ml) = (|-) (m| - m|) . (60) 

Furthermore, the supersymmetric algebra insures that (provided \ml — ml\ ^ 1 
TeV) the quadratic divergence to the Higgs mass is cancelled at all orders of 
perturbation theory. The algebra of supersymmetry naturally guarantees the 
existence of new particles, with the required properties, associating to all of the 
particles of the Standard Model superpartners with the same mass and opposite 
spin- type (boson or fermion). 

Another reason for interest in supersymmetric theories comes from the unifi- 
cation of gauge couplings at a scale Mjj ^ 2x 10^® GeV (see Fig. llU|l . Although 
extrapolation of the coupling constants using only Standard Model particles 
fails to unify them to a common value (left frame of Fig. I10|l . by introducing 
supersymmetry at the TeV scale, it was shown [55] that these forces naturally 
unify at a scale Mu ~ 2 x 10^^ GeV (right frame of Fig. EJ. This has been 
taken as a strong hint in favor of a Grand Unified Theory (GUT) which predicts 
gauge coupling unification below the Planck scale. 

The new generators introduced with supersymmetry change fermions into 
bosons and vise versa, i.e. 

Qlfermion) = jboson); Qjboson) — jfermion). (61) 

Because of their fermionic nature, the operators Q must carry spin 1 /2, which 
implies that supersymmetry must be a spacetime symmetry. The question then 
arises of how to extend the Poincare group of spatial translations and Lorentz 
transformations to include this new boson/fermion symmetry. The structure 
of such a group is highly restricted by the Haag-Lopuszanski-Sohnius extension 
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of the Coleman and Mandula theorem cited above. For rcahstic theories, the 
operators, Q, which we choose by convention to be Majorana spinors, must 
satisfy 



{Qa,Pf^} = 



(62) 
(63) 
(64) 



where 



a 



(65) 



and 




(66) 



are the structure constants of the theory. 

Just as Lorentz invariance is manifest in Minkowski space-time, supersym- 
metry is manifest in the so-called superspace formalism, where a superspace is 
defined as a set of coordinates {x, 9, 6}, where x — x^^ are the usual coordinate 
of Minkowski spacetime, and 9, 9 are anti-commuting Weyl spinors. 

A superfield is then a function, $(a;, 9, 9), defined on a superspace; it is com- 
mon to introduce chiral fields representing matter and vector fields representing 
gauge fields. 

3.2.2 Minimal supersymmetric Standard Model 

To continue our brief introduction to SUSY, we consider the minimal supersym- 
metric extension of the Standard Model (MSSM, for Minimal Supersymmetric 
Standard Model). The MSSM is minimal in the sense that it contains the small- 
est possible field content necessary to give rise to all the fields of the Standard 
Model. This can be done as follows: 

• We associate fermionic superpartners to all gauge fields. Gluons, and 
B bosons then get fermionic partners called gluinos {g), winos {W^) and 
hinos {B), respectively. The common name for all partners of gauge fields 
is the gaugino. 

• We associate scalar partners to the fermions, i.e. quarks and leptons get 

scalar partners called squarks and sleptons. 

• We introduce one additional Higgs field (for a total of two Higgs doublets, 
corresponding to five physical Higgs states) and associate one spin 1/2 
Higgsino to each Higgs boson. This is done to give masses to both up and 
down-type quarks upon electroweak symmetry breaking and also preserve 
supersymmetry (therefore, we cannot use the conjugate of the Higgs as is 
done in Standard Model). Introducing another Higgs doublet also makes 
the theory anomaly free. 

The resulting particle content of the theory is shown in tables 3 and 4. 
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Superfield 


SM particles 


Spin 


Superpartners 


Spin 


Q 
D" 


(S) 

dR 


1/2 

1/2 
1/2 


m 

d* 








L 


(S) 


1/2 
1/2 


3^ 






Hi 
Hi 


Hi 
H2 






Hi 
62 


1/2 
1/2 


G" 
W, 
B 


9 

W, 
B 


1 
1 
1 


a 

B 


1/2 
1/2 
1/2 



Table 3: Field content of the MSSM. 



The MSSM is then specified through the superpotential, defined as 

W = [veHiVE' + VdHiQW + VuHlQW] + e,j/iHji/| (67) 

where i and j are SU(2) indices, and y are Yukawa couphngs. Color and gener- 
ation indices have been suppressed in the above expression. The superpotential 
represents a supersymmetrization of the Standard Yukawa couplings plus a bilin- 
ear Higgs term. The superpotential enters the Lagrangian of the theory through 
the terms 

-^^SUSY = -\{W''^^^o + W^* V'^V^'^) - W^W* (68) 

where we have used = W* = dW/d(j)'*, and W'^ = d^W/d^id4>j . 

4>i and Tpi are scalar and fermion fields, respectively. 

One additional ingredient of the MSSM is the conservation of i?-parity. im- 
parity is a multiplicative quantum number defined as 

i?= (-1)3^+^+2^ . (69) 

All of the Standard Model particles have i?-parity R = 1 and all sparticles 
{i.e. superpartners) have R = —1. Thus, as a consequence of i?-parity conserva- 
tion, sparticles can only decay into an odd number of sparticles (plus Standard 
Model particles). The lightest sparticle (dubbed the LSP, for Lightest Super- 
symmetric Particle) is, therefore, stable and can only be destroyed via pair 
annihilation, making it an excellent dark matter candidate 25f| I213| . Note 
that this not the original motivation for R-parity. In fact, R-parity was first 
introduced to suppress the rate of proton decay |228[ 11831 1495L I185| . 

The nature of the LSP in the MSSM is constrained by many observations. 
It cannot have a non-zero electric charge or color, or it would have condensed 
with baryonic matter to produce heavy isotopes, in conflict with observations. 
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Standard Model particles and fields Supersymmetric partners 







Interaction eigenstates 


Mass eij 
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Name 
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q = d,c, b, u, s, t 


quark 


qL, qR 


squark 


<?1. 92 


squark 


I = e, fi,T 


lepton 




slepton 


h, h 


slepton 


V = Ve., V^, Vr 


neutrino 


V 


sneutrino 


V 


sneutrino 


9 


gluon 


a 


gluino 


a 


gluino 




VF-boson 




wino 1 


\ Xt2 




H- 


Higgs boson 




higgsino ' 


chargino 


H+ 


Higgs boson 




higgsino J 






B 


B-field 


B 


bino > 








Vy3-field 




wino 


1 ^1,2,3,4 




HO 


Higgs boson 


m 


higgsino 1 


neutralino 




Higgs boson 




HI 


Higgs boson 


higgsino J 







Table 4: Standard Model particles and their superpartners in the MSSM (adapted 
from Ref. HOSI)- 



Among the neutral candidates, a possibile LSP could be the sneutrino. Sneu- 
trino LSPs have, however, been excluded by direct dark matter detection ex- 
periments (see sections lOl and 1^ . Although axinos and gravitinos cannot be 
a prori excluded, they arise only in a subset of supersymmetric scenarios and 
have some unattractive properties (see section I5!T|) . In particular, gravitinos 
and axinos have very weak interactions and would be practically impossible to 
detect, making them less interesting from a phenomenological perspective. The 
lightest neutralino remains an excellent dark matter candidate, and is further 
discussed in the next section. 

To determine the identity of the LSP (or other characteristics) in a given 
supersymmetric scenario, we have to specify how supersymmetry is broken. If 
supersymmetry were not broken, then each superpartner would have a mass 
identical to its Standard Model counterpart, which is clearly not the case. Thus, 
new terms which break supersymmetry must be added to the Lagrangian. These 
terms, however, should be added carefully, in order not to destroy the hierarchy 
between Planck and electroweak scales. The possible forms for such terms are 

-^(BMY'Mj - ^{Ayy^Mj^k + h.c, (70) 

where the are gaugino masses, are soft scalar masses, B is a bilinear 
mass term, and A is a trilinear mass term. We will discuss some specific super- 
symmetry breaking scenarios later in this section. 
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3.2.3 The lightest neutralino 

In the MSSM, the superpartners of the B, W3 gauge bosons (or the photon and 
Z, equivalently) and the neutral Higgs bosons, and H2, are called binos 
{B), winos (1^3), and higgsinos (iJ^ and H2), respectively. These states mix 
into four Majorana fermionic mass eigenstates, called neutralinos. The four 
neutralino mass eigenstates are typically labelled Xi; X2j X3 ^.nd X41 ordered 
with increasing mass. In the following we will refer to Xij the lightest of the 
four neutralinos, as the neutralino, and denote it simply as, x = Xi- 

In the basis (i3, IV3, if", i/^ )7 the neutralino mass matrix can be expressed 

as 

I Ml -Mz cos 13 sinOyy Mz sin (3 sinOyy \ 

M2 MzCOsPcosOy^ —Mz sin (3 cosO^f 

-Mz cos (3 sinOyy Mz cos [3 cosO^, -jjL 

\ Mz sin (3 sin 9-^ —Mz sin [3 cosOyy -fi J 

(71) 

where Mi and M2 are the bino and wino mass parameters, respectively, 9w 
is the Weinberg angle and tan /3 is the ratio of the vacuum expectation values 
of the Higgs bosons, fi is the higgsino mass parameter. As we have seen, the 
(lightest) neutralino is a linear combination of B, W3, Hi and H2, 

X - NiiB + N12W3 + Nr:iH° + N^Hl (72) 

We then define the gaugino fraction, fc, and the higgsino fraction, fn, as 

/g = Nfi + Nf2 (73) 

and 

fH = iV^g + Ni^ . (74) 

For the analytic expressions used to diagonalize the neutralino mass matrix, see 
Appendix 1X1 

The neutralino interactions most relevant for the purposes of dark matter are 
self annihilation and elastic scattering off of nucleons. Neutralinos are expected 
to be extremely non-relativistic in the present epoch, allowing us to safely keep 
only the a-term in the usual expansion of the annihilation cross section, 

av^a + bv^ + {v^) . (75) 

The 5-term must be included in performing calculations of the neutralino relic 
density, however. 

At low velocities, the leading channels for neutralino annihilation are anni- 
hilations to fermion-antifermion pairs (primarily heavy fermions, such as top, 
bottom and charm quarks and tau leptons), gauge bosons pairs (M^"*"!^" and 
Z^Z^) and final states containing Higgs bosons. In appendix IbI we give the 
most important neutralino annihilation diagrams, amplitudes and cross sec- 
tions (in the low velocity limit). For a complete list of all tree level processes, 
diagrams, amplitudes and cross sections, see the excellent review of Jungman, 
Kamionkowski and Griest |319j . 
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3.2.4 Supersymmetric models 

Although relatively simple in many respects, the MSSM has a huge number of 
free parameters. Most of these parameters represent masses and mixing an- 
gles, much as in the case of the Standard Model. To allow for the practical 
phenomenological study of the MSSM, the number of parameters which are 
considered must be reduced. This can be done by making (theoretically well 
motivated) assumptions which reduce the free parameters from more than 100 to 
a more tractable quantity. Depending on the assumptions used, one obtains dif- 
ferent supersymmetric models. In the following section, we will describe a few of 
the most widely considered supersymmetric scenarios, including mSUGRA (of- 
ten called the constrained MSSM) and a phenomenologically simplified MSSM 
(called the phenomenological, or, pMSSM). We also discuss the phenomenolog- 
ical features of the MSSM in anomaly, gauge and gaugino mediated scenarios. 

mSUGRA The mSUGRA, or constrained MSSM, scenario is a simple phe- 
nomenological model based on a series of theoretical assumptions (see e.g. Kane 
et al. \327\ ). The number of free parameters is reduced in this scenario by 
assuming that the MSSM parameters obey a set of boundary conditions at the 
Grand Unification scale: 

• Gauge coupling unification 

ai{Mu) = a2{Mu) ^ asiMu) = au (76) 

with a.i = /47r 

• Unification of the gaugino masses 

Ml {Mu)= M2 {Mu ) = Ms {Mu ) = mi/2 (77) 

• Universal scalar [sfermion and Higgs boson] masses 

MQiMu) = A'h,{Mu) = M^jjMu) = Ah{Mu) = M^JM^) 

= Mff„ {Mu) - Mh, {Mu) = mo (78) 

• Universal trilinear couplings: 

Au{Mu) ^ AdiMu) ^ Ai{Mu) = Ao (79) 

By requiring the minimization of the Higgs potential (in order to recover 
electroweak symmetry breaking), we are left with five (four continuous and one 
discrete) free parameters: 

tan/3 , mi/2 , mo , ^0 , sign(/i), (80) 

where tan (3 is the ratio of the vacuum expectation values of the two Higgs fields 
and fi is the higgsino mass parameter. 
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tanB = 10,u>0 tanP = 10,n<0 
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Figure 11: The (mi/2,n^o) planes for mSUGRA with (a) tan/3 = 10, p > 0, (b) 
tan /3 = 10, /i < 0, (c) tan /3 = 35, < 0, and (d) tan /9 = 50, /i > 0. In each panel, the 
region allowed by the older cosmological constraint 0.1 < H^ft^ < 0.3 has cyan shading, 
and the region allowed by the newer cosmological constraint 0.094 < fl^h^ < 0.129 
has dark blue shading. For more details, see Ref. |217| . 

A recent study of mSUGRA parameter space in light of the WMAP mea- 
surement of the dark matter relic density can be found in Ref. [21 7j . We show in 
Fig. II Hand Fig.^jthe regions of the (rni/2, mo) plane consistent with CMB and 
accelerator data. It is worth mentioning that neutralino models with relic den- 
sities lower than the WMAP measurement are not ruled out, although evidently 
they cannot make up all the dark matter. 

In addition to constraints on models in mSUGRA which come from the 
WMAP measurements, strong constraints can also be placed by collider data. 
In particular, constraints arise from the absence of new particles at LEP below 
« 100 GeV and the agreement of 5 — > 57 decays with predictions of the Standard 
Model. Measurements of the anomalous magnetic momentum of the muon, — 
2, also provide a possible constraint. These constraints have been studied in the 
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Figure 12: Regions of the (mi/2, Jtiq) plane in mSUGRA that are compatible 
with 0.094 < Q-^h^ < 0.129 and laboratory constraints for fi > and tan/3 = 
5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55. The parts of the strips compatible with - 2 
at the 2-a level have darker shading. From Ref. |217| . 



context of mSUGRA in great detail gTlEIHlllSni- The interested reader will find 
a discussion of mSUGRA parameters and the definition of SUSY benchmarks 
points in Ref. 65 . For more on collider constraints, see section [3. 51 



The phenomenological MSSM The scenario we present in this section is 
not necessarily motivated by any theoretical arguments, but rather is justified 
by focusing on the aspects of supersymmetric phenomenology which are the 
most interesting for neutralino dark matter. The phenomenological MSSM, or 
pMSSM, is an adaptable framework which can be described by as many as tens 
of parameters, or as few as five or seven. It is NOT a model, but rather a 
convenient description of the phenomenology most relevant to supersymmetric 
dark matter. Common choices in defining a phenomenological MSSM include 
a) no new sources of CP violation (all the phases in the SUSY breaking terms 
are set to zero), b) no fiavor-changing neutral currents and c) first and second 
generation universality. 

One example of a phenomenological MSSM is used in the DarkSusy program 
package |264| . In this scheme, in addition to the common features described 
above, gaugino unification is assumed (similar to Ea.l77|l. The remaining inputs 
are defined by seven free parameters: 

/i, M2,tan/3, Ma, mo, Ah and At (81) 

where Ma is the mass of the pseudo-scalar Higgs boson, tuq is the common 
scalar mass, and Ab^t are trilinear couplings appearing in SUSY breaking terms. 
Unlike in the case of the mSUGRA scenario, the input parameters are chosen at 
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the electroweak scale without making use of renormahzation group equations. 
The inputs used in DarkSusy can be expanded beyond these seven to include 
other parameters, thus representing a more general MSSM. 

We note that these scenarios are less theoretically motivated in comparison 
to mSUGRA. Various combinations of theoretically and phenomenologically- 
based descriptions for supersymmetry are often considered in the literature, 
often maintaining some of the theoretically motivated constraints of mSUGRA 
while relaxing other requirement (for example, see Refs. |9(JI llU9j ). 

The focus point region of mSUGRA In most of the parameter space 
of mSUGRA or other similar scenarios, the lightest neutralino is a gaugino- 
likc neutralino with a mass of a couple hundred GeV or less. In the so-called 
"focus point" region of mSUGRA, however, the lightest neutralino can have a 
considerable higgsino content, and be significantly more heavy |194[ 11981 1229L 
1230, 231 . 

In the focus point region, very large scalar masses are possible without vi- 
olating naturalness constraints. This occurs because the soft masses squared 
of the Higgs boson, mj^ , have pseudo fixed-point behavior, and can start with 
a wide range of input values and run to a similar negative value at the low 
scale. This is interesting because it indicates that, in the focus point region, 
electroweak symmetry breaking does not require fine-tuning in the high energy 
input values. 

A typical feature of the focus point region are large scalar masses (usually ~ 
TeV). The main reason for a larger higgsino content in the LSP is the larger input 
value of the soft scalar mass. The tree level electroweak symmetry breaking 
condition gives 

^ml^-mj,^^^ (82) 

In the typical mSUGRA scenarios, rnj^ is driven to some large negative value 
due to the running of the renormahzation group equations. This requires a 
large value of fi to give the correct Z mass. In the focus point region, however, 
it is possible that the large input value of the scalar soft mass makes rnj^ less 
negative. Hence, a smaller value of is possible, which leads to a larger higgsino 
content in the LSP. 

Anomaly mediated SUSY breaking Anomaly Mediated Supersymmetry 
Breaking (AMSB) is an attractive alternative to general gravity mediated sce- 
narios as it provides an elegant solution to the so-called flavor problem through 
an elegant decoupling mechanism. The resulting soft parameters are "UV insen- 
sitive" . In this scenario, the SUSY parameters can all be expressed in terms of 
low energy parameters such as the Yukawa and gauge couplings. Although the 
details of AMSB are quite technical, and are beyond the scope of this work, we 
will here describe some of the phenomenological features of this scenario which 
are most relevant for dark matter. 

In AMSB, the gaugino spectrum is given by 
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(83) 



where P are beta-functions computed in the supersymmetric limit |424l I255L 
I250j . The proportionaUty of the /3-function to the low energy masses leads to a 
specific relationship between the gaugino masses: Mi : M2 : M3 = 2.8 : 1 : 7.1. 
This is very different than would be predicted by the GUT relations used in 
mSUGRA, for example, with the wino and bino mass hierarchy reversed. When 
the neutralino mass matrix is diagonalized in AMSB, these ratios result in an 
LSP which is almost purely (neutral) wino EHZIIlTni. Additionally, a charged 
wino, with a mass only a few hundred MeV heavier than the LSP, is predicted. 
This leads to a long lived chargino with distinctive collider signatures. Also, in 
AMSB, with such large values of AI3, the gluino and squarks are predicted to 
be considerably heavier than in mSUGRA or other similar scenarios. 



The heterotic orbifold model The weakly coupled heterotic string with 
orbifold compactification is among the earliest and best understood string mod- 
els that can accomodate in four dimensions the Standard Model gauge group, 
three generations of squarks and a coherent mechanism of supersymmetry break- 
ing. These models show a behavior that interpolates between the phenomenol- 
ogy of mSUGRA and models dominated by superconformal anomalies (AMSB) 

Recently, the full one loop soft supersymmetry breaking terms in a large class 
of superstring effective theories have been calculated (10"2| based on orbifold 
compactifications of the weakly coupled heterotic string (including the so-called 
anomaly mediated contributions). The parameter space in this class of models 
has already been severely constrained by taking into account accelerator and 
relic density constraints, as well as direct and indirect searches (see Refs. |95l 

nnniEniESiiinHi- 



Gauge mediated SUSY breaking Another alternative SUSY breaking mech- 
anism is mediated by gauge interactions 186, 391, 28, 187, 188", 189\ In Gauge 
Mediated Supersymmetry Breaking (GMSB), we have the following approxi- 
mate relationship between the low energy SUSY masses and the gravitino mass 

"^^/^ - ' «1, (84) 



WSUSY Ola Mpi 

where Ms is some typical supersymmetry breaking scale. Therefore, generically, 
we will have a very light gravitino as the LSP |184| . Such a scenario provides 
a dark matter candidate which is very difficult to observe. We will not discuss 
gravitino dark matter or GMSB further for this reason. 



Gaugino mediated SUSY breaking Gaugino mediated supersymmetry 
breaking |330[ I146j represents another class of SUSY breaking mediation moti- 
vated by the brane-world scenario. In this scenario, unwanted supersymmetry 
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breaking effects, such as the flavor violating couplings, are suppressed by the 
separation of the observable and hidden sectors via the separation of their re- 
spective branes. Gauginos are allowed to propagate off of the branes (in the 
bulk) in this scenario, communicating SUSY breaking from the hidden sector. 

The most important phenomenological feature of this mechanism is that the 
sfermion masses are suppressed relative to the gaugino masses. This is because 
sfermion masses can only be generated from the 1-loop diagrams in which a 
gaugino is emitted, travels through the bulk to the supersymmetry breaking 
brane, gets the information of SUSY breaking and then returns to join the 
sfermion propagator again. Generically, these masses are suppressed relative to 
the gaugino mass by a loop factor, TOj ~ M|/(167r^). 

3.3 Extra Dimensions 

Although our world appears to consist of 3+1 (three space and one time) di- 
mensions, it is possible that other dimensions exist and appear at higher energy 
scales. 

From the physics point-of-view, the concept of extra dimensions received 
great attention after the idea of Kaluza, in 1921, to unify electromagnetism with 
gravity by identifying the extra components of the metric tensor with the usual 
gauge fields. More recently, it has been realized that the hierarchy problem 
(see section r}.2.1f) could be addressed, and possibly solved, by exploiting the 
geometry of spacetime. 

In many extra-dimensional models, the 3+1 dimensional spacetime we expe- 
rience is a structure called a brane, which is embedded in a (3 + (5+ 1) spacetime 
called the bulk. The hierarchy problem can then addressed by postulating that 
all of the extra dimensions are compactified on circles (or other topology) of 
some size, R, as has been done in the Arkani-Hamed, Dimopoulos and Dvali 
(ADD) scenario 02]) thus lowering the fundamental Planck scale to an energy 
near the electroweak scale. Alternatively, this could be accomplished by intro- 
ducing extra dimensions with large curvature (warped extra dimensions) as has 
been suggested by Randall and Sundrum j423| . The extra dimensional scenario 
which we will focus on throughout the remainder of this review (universal extra 
dimensions) does not share the features of the ADD or RS scenarios. Rather, it 
introduces flat extra dimensions which are much smaller than those in the ADD 
framework. 

In addition to the hierarchy problem, motivation for the study of theories 
with extra dimensions comes from string theory and M-theory, which today 
appear to be the best candidates for a consistent theory of quantum gravity 
and a unifled description of all interactions. It appears that such theories may 
require the presence of six or seven extra-dimensions. 

A general feature of extra-dimensional theories is that upon compactifica- 
tion of the extra dimensions, all of the fields propagating in the bulk have their 
momentum quantized in units of ~ The result is that for each bulk 

field, a set of Fourier expanded modes, called Kaluza-Klein (KK) states, ap- 
pears. From our point of view in the four dimensional world, these KK states 
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appear as a series (called a tower) of states with masses m„ = n/R, where n 
labels the mode number. Each of these new states contains the same quantum 
numbers, such as charge, color, etc. 

In many scenarios, the Standard Model fields are assumed to be confined 
on the brane, with only gravity allowed to propagate in the bulk. Nevertheless, 
if the extra-dimensions are small, it would be possible for all fields to freely 
propagate in the extra dimensions. Such is the case in models with universal 
extra dimensions, which we discuss in the next section. 

3.3.1 Universal extra dimensions 

Scenarios in which all fields are allowed to propagate in the bulk are called 
Universal Extra Dimensions (UED) Following Ref. |446) . we note that 

there is significant phenomenological motivation to have all Standard Model 
fields propagate in the bulk, including: 

• Motivation for three families from anomaly cancellation 

• Attractive dynamical electroweak symmetry breaking 

• Prevention of rapid proton decay 

• Provides a viable dark matter candidate. 

In the case of one extra dimension, the constraint on the compactification 
scale in UED models from precision electroweak measurements is as low as 
^ 300 GeV [33]. Recently, it was shown that this bound can be weakened 
to ^ 280 GeV if one allows a Higgs mass as heavy as mn k, 800 GeV 
|36) . This is to be contrasted with another class of models where Standard 
Model bosons propagate in extra dimensions while fermions are localized in 4 
dimensions. In such cases, the constraint on the compactification scale is much 
stronger, requiring R^^ ^ several TeV |151| . 

The prospect of UED models providing a viable dark matter candidate is 
indeed what motivates us in our discussion here. The existence of a viable dark 
matter candidate can be seen as a consequence of the conservation of momen- 
tum in higher dimensional space. Momentum conservation in the compactified 
dimensions leads to the conservation of KK number. This does not stabilise the 
lightest KK state, however. To generate chiral fermions at the zero mode, the 
extra dimensions must be moded out by an orbifold, such as S/Z2 for one extra 
dimension or T"^ /Z2 for two. This orbifolding results in the violating of KK 
number, but can leave a remnant of this symmetry called KK-parity (assuming 
that the boundary terms match). All odd- level KK particles are charged under 
this symmetry, thus ensuring that the lightest (first level) KK state is stable. In 
this way, the Lightest Kaluza-Klein Particle (LKP) is stabalized in a way quite 
analogous to the LSP in R-parity conserving supersymmetry. 

In the next section, we will discuss some of the characteristics of the LKP 
in models of UED. 
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3.3.2 The lightest Kaluza Klein particle 

The study of the Lightest Kaluza-Klein Particle (LKP) as a dark matter candi- 
date dates back to the work of Kolb and Slansky in 1984 |338| . where the KK 
excitations were referred to as pyrgons, from the Greek Trvpjo<; for "scale" or 
"ladder" . The LKP has since been reconsidered in the framework of universal 
extra dimensions, in which it is likely to be associated with the first KK excita- 
tion of the photon, or more precisely the first KK excitation of the hypercharge 
gauge boson |149| . We will refer to this state as B^^\ 

A calculation of the B^^'^ relic density was performed by Servant and Tait |44f)j . 
who found that if the LKP is to account for the observed quantity of dark mat- 
ter, its mass (which is inversely proportional to the compactification radius R) 
should lie in the range of 400 to 1200 GeV, well above any current experimental 
constraint. 

We show in Fig. ^| the relic density of the -B*^^-* particle versus its mass, 
including coannihilations (see section [T3l with the next-to-lightest KK particle, 
which in the case shown is e^^\ the first KK excitation of the right-handed 
electron. This figure is a new version of Fig. 3 in Ref. |446| , updated to include 
the new WMAP constraints on the cold dark matter relic density 

Note that the results of the LKP relic density calculation can vary depending 
on the spectrum of other first level KK states. Unlike in the case of supersym- 
metry, the density of KK dark matter is increased through coannihilations with 
other KK particles. 

This is due to the fact that in the case of neutralinos, the cross section 
for the interaction between neutralinos and the NLSP is much larger than the 
neutralino self-annihilation cross section, which implies that DM particles are 
kept longer in thermodynamic equlibrium, thus decoupling with a lower relic 
density. In constrast, the interactions between the i?*^^^ and e^'' are comparable 
with the B^^^ self-interaction. Decoupling in presence of coannihilations thus 
happens essentially at the same time as in the case with no coannihilations, and 
the B^^^ relic density becomes larger since the after decoupling at the same 
time, decays in the B^-^\ 

The spectrum of first level KK states has been calculated to one loop by 
Cheng, Matchev and Schmaltz |149| . although higher dimensional operators 
localized at the boundary may change the details of the spectrum (without 
affecting KK parity). Variations in this spectrum can result in variations for 
the predicted LKP relic abundance. 

The B^^^ annihilation cross section has been studied in Ref. |446| . and is 
given by 

95gt 0.6 pb 

(TV — ^ ~ ^ ; . (85) 

3247rm^,,, "i^(i,[TeV] ^ ' 

The branching ratios for B^^^ annihilation (see table 5) are almost independent 
of the particle mass. Unlike in the case of supersymmetry, the bosonic nature of 
the LKP means that there will be no chirality suppression in its annihilations, 

*see section [2. 31 for a discussion of the CMB and, in particular, the recent WMAP data 
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Channel 


Branching ratio 


quark pairs 


35% 


charged lepton pairs 


59% 


neutrino pairs 


4% 


Higgs bosons 


2% 



Table 5: Branching ratios for the annihilation of the B'^' particle. Note that small 
variations from these results can occur with variation in the KK spectrum. 




Figure 13: Relic density versus mass of tlie B^^\ The solid line is the case for B^^^ alone, 
dashed and dotted lines are for one (three) flavors of nearly degenerate e^'. For each case, 
black curves (upper of each pair) are for A = 0.01 and red curves (lower of each pair) for 
A = 0.05. Figure kindly provided by G. Servant. 



and thus can annihilate efficiently to fermion-fermion pairs. In particular, since 
the annihilation cross section is proportional to hypercharge'* of the final state, 
a large fraction of LKP annihilations produce charged lepton pairs. 

Direct detection of the LKP via its elastic scattering with nuclei was investi- 
gated in Refs. , 150. , 447^ . It was emphasized in Ref. 447 that a one-ton detector 
is needed to probe the expected heavy masses as indicated by the relic density 
calculation |446j of the LKP. One must, therefore, wait for the next generation 
of direct detection experiments such as GENIUS ji33 or XENON 37 (see sec- 
tion Simultaneously, the LHC should probe most of the relevant KK mass 
parameter space (up to R^^ ~ 1.5 TeV [148) ') and confirm or rule out UED at 
the TeV scale. 
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3.4 Superheavy Candidates 

Dark matter particles are usually assumed to be relatively "light", meaning 
lighter than a few hundred TeV. This "limit" is a consequence of the existence 
of a maximum annihilation cross section, av, for a particle of a given mass, 
niDM, set by the so-called unitarity hound (see e.g. Ref. 14931 1. Griest and 
Kamionkowski |277| applied this bound and the constraint on the relic density 
to infer an upper limit on the dark matter particle mass: 

moM^ 340 TeV . (86) 

We note that nowadays, using the WMAP constraint on fluMh^, such a con- 
straint can be made ten times stronger, 

moM ^ 34 TeV . (87) 

The assumption behind this argument is that the dark matter particle is a 
thermal relic of the early Universe, otherwise we could not have applied the 
relation between fljyMh^ and av. In this section, we consider superheavy dark 
matter candidates, defined as candidates with mass raoM > 10^° GeV, that we 
call generically wimpzillas |339l I144| . Thus the first condition for this scenario 
is that wimpzillas must not have been in thermal equilibrium during freeze- out. 
Since they are not in thermal equlibrium during freeze-out, their relic abundance 
does not depend on their annihilation cross section, but rather is a function of 
the wimpzilla's production cross section. Furthermore, we want them to be 
sufficiently stable against decay and annihilation to significantly contribute to 
the present day matter density. 

There are many ways to produce wimpzillas in the early Universe. Among 
the most studied is gravitational production at the end of inflation, resulting 
from the expansion of the background spacetime (for details on this and other 
scenarios see e.g. Ref. |339| and references therein). 

Natural mass scales for wimpzillas include the inflaton or grand unified 
masses, which are usually assumed to be roughly 10^^ GeV or 10^^ GeV, re- 
spectively. Alternatively, D-matter provides a good candidate for wimpzillas 
with a somewhat larger mass 449 . The interaction cross sections with ordi- 
nary matter for such particles can vary from very weak to strong (in the latter 
case supermassive particles are sometimes called simpzillas). 

A common motivation for superheavy dark matter comes from the observa- 
tion of cosmic rays at ultra- high energies (l 6). above the so-called GZK (Greisen- 
Zatsepin-Kuzmin) cutoff [2731 I5()6| . Above this cutoff, which occurs at ~ 
5 X 10"'^^ eV, protons interact at resonance with CMB photons with a center- 
of-mass energy nearly equal to the mass of the A-hadron (1.232 GeV). The 
cross section for this interaction is quite large, thus making the Universe opaque 
to ultra- high energy protons over cosmological distances ( > 50 Mega-parsecs) . 
Since no astrophysical sources of ultra-high energy protons are known within this 
range, more exotic scenarios have been developed to account for these observed 
events. Such scenarios include ultra-high energy cosmic-ray production via 
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the decay or annihilation of superheavy dark matter particles, called top-down 
cosmic-ray models (see, for example, Refs. fTT^ 17^ BI^I^ITTTllKnHl 111)11 1^ ). 



3.5 Collider Constraints 

The constraints which can be placed on a dark matter candidate from collider 
experiments are highly model dependent in nature. It is, unfortunately, impos- 
sible to completely or simply describe the reach of colliders in their search for 
dark matter in any kind of general way. We will here, rather, review several 
of the most important collider searches which have been carried out for dark 
matter particles and for particles associated with a dark matter candidate. 

3.5.1 Current collider constraints 

• Invisible Z width 

If a dark matter candidate is sufficiently light, Z bosons may decay in- 
visibly to such particles with a non-zero branching fraction. Of course, 
there is a substantial background to such events, namely Z ^ vv decays. 
Presently, to contribute less than one standard deviation to the measured 
neutrino contribution, the analysis of LEP2 finds that a decay width of 
^ z~>xx < 4.2 MeV is required {X denotes a dark matter particle). 

Similarly, single photon events can be an interesting search channel for 
light dark matter particles. At LEP2, the Standard Model background 
process for this signature is e^e^ Z vi? with an additional photon 
radiated off of the initial state. The total cross section for 7Z production 
at LEP2 is less than 31 pb with a minimum 1 GeV transverse momentum 
cut for the photon. The contribution to this final state (single photon) 
from particle dark matter, in addition to Z ^ XX with an additional 
photon radiated off of the initial state, are t-channel dark matter produc- 
ing processes in which a photon is radiated off of a charged propagator (a 
selectron in supersymmetry, for example). 

• Searches for new charged particles 

LEP2 has placed very stringent bounds on charged particles lighter than 
about 100 GeV. In e+e" colliders, cross sections for the direct pair pro- 
duction of charged particles are quite large, allowing for limits to be placed 
at or slightly below half of the center-of-mass energy of the collision. For 
LEP2, which reached a center-of-mass energy of 209 GeV, limits of 87- 
103 GeV have been placed for such particles; in particular for charginos 
(m^± > 103 GeV) and charged sleptons {mg > 99 GeV, to^ > 96 GeV, 
irif > 87 GeV) in supersymmetric models. If the LSP is only slightly less 
massive than the charged particle, however, this limit may be substantially 
lower PII^ITTII^. 

Limits on charged particles can only indirectly constrain dark matter, 
however. In supersymmetry, chargino masses and neutralino masses are. 
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in some models, related by the unification of gaugino masses. Although 
such a relationship is often assumed, it is quite possible that the pattern 
of gaugino masses is not so simple. If gaugino mass unification is assumed 
chargino limits can translate to neutralino mass limits of about half of the 
chargino mass limit {m^o > 50 GeV). Without such a relationship, the 
LSP could be much lighter [TOIT^ . 

Sneutrino limits 

Limits for charged sleptons can be used to indirectly limit the possible 
masses for sneutrinos beyond the invisible Z width constraints. Such a 
bound is the result of a basic SU(2) symmetry between the supersymmetry 
breaking masses of the left handed slepton and the sneutrino of a given 
lepton flavor. Limits somewhat lower than for charged particles {mci ^ 85 
GeV) can be placed on sneutrinos if such theoretical assumptions are made 

oiniEnsi 

Searches for colored particles 

Hadron colliders, such as the Tevatron, can place the strongest limits on 
colored particles (squarks and gluinos in supersymmetry or KK excita- 
tions of quarks and gluons in models with universal extra dimensions, for 
example). Such particles would most likely undergo a series of cascades 
upon decaying, possibly producing dark matter candidates among other 
particles. Combinations of squarks and gluinos are searched for using jets 
and missing energy signatures. This leads to exclusion contours in the 
squark and gluino mass plane. 

In supersymmetry, the spectrum of neutralinos, charginos and sleptons 
lighter than the decaying squark/gluino is very important in placing limits 
on squark and gluino masses. Similar ambiguities are present in other 
models as well, such as universal extra dimensions, etc. Typically, limits 
of ~200 GeV are obtained for new colored particles, unless there exists an 
invisible final state particle with a mass close to the new colored particle's 
mass HElliElllini. 

New gauge bosons 

Heavy gauge bosons appear in many models of particle physics beyond the 
Standard Model. Heavy charged gauge bosons (called W 's) and heavy 
neutral gauge bosons (called Z 's) have been excluded below about 600- 
800 GeV, depending on the details of the analysis |14[ |5] . These limits 
assume that these particles have couplings equal to their Standard Model 
counterparts. If their couplings were smaller, the resulting limits could 
be considerably weaker. Electroweak precision measurements can also 
constrain heavy gauge bosons considerably (see below). 

Higgs searches 
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In supersymmetric models, the Higgs mass is increased from its tree level 
mass (below mz) by loop processes involving superparticles, most im- 
portantly top squarks. Current bounds on the (lightest) Higgs mass, 
therefore, constrain the masses of top squarks and other superparticles. 
Furthermore, if supersymmetry is manifest below 1 TeV, as is normally 
expected, the Higgs mass must be less than about 130 GeV, not very far 
above current limits from LEP2 {mn < 114.1 GeV) |^ 1^ 1^ 1^ 
121 El • Note that this bound is somewhat lower for cases with very large 
values of tan/3. 

Searches for charged Higgs bosons can also provide constraints on models 
of physics beyond the Standard Model. More sensitive, however, may be 
the impact of charged Higgs bosons in the branching fraction of 6 ^ 57. 

Flavor changing neutral currents 

Many models of physics beyond the Standard Model introduce flavor 
changing neutral currents, often at tree level. To avoid the correspond- 
ing flavor constraints, either the masses of new particles involved must be 
quite large, or symmetries must be imposed to solve the "flavor problem" . 
For example, the squark and slepton mass matrices are flavor diagonal in 
the constrained MSSM (mSUGRA) scenario, thus suppressing such pro- 
cesses. Flavor changing neutral currents in models with universal extra 
dimensions have also been explored |lH4j . 

6 57 

The branching fraction for & ^ 57 (221 U\ I359L US] , measured at 

CLEO and BELLE, is of particular interest for supersymmetry and other 
beyond the Standard Model phenomenology. In many scenarios, the con- 
tributions to this process from new physics can add substantially to the 
Standard Model prediction. In particular, light charged Higgs bosons 
and/or charginos can be quite important for this decay 138, 401, 402]. In 
supersymmetry, the constraint is considerably stronger if < 0, but also 
relevant for n > 0, especially for large values of tan/3. 6 ^ 57 is also an 
important constraint in models of universal extra dimensions |134) . 

The branching fraction for Bs — » ^J-^ fJ'^ is quite small in the Standard 
Model (~ 3.5 x 10~^) j^^. The contribution from supersymmetry scales 
as tan^ f3, and thus becomes quite large for models with large values of 
tan p. In run I of the Tevatron, a value consistent with the Standard Model 
was found. The sensitivity of run II of the Tevatron to this quantity will 
be considerably greater. 

The anomalous magnetic moment of the muon, — 2 

In 2001, the E821 experiment at the Brookhaven National Laboratory 
reported a measurement of the muon's magnetic moment which was 2.6 
standard deviations from the Standard Model prediction JHOf • Since then, 



54 



however, an important error in the theoretical calculation was discov- 
ered which reduced the significance of this anomaly to about 1.6 stan- 
dard deviations p36, 337, 290) [113 CSl ■ With the reduction of statisti- 
cal error which has been achieved more recently, the deviation from the 
Standard Model prediction of this measurement is again about 3a using 
e+e~ data (although the significance is somewhat less using r decay data) 
[701 1174L I283L 13161 . These measurements, although somewhat difficult to 
interpret, can be used to constrain TeV-scale particle physics beyond the 
Standard Model. 

• Electroweak precision measurements 

In addition to the useful direct particle searches at LEP2, the Tevatron 
and other experiments, impressively accurate electroweak measurements 
have been made. Various limits on the scale of new physics and associated 
particle masses have been inferred from these measurements. Given these 
constraints, models with universal extra dimensions are limited to the scale 
~ 300 GeV or higher These measurements also yield particularly 

important bounds for models without a custodial SU(2) symmetry, such 
as many little Higgs models |ltj9| . 

Together, these constraints can be very powerful, often providing very tight 
bounds for specific models. For example, in Fig. 1141 we show the impact 
of collider and cosmological constraints on the constrained Minimal Super- 
symmetry Standard Model (or mSUGRA). We find that over the parameter 
space shown, constraints from LEP2 searches (Higgs, charginos and selectrons), 
along with — 2 and relic density constraints, leave only a small region near 
mi/2 = 300 - 400 GeV and mo = 80 - 150 GeV. Although the power of these 
and other constraints is quite model dependent, they are often very useful in 
supersymmetry and other classes of models. 

3.5.2 The reach of future coUider experiments 

• Future reach of the Tevatron 

The reach of the Tevatron extends to higher energies than any other accel- 
erator until the time at which the Large Hadron Collider (LHC) becomes 
operational. The range of masses which can be searched for colored par- 
ticles (squarks, gluinos and KK quarks, for example), heavy gauge bosons 
and other new physics will be increased significantly at the Tevatron lib 

[H23ini 

• The Large Hadron Collider 

The Large Hadron Collider (LHC) is expected to begin operation around 
2007 with proton-proton collisions at 14 TeV center-of-mass energy. A 
luminosity of 300 inverse femtobarns is expected to be achieved, making 
the prospects for discovering new physics at the LHC excellent. Numerous 
classes of models which provide interesting dark matter candidates will be 
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Figure 14: An example ol the impact ol collider (and cosmological) constraints on 
a model of particle dark matter I403| . The model shown is the constrained Minimal 
Supersymmetric Standard Model (cMSSM or mSUGRA) with tan P = 10, fj, > and 
Ao = 0. The almost vertical lines represent the limits on chargino (left) and Higgs 
(right) masses from LEP2. The blue dot-dash curve in the bottom left corner follows 
the 99 GeV selectron mass contour, excluded by LEP2. In the dark red region in the 
lower right, the LSP is a stau and is not, therefore, a viable dark matter candidate. 
The green region in the lower left corner is excluded by the b s'y constraint. The 
long and often narrow turquoise region provides a relic density of 0.1 < Qh^ < 0.3, 
near the observed quantity. The pink region extending over much of the lower left is 
the region within the 2a range for — 2. The two dashed curves within this region 
are the la bounds. 
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tested at this very important experiment, searching at scales of up to 
several TeV. In addition to the Higgs boson(s), the LHC will be sensitive 
to most supersymmetry scenarios, models with TeV-scale universal extra 
dimensional, little Higgs models, etc. 

For a few examples of studies which discuss the sensitivity of the LHC to 



new physics, see Refs. |SI USHl EHI E3 EEl EHl EHZl EOI • In Fig. M 



an example of such a study is shown 0^1 • It is interesting to note that in 
the region of the MSSM which is the most difficult to probe at the LHC, 
direct dark matter detection rates are very high |48) . 

• Beyond the LHC 

After the LHC, other collider experiments are likely to follow. Although 
no specific post-LHC program is certain at this time, a 500-1000 GeV 
linear collider is a possibility, perhaps followed by a Very Large Hadron 
Collider (VLHC). These or other post-LHC colliders will, of course, have 
great value to particle dark matter studies |492[ IH7| . 



4 Experiments 

4.1 Direct Detection Experiments 

Direct detection experiments appear today as one of the most promising tech- 
niques to detect particle dark matter. The idea is very simple: if the galaxy is 
filled with WIMPs, then many of them should pass through the Earth, making 
it possible to look for the interaction of such particles with matter, e.g. by 
recording the recoil energy of nuclei, as WIMPs scatter off them |2()()ll2fif)ll487| . 

The key ingredients for the calculation of the signal in direct detection ex- 
periments are the density and the velocity distribution of WIMPs in the solar 
neighborhood and the WIMP-nucleon scattering cross section. With this infor- 
mation, it is then possible to evaluate the rate of events expected in an exper- 
iment {i.e. WIMP-nucleon scattering events) per unit time, per unit detector 
material mass. 

The rate is approximately given by 




(88) 



where the index, i, runs over nuclei species present in the detector 



Detector mass 



Atomic mass of species i 



is the number of target nuclei in the detector. 



WIMP energy density 



n 



'X — 



WIMP mass 
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Figure 15: An example of the reach of the Large Hadron Collider (LHC) to new 
TeV-scale physics |49|. As a function of mo and in the mSUGRA (or constrained 
MSSM) scenario, with tan/3 = 10, Ao = Q and positive ^, the reach is shown for a 
variety of channels: zero leptons (OZ), one lepton (11), leptons with opposite charge 
[OS), leptons with the same charge (SS), three leptons (3Z), four or more leptons 
(> 41), any number of leptons plus a photon (7), at least two opposite sign leptons 
with the invariant mass within an optimized interval around the Z mass [Z l^l~) 
and an "inclusive" missing transverse energy channel. Also shown are the 2 TeV up 
squark and 2 TeV gluino mass contours. The red regions are excluded by theoretical 
constraints, while the magenta region is excluded experimentally. lOOfb"^ of integrated 
luminosity is assumed. 
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is the local WIMP density and < cr^^ > is the cross section for the scattering of 
WIMPs off nuclei of species i, averaged over the relative WIMP velocity with 
respect to the detector. For a more through discussion see, e.g. , Ref. | 319| . 

4.1.1 Scattering classifications 

The type of scattering processes considered can be classified by two impor- 
tant characteristics: elastic or inelastic scattering and spin-dependent or spin- 
independent scattering. 

• Elastic and Inelastic Scattering 

The elastic scattering of a WIMP off of a nucleus in a detector is simply 
the interaction of the WIMP with a nucleus as a whole, causing it to 
recoil, ideally often enough to measure the recoil energy spectrum in the 
target. With a Boltzman velocity distribution of WIMPs, centered at 
270 km/s, the spectrum of recoils is exponential with typical energies of 
< E >~ 50keV. Current experiments can detect recoils of considerably 
lower energy, as low as 1-10 keV. 

Inelastic scattering, on the other hand, is not observed by the recoil of 
a target nuclei. Instead, the WIMP interacts with orbital electrons in 
the target either exciting them, or ionizing the target. Alternatively, the 
WIMP could interact with the target nuclei leaving it in an excited nuclear 
state. This process leaves the signature of a recoil followed by the emission 
of a photon a nanosecond, or so, later |210| . Such signatures have to 
compete with backgrounds of natural radioactivity, however. 

• Spin-Dependent and Spin-Independent Scattering 

WIMP scattering off of nuclei is commonly discussed in the context of 
two classes of couplings. First, axial-vector (spin dependent) interactions 
result from couplings to the spin content of a nucleon. The cross sections 
for spin-dependent scattering are proportional to J(J-t-l) rather than the 
number of nucleons, so little is gained by using heavier target nuclei. For 
scalar (spin-independent) scattering, however, the cross section increases 
dramatically with the mass of the target nuclei, and typically dominates 
over spin-dependent scattering in current experiments which use heavy 
atoms as targets. 

It should be pointed out that a WIMP which is not a Majorana parti- 
cle could also scatter by vector interactions. Heavy Dirac neutrinos or 
MSSM sneutrinos are examples of particles which would scatter in this 
way. Neutralinos and Kaluza-Klein dark matter do not have such cou- 
plings, however. 

For more on scalar, axial- vector and vector WIMP-nucleon scattering, see 
Appendix [Cl 
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4.1.2 Experimental Efforts 

More than 20 direct dark matter detection experiments are either now oper- 
ating or are currently in development. In these many experiments, numerous 
techniques have been developed to measure the nuclear recoil produced by dark 
matter scattering. Some of these methods include the observation of scintillation 
(used by DAMA, ZEPLIN-I, NAIAD, LIBRA), photons (used by CREST and 
CUORICINO) and ionization (used by EDMS, GENIUS, IGEX, MAJORANA 
and DRIFT). Some experiments use multiple techniques, such as CDMS and 
Edelweiss which use both ionization and photon techniques, CRESST-II and 
ROSEBUD which use both scintillation and photon techniques and XENON, 
ZEPLIN-II, ZEPLIN-III and ZEPLIN-MAX, which use both scintillation and 
ionization techniques. 

The use of such a large array of techniques and technologies is important 
not only to accelerate the progress of the field, but also to vary the systematic 
errors from experiment to experiment, allowing for a critical assessment of a 
positive signal. 

Some experiments are also attempting to separate WIMP signatures from 
background by looking for an annual modulation in their rate. Such an effect 
would arise due to the Earth's annual motion around the Sun, resulting in a 
relative velocity relative to the galaxy's frame of reference |199| . Under this 
effect, the Earth's velocity is given by 



where tm is roughly the beginning of June and the times are in units of years. 
The result of this effect is a = 7% variation in the WIMP flux and direct 
detection rate over the course of the year. Since this variation is quite small, 
many events are needed to identify such a signature. For more on this technique 
and the status of direct detection techniques, see sectional 

4.2 Gamma-Ray Experiments 

In addition to detecting WIMPs directly, efforts are underway to attempt to 
observe the products of WIMP annihilations in the galactic halo, the center 
of the Sun or other regions. These annihilation products include neutrinos, 
positrons, anti-protons and gamma-rays. 

To observe cosmic gamma-rays directly, observations must be made from 
space. This is because in the energy range we are most interested (GeV to 
TeV), photons interact with matter via e"'"e~ pair production, which leads to 
an interaction length of approximately 38 g cm~^, which is much shorter than 
the thickness of the Earth's atmosphere (1030 g cm^^). Thus, at the energies we 
are considering, gamma-rays cannot reach ground based telescopes. Efforts have 
been developed, nevertheless, to observe gamma-rays indirectly with ground 
based experiments. In this section, we discuss the status of both ground and 
space-based gamma-ray telescopes. 




(89) 
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4.2.1 Ground-based telescopes 

When photons interact in the atmosphere, they produce an electromagnetic cas- 
cade and thus a shower of secondary particles, allowing ground-based telescopes 
to indirectly observe gamma-rays through the detection of secondary particles 
and the Cerenkov light originating from their passage through the Earth's at- 
mosphere. 

It was P. Blackett (winner of the Nobel prize in 1948) who first realized the 
possibility of detecting Cerenkov light from cosmic air showers. This realization 
was experimentally confirmed by W. Galbraith and J. Jelley (1953). Cosmic 
gamma-rays can be difficult to observe in this way, however, as most of the 
observed Cerenkov light is due to cosmic-ray induced showers with isotropic ar- 
rival directions. For detecting gamma-ray showers, an excess above the isotropic 
background of cosmic rays must be seen in the direction of a source. To accom- 
plish this, the rejection of cosmic ray showers is of crucial importance. 

To distinguish between cosmic ray and gamma-ray induced air showers, the 
observed Cerenkov light is compared with numerical simulations of atmospheric 
showers (see Fig. ll6|) . Apart from the difficulties in the treatment of interactions 
at very high energies, numerical simulations are complicated by uncertainties 
associated with the density profile of the atmosphere and the Earth's magnetic 
field. Nevertheless, reliable codes for simulating atmospheric showers exist on 
the market, for example CORSIKA and AIRES (see e.g. Knapp et al. \6'6b\ ). 




'300 -200 -100 100 200 300 -300 -200 -100 100 200 300 



Figure 16: Simulations of Cerenkov light from electromagnetic cascades initiated by a 
1 TeV photon (left) and a 1 TeV proton (right). The figures show the distribution of 
light to the ground in a 600 x 600 m^ area. Figures kindly provided by I. Perez. 

The methods of collecting Cerenkov light are quite varied, ranging from 
telescopes (and array of telescopes) to converted solar arrays. We show in table 
6 a list of existing ground based experiments for the detection of gamma-rays. 
These include imaging and non-imaging Air Cerenkov telescopes, reconverted 
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Figure 17: Sky Map for sources of very high-energy (TeV) gamma-rays, as of January 
2003 (Ref. Era]). 

solar arrays as well as experiments which detect secondary particles produced 
in showers. 

The first observation of Cercnkov light due to gamma-ray emission from 
an astrophysical source was the detection of the Crab Nebula (which today is 
regarded as the "standard candle" at these energies) with the Whipple Obser- 
vatory 10m reflector 490 . Currently, only six TeV gamma-ray sources have 
been confirmed, above 10 GeV, having been detected by multiple experiments 
at a high significance level (red symbols in Fig. I17|) . Eight sources are probable, 
i.e. detected at high significance by at least one group (blue symbols), and two 
are "possible" (light blue symbols, see Ref. |406j for more details). 

Although only a few TeV gamma-ray sources have been confirmed, many 
more could be detected in next generation experiments. Among these experi- 
ments: 

is a 17m imaging air Cerenkov telescope recently completed on 
the island of La Palma |511| . It has already started taking data. 

• "CANGAROO-III is an array of four 10 m Cerenkov telescopes being con- 
structed in Woomera, Australia j512| . It should start taking data in 2004. 

• HESS consists of four 12 m diameter Cerenkov telescopes, at a site in 
the Gamsberg area of Namibia |513| . The telescopes are operational and 
started taking data. 
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Table 6: Atmospheric Cerenkov Imaging Observatories area September 2003. 
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Figure 18: Sensitivity of present and future detectors in gamma-ray astrophysics (from 
Ref. 



IVERITAEI is an array of 7 telescopes in construction on Kitt Peak in 
Arizona, USA |514j . A preliminary version, VERITAS-4, with 4 telescopes 
should be operative around 2006. 



4.2.2 Space-based telescopes 

The first high-energy (above GeV) gamma-ray space telescopes was EGRET 
(the Energetic Gamma- Ray Experiment Telescope), onboard the Compton gamma- 
ray observatory. Launched in 1991, EGRET has observed the universe in a range 
of energies extending up to approximately 30 GeV, amassing a large catalog 
of observed gamma-ray sources, although around 60% of these sources remain 
unidentified. 

The next space-based gamma-ray observatory will be GLAST (Gamma-ray 
Large Area Space Telescope), which is scheduled for launch in 2007. As for its 
predecessor, GLAST will detect gamma-rays by recording the characteristics of 



64 



e+e" pairs produced in the interaction of the incident gamma-ray with a dense 
layer of tungsten. GLAST's effective area to gamma-rays will be a full square 
meter, considerably larger than with EGRET. GLAST will have an angular 
resolution on the order of arcminutes, compared to the degree level with EGRET 
(energy resolution varies with energy for both experiments). Unlike EGRET, 
GLAST will be sensitive to gamma-rays up to several hundred GeV in energy. 

GLAST is expected to be complementary to ground-based telescopes due to 
the lower range of energies observed, larger field of view, and higher duty cycle. 
We show in Fig^| the sensitivity of some of the present and next generation 
ground-based and space-based gamma-ray experiments. We will use such in- 
formation when discussing the prospects for indirect detection of particle dark 
matter. 



4.3 Neutrino Telescopes 

In addition to gamma-rays, neutrinos can be produced in the annihilations of 
dark matter particles. In this section, we review the status of high-energy 
neutrino telescopes, in particular, large volume Cerenkov detectors such as 
AMANDA, ANTARES and IceCube. 

Neutrinos arc considerably more difhcult to observe than gamma-rays due 
to their weak interactions with ordinary matter. Neutrinos are not easily ab- 
sorbed, however, allowing for their observation in underground, low background, 
experiments. In the GeV-TeV energy range, neutrinos are most easily observed 
by their "muon tracks" produced in charged current interactions inside of or 
nearby the detector volume. These muons travel through the detector emitting 
Cerenkov light which allows their trajectory to be reconstructed. 

For a cosmic neutrino flux, d^i,/dE^, the rate of muon tracks in a detector 
is given by 

rate= / dE, / dy A(E^)P^(E„ y; E^'") (90) 

JEj,h. Jq at I, 

where E^^'^ is the muon threshold energy of the experiment (generally between 
10 and 100 GeV), ^(-E^) is the effective area of the detector, typically in the 
range of 0.01 to 1.0 square kilometers (for further details and numerical values, 
see e.g. Ref. (S^) and Pfj_{E^, y; Ej^'^) is the probability that a neutrino of energy 
El, interacts with a nucleon producing a muon of energy i?^ = (1 — y)Ei, above 
the detector threshold energy. As one would expect, this probability depends 
on the muon range, i?(£'^, Ej^'-'), i.e. the distance travelled by muons before 
their energy drops below Ej}^''. The fmiction P^{Eu, y; Ej}^'^) is thus given by 

P^iE,., y; Efn = Na R{E^, Ef^^) '^''cciE.^y) ^g^^ 

where Na — 6.022 x 10^'^.g^^ is Avogadro's number and da^^{Eiy,y)/dy is the 
differential cross section for ncutrino-nucleon charged-current scattering. 
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The cross-section used in Eq. |^ is described in high-energy physics text- 
books, but carries uncertainties due to our hmited knowledge of parton densities. 
It can be expressed as 



with 

(-) 



= ^ {l + xys/M^)-'[{l-y)F,^'+y'xFl^'±y{l^l)xF,'^'] (93) 

where Fi = Fi{x,Q'^ xys) are the structure functions, s — 2Afjv£'(-) and 

Gp = 1.1663 X 10~^ GeV~^. For details regarding the calculation of high-energy 
neutrino-nucleon interactions including structure functions, see Ref. (246j ). 



The muon range, R{Efj^, i?™'"), appearing in Eg. 19 II follows from the energy- 
loss equation j281| 

- dE^/dX = a^iE^) + P^{E^) E^, (94) 

with X being the thickness of matter traversed by the muon, and the quanti- 
ties afj_{Ef^) and I3^{E^) are the ionization loss and the fractional energy loss 
coefficients, respectively. Integrating this result, we get the muon range 

RiE^,, Kf") = X(R7") - X{E^ = In _ (gg) 

We adopt here the following values for the coefficients: — 2.0 x 10^"^ GeV 
(cm we)^^ (cmwe = g/crp?) and = 6.0 x 10^^ (cm we)^^ |253| . 

We pass now to a brief description of existing and future neutrino telescopes, 
focusing on kilometer-scale experiments. The key idea is to detect muons, 
originating from neutrino fluxes as discussed above, building large arrays of 
photo-multipliers deep in the ice, in a lake or in the sea, to search for the 
Cerenkov light they are expected to emit as they move through these media. 

The early pioneering effort made by the DUMAND collaboration |274j was 
followed by the deployment of the Lake Baikal experiment 56 and of AMANDA 
1211 EH at the South Pole. Although these experiments have observed neutrinos 
produced in the Earth's atmosphere, they have not, thus far, identified any 
extra-terrestrial neutrinos. 

AMANDA, with approximately 50,000 square meters of effective area (at 
trigger level) and a 30 GeV muon energy threshold, has been taking data for 
several years in its current configuration. ANTARES |43II112| . currently under 
construction in the Mediterranean, will have a similar effective area and a lower 
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Figure 19: A muon neutrino event in AMANDA. Shown is the central part of the 
detector. The colorscale and symbol size correspond to hit time and amplitude |499| . 
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energy threshold of about 10 GeV. Unhke with experiments at the South Pole, 
ANTARES will be sensitive in the direction of the galactic center. 

IceCube |l8l I19| . beginning construction in 2005, and scheduled for comple- 
tion in 2010, will be the first kilometer scale high-energy neutrino telescope. 
Using technology similar to AMANDA, IceCube will be considerably more sen- 
sitive to dark matter annihilations than current experiments. Even larger, and 
perhaps lower threshold, experiments may be needed beyond IceCube to further 
search for evidence of dark matter annihilations. For a review of high-energy 
neutrino astronomy, see Refs. [2841 1353) . 

4.4 Positron and Anti-Proton Experiments 

Evidence for dark matter annihilations may also be observed in the spectra of 
cosmic positrons or anti-protons. Unlike gamma-rays and neutrinos, however, 
these charged particles do not point to their source due to the presence of galactic 
magnetic fields. Here we describe some of the experiments most important to 
these measurements. 

The HEAT (High-Energy Antimatter Telescope) experiment made its first 
balloon flight in 1994-1995, measuring the spectrum of positrons between 1 and 
30 GeV J54.. The results of this flight were very interesting, as they indicated 
an excess in the positron flux peaking at about 9 GeV and extending to higher 
energies. This excess could be a signature of dark matter annihilation in the 
local galactic halo (see section IH^jl . A second HEAT flight in 2000 confirmed 
this observation [1621 1163) . 

The BESS (Balloon borne Experiment Superconducting Solenoidal spec- 
trometer) experiment has had several successful balloon flights since 1993, pro- 
viding the most detailed measurements of the cosmic anti-proton spectrum to 
date in the range of about 200 MeV to 3 GeV gSHlEZni- Above this energy, 
up to about 40 GeV, the CAPRICE experiment provides the best anti-proton 
measurements |121| . There appears to be a mild excess in the anti-proton spec- 
trum in the hundreds of MeV range, although it is very difficult to assess this 
result with any certainty. 

In the future, the experimental sensitivity to the cosmic positron and anti- 
proton spectra is likely to improve a great deal. Perhaps as early as 2005, 
the satellite borne PAMELA experiment, will begin its mission, measuring the 
spectra of both cosmic positrons and anti-protons with considerably improved 
precision. The primary objective of PAMELA is to the measure the cosmic anti- 
proton spectrum in the range of 80 MeV to 190 GeV and the cosmic positron 
spectrum in the range of 50 MeV to 270 GeV, far beyond the energies measured 
by HEAT, BESS or CAPRICE. In Fig.EOI we show the projected sensitivity of 
PAMELA to cosmic positrons (left) and anti-protons (right). The results are 
shown assuming a contribution from annihilating neutralino dark matter. It is 
clear that PAMELA will measure these spectra to far greater precision than 
previous experiments, especially at high energies (above ~10 GeV). 

AMS (the Alpha Magnetic Spectrometer) will considerably refine the mea- 
surement of the positron spectrum in its next manifestation, called AMS-02, 
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Figure 20: The projected ability of the PAMELA experiment to measure the spectra 
of cosmic positrons (left) and anti-protons (right). A contribution from annihilating 
neutralino dark matter is included in the spectra shown. Notice, in comparison to the 
measurements made by HEAT, CAPRICE, BESS and other experiments, the dramatic 
improvement in precision. Also note the reach to higher energies made possible with 
PAMELA. From Ref. |IT7| . 

onboard the International Space Station AMS-02, with a 5000 cm^ sr aper- 
ture and a 1000 day duration, will provide exceptional precision in measuring 
the spectrum of cosmic positrons. 



4.5 Observations at Radio Wavelengths 

Observations at radio wavelengths belong to the realm of "classical" astronomy. 
Radio emission from the galactic halo, particularly from the galactic center, can 
provide a valuable probe of particle dark matter. 

Electrons and protons produced in dark matter annihilations in the Galactic 
Center region will emit synchotron radiation (at radio wavelengths) as they 
propagate through the galactic magnetic fields (see section IH^ . 

The observed Sgr A* (Galactic Center) radio emission could be explained in 
terms of synchrotron radiation emitted by shock-accelerated electrons (for more 
details see Ref. [^ffiOj and references therein). 

Rather than reviewing the subject of radio observations of the galactic halo, 
we refer to Cane (Ref. \1'67\ ) and to Brown (Ref. |131j ). which also includes an in- 
teresting discussion of the absorption of radio emission at different wavelengths. 
A complete catalog of observations of the Galactic center at all frequencies, 
and in particular at radio wavelengths, can be found in Ref. |392| . Additional 
information on specific measurements can be found in Refs. [ 1751 131L I352 | . 
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Figure 21: Current (left) and future (right) sensitivities of direct detection experi- 
ments. In the left frame, from top to bottom along the right side of the figure, the cur- 
rent limits from the CDMS, ZEPLIN-I and Edelweiss experiments are shown. The filled 
region near 30-100 GeV and 10~^^ cm'^ is the parameter space favored by the DAM A 
experiment. In the right frame, from top to bottom along the right side of the figure, 
the projected reach of the GENIUS test facility (solid), CRESST-II (dots), CDMS- 
Soudan (solid), Edelweiss-II (dashed) and ZEPLIN-MAX (dots) are shown. In each 
frame, as filled regions, the space of models predicted by supersymmetry are shown 
|258| . The narrow region along the right side of the figure represents higgsino-Iike mod- 
els, the region that reaches to the top of the figure represents mixed higgsino-gaugino 
models and the largest region represents gaugino-like models. These figures were made 
using the interface found at http://dendera.berkeley.edu/plotter/entryform.html. 



5 Direct Detection 

Many direct detection experiments have already produced quite strong limits 
on the elastic scattering cross section with protons or neutrons of potential 
dark matter candidates. Furthermore, experiments in the coming years will 
improve on current limits by several orders of magnitude making the prospects 
for discovery very great. 

Presently, the best direct detection limits come from the CDMS, Edelweiss 
and ZEPLIN-I experiments, shown in the left frame of Fig. [2] These limits are 
for spin-independent (scalar) interactions. With modern experiments, which 
use very heavy target nuclei, spin-dependent scattering experiments are not as 
sensitive to most dark matter candidates. 

Also shown in the left frame of Fig. |^ is the region in which the DAMA 
experiment claims a discovery (see e.g. Ref. for a recent review). DAMA, 
located at the INFN laboratories under the Gran Sasso mountain in Italy and 
consisting of high purity Nal crystals, has reported an annual modulation of 
their event rate consistent with the detection of a WIMP with a mass of ap- 
proximately 60 GeV and a scattering cross section on of the order of 10"^^ 
cni'^. 
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other experiments, such as EDELWEISS [HI and CDMS ED have explored 
the parameter space favored by DAMA without finding any evidence of dark 
matter. A recent model independent analysis has shown that it is difhcult to 
reconcile the DAMA result with other experiments f3W\ (see also Ref. |478| ') 
although it may still be possible to find exotic particle candidates and halo 
models which are able to accommodate and explain the data from all current 
experiments (for example, see Refs. |422l 12471 [4561 1472p . 

Theoretical and experimental results on direct detection are usually obtained 
under some simplifying assumptions on the dark matter profile. In particular, 
an isothermal profile is often assumed, with p oc (thus, with a flat rotation 
curve), a local density of po = 0.3 GeV cm~^, and a Maxwell-Boltzmann velocity 
distribution with a characteristic velocity of vq = 270 km s~^. Uncertainties in 
the density and velocity distribution of dark matter lead to the enlargement of 
the allowed region in the cross section-mass plane shown in Fig. [21 however, 
extending the mass range up to ~ 250 GeV and the cross section range down 
to (Tx-« lO-'^pb IHi [nni- If this results were confirmed, it could explain 
the discrepancy between the observational findings of different experiments. 
Unfortunately, subsequent analysises lead to different results (see in particular 
Refs. |16UI 1272] ). leaving the experimental situation unclear. 

Nevertheless, the DAMA collaboration (whose raw data are not publicly 
available) insists on the compatibility of their result with null searches of other 
experiments |88| . questioning specific experimental issues like rejection proce- 
dures and energy scale determination. 

Next generation experiments should clarify the experimental situation, thanks 
to the large improvement expected in sensitivity, around two orders of magni- 
tude in scattering cross section for EDELWEISS II and even more for ZEPLIN- 
MAX (See Fig. 1211 generated using the dark matter limit plot generator at 
^ttp://dendera. berkeley.edu/plotter/entryform. html). 

In Fig. 1221 we show the current constraints on neutralino dark matter in 
different supersymmetric scenarios. Shades paler than in the legend denote a 
value of the muon anomalous magnetic moment outside of the 2a range |283 | 
8.1 X 10-1° < (5^""y = J^'^P - Sf^ < 44.1 X 10-1° (see Ref. El for more details). 
Future generation detectors will probe a wide portion of the supersymmetric 
parameter space and will give important insights into the nature of dark matter 
particles. The B^^'' particle (KK dark matter) should have a scattering cross 
section with nucleons in the range of 10^i°-10-i^ pb, depending on its mass 
and on the mass difference with KK quark states j447| . 

6 Indirect Detection 

Indirect detection of dark matter is the technique of observing the radiation 
produced in dark matter annihilations. The flux of such radiation is propor- 
tional to the annihilation rate, which in turn depends on the square of the dark 
matter density, F^ oc p'j^M- Therefore, the "natural" places to look at, when 
searching for significant fluxes, are the regions where large dark matter densities 
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Figure 22: Current experimental sensitivity of WIMP direct detection experi- 
ments. Shown are the upper limits set by EDELWEISS yi and CDMS at Soudan 
|20|. For comparison we also show predictions for different supersymmetric scenarios 
(see text). Figure kindly provided by E. Nezri. 



accumulate. We will also refer to these regions or objects as amplifiers. 

Dense regions of the galactic halo, such as the galactic center, may be ex- 
cellent amplifiers for the purposes of detecting gamma-rays or neutrinos. Other 
astrophysical objects, such as the Sun or the Earth, could also as act as ampli- 
fiers for dark matter annihilations by capturing dark matter particles as they 
lose energy through scattering with nucleons in the interiors of these objects. 
Only neutrinos can escape these dense objects, however. Annihilation products 
which are charged move under the infiuence of magnetic fields making it impos- 
sible to consider point sources of such radiation. Despite this, observations of 
cosmic positrons and anti-protons can be valuable tools in searching for particle 
dark matter. 

In this section, we compare the predictions for gamma-ray, synchrotron, 
neutrino, positron and anti-proton fluxes from dark matter annihilation with 
current experimental data and with the expected sensitivities of future experi- 
ments. We will show that by using these techniques, it is possible to constrain 
dark matter models and, in the future, potentially detect the presence of particle 
dark matter. 
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6.1 Gamma-rays and neutrinos from the Galactic center 

One of the most interesting regions for the indirect detection of dark matter is 
the galactic center, where, according to the results of numerical simulations, the 
dark matter density profile is expected to grow as a power-law, p oc r~". The 
possible values of a, as well as alternative density profiles, have been discussed 
in chapter 121 Also recall from chapter [3 that an additional enhancement of the 
density in this region could result from the process of adiabatic accretion onto 
the supermassive black hole at the galactic center. 

Gamma-ray emission from the Galactic center has been discussed in the 
past by numerous authors (see e.g. Bouquet, Salati and Silk |128| . Stecker |461| . 
Berezinsky et al. jTS], Bergstrom, UUio and Buckley [SUI for neutralinos, Bertone, 
Servant and Sigl for B^^'> particles, Bertone and Sigl and Silk for the case 
of a density spike at the galactic center) . Here, we will review these calculations 
and arguments for evaluating the prospects for the indirect detection of dark 
matter near the Galactic center with present and next-generation experiments. 

The flux of dark matter annihilation products is proportional to the number 
of annihilations per unit time, per unit volume, (x av •n?{r) = avp'^{r)/m'^-^, 
where n(r) and p(r) are the number and the mass density of a dark matter 
particle, respectively. mDM is the dark matter particle's mass and av is its 
annihilation cross section multiplied by velocity, r is the distance from the 
galactic center. The flux is also proportional to the spectrum of secondary 
particles of species, i, per annihilation, dNi/dE. The flux observed is found 
by integrating the the density squared along the line-of-sight connecting the 
observer (the Earth) to the Galactic center . Including all factors, the observed 
flux can be written as 

Mi>,E)^av^-^^ [ dV(Ks»), (96) 

dE Airml^ Jli^c oi sight 

where the index i denotes the secondary particle observed (in this section, 
7— rays and neutrinos) and the coordinate s runs along the line of sight, in 
a direction making an angle, ip, from the direction of the galactic center. If the 
dark matter particle is not its own anti-particle (particle-antiparticle annihila- 
tion), Eq. I^^lis reduced by a factor of 2. 

In order to separate the factors which depend on the halo profile from those 
which depend only on particle physics, we introduce, following Ref. [8I| , the 
quantity J{ip): 

■^W = F^f ..^L 3 ) I dsp^ris.i;)) . (97) 



8.5 kpc V 0.3 GeV/cm'^ / Jline of sight 

We define J{Ail) as the average of J{ip) over a spherical region of solid angle, 
Afi, centered on -0 = 0. The values of J(Ari = 10~^str) are shown in the last 
column of table 7 for the respective density profiles. 

We can then express the flux from a solid angle, Afi, as 



$i(A17,£;) ~ 5.6 X 10~'^ — -[ — ][ J (Arj) Arj cm-^s~^ (98) 

dE Vpb/ V"^DM/ 
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a 


/5 


7 


R (kpc) 


7(10-3) 


Kra 


2.0 


3.0 


0.4 


10.0 


2.166 X 10^ 


NFW 


1.0 


3.0 


1.0 


20 


1.352 X 10^ 


Moore 


1.5 


3.0 


1.5 


28.0 


1.544 X 10^ 


Iso 


2.0 


2.0 





3.5 


2.868 X 10^ 



Table 7: Parameters of some widely used density profiles and corresponding value of 
J(10~'^str). For more on halo profiles, see cliapter|21 



6.1.1 Prospects for Neutralinos 

To study the detectability of gamma-ray fluxes from neutralino annihilations, 
extensive scans of the MSSM are conducted, retaining only the small minority 
of models which are consistent with accelerator and cosmological constraints. 
We show in Fig. the expected gamma-ray fluxes from the Galactic center 
for neutralino dark matter, considering an NFW halo proflle. To adapt this 
spectrum to another proflle, simply scale the flux by the value of J(Ar2 = 
IQ-"^ str) found in table 7. In Fig. [23 all continuum processes are included 
(typically dominated by annihilations to heavy quarks and gauge bosons for 
neutralino annihilation). Shown for comparison are the limit from EGRET and 
the projected reach of GLAST. For fairly heavy neutralinos, ACTs can also be 
effective. 

In addition to continuum gamma-ray emission, neutralinos can annihilate to 
mono-energetic gamma-ray lines via the processes XX ~^ 77 and XX ~^ 1^ |433) . 
Such a line, if observed, would be a clear signature for dark matter annihilation 
(a "smoking gun"). The flux of gamma-rays from such process is quite small, 
however, as no tree level feynman diagrams contribute to the process. For 
the loop-level feynman diagrams which lead to gamma-ray line emission, see 
appendix IB. 41 

The gamma-ray fluxes predicted from the Galactic center can be consid- 
erably enhanced if a density spike is considered. In Fig. [SBJ the (continuum) 
gamma-ray flux from the Galactic center is shown for four different values of 7, 
the slope of the inner halo proflle, in the presence of a density spike. For most 
of the models, a value of 7 between 0.05 and 0.1 can reproduce the scale of the 
gamma-ray flux observed by EGRET. 

6.1.2 Prospects for Kaluza-Klein dark matter 

Using the expression for the i?'"'^' annihilation cross section found in Sec. 13.3.21 
the flux of annihilation products from the Galactic center can be simplified to 

$,(A17) ~ 3.4 X 10-^^^ ( J {An) An cm-2 . (99) 
dE \mB(i) J 
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Figure 23: The flux of gamma-rays above 1 GeV per square meter per year from the 
Galactic center from annihilations of neutralino dark matter. A NFW halo profile 
has been used. For each point, the thermal relic density is below the maximum value 
allowed by WMAP. The solid and dashed lines are the limit from the EGRET exper- 
iment and predicted sensitivity for GLAST, respectively lljU2l . The various shadings 
refer to different scenarios of supersymmetry breaking. For more information, see 
Ref. EIO]. 
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Figure 24: As in the previous figure, but from tlie annifiilation of neutralino dark mat- 
ter to a 77 line. The dashed hue represents the predicted sensitivity for GLAST. The 
flux for neutralino annihilation to 7Z is similar. For more information, see Ref. )31U| . 
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Figure 25: Predicted gamma-ray fluxes for a large set of supersymmetric models and 
halo proflles: 7=0.05(triangles), 7=0.12(diamonds), 7=0.2 (dots), 7=1.0 (squares). 
The flux observed by EGRET |392| is shown as grey boxes. Also shown are the 
projected sensitivities of for GLAST (1 month observation time) and MAGIC (50 
hours). 

In Fig. 1261 the predicted 7— ray flux from KK dark matter annihilations in the 
Galactic center is shown. Results for LKP masses of 0.4, 0.6, 0.8, and 1 TeV are 
shown. A halo profile with J (lO~^) = 500 has been used, although the effect 
of this choice is easily scaled with the values shown in table 7. 

Unlike in the case of supersymmetry, with Kaluza-Klein dark matter, there 
are few free parameters in calculating the gamma-ray spectrum from the Galac- 
tic center (m^ci) and J{AVl)). We can, therefore, easily place limits on the 
halo profile as a function of the LKP mass. We show in Fig.|23the constraints 
on these parameters based on the expected sensitivity of GLAST, MAGIC and 
HESS. For example, with an NFW profile, LKP masses below about 600 GeV 
will be excluded if MAGIC does not observe any radiation from the galactic 
center. 

Neutrino telescopes will also be capable of searching for signals of dark mat- 
ter annihilation in the Galactic center (see Sec. I4.3|l . although these prospects 
are considerably poorer. In Fig. 1281 we plot the integral fiux of muon neutrinos 
above 50 GeV (solid line) as a function of the B^^^ mass. This result is obtained 
by adding the neutrino fiuxcs from three different channels: 
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Figure 26: Expected 7— ray fluxes for (top to bottom) rrigd) = 0.4, 0.6, 0.8, and 1 
TeV and J(l0"^) = 500. For comparison shown are typical 7— ray fluxes predicted 
for neutralinos of mass ~ 200 GeV, (shadowed region) as weU as EGRET |377| data 
and expected sensitivities of the future GLAST|435j, MAGICEig and HESSEHSl 
experiments. 
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Figure 27: Value of J = J(10 ^) required to produce 7 fluxes observable by the 

future GLAST, MAGIC and HESS experiments, as a function of the S'^' mass. For 
comparison we show the values of J for some profiles discussed in the text. 
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Figure 28: Integral flux of 7— rays (sloped dotted line) and muon neutrinos (solid) 
above 50 GeV, for J (lO~^) — 500 . The dashed line shows the contribution of direct 
B'^-* annihilation into neutrinos. Horizontal lines are sensitivities of present and future 
experiments for 7— rays (dotted horizontal lines) and neutrinos (upper solid line). 

• Neutrinos produced directly in B^^^ annihilations (dashed line), their spec- 
trum being a line at energy = mg(i) . 

• Secondary neutrinos from decay of charged pions. This spectrum can 
be evaluated using the expressions for the charged pion decay found in 
Ref. [HSl- 

• Secondary neutrinos from "prompt" scmi-leptonic decay of heavy quarks 
(solid line). This spectrum is given, for example, in Ref. |818| . 

We show in the same figure an estimate of the sensitivity of the neutrino 
telescope ANTARES (upper solid line). To estimate this sensitivity, we first 
evaluated the rate of muons in ANTARES from the direction of the galactic 
center, which depends (see Eq. I90() on specific experimental quantities, such 
as the detector effective area and the threshold energy for the detection of 
muons. The rate is higher for more energetic neutrinos, being proportional to the 
muon range and to the neutrino-nucleon cross section, which are both increasing 
functions of energy. See Sec. 14.31 for a discussion of neutrino telescopes. 

Of these neutrino producing channels, ANTARES is most sensitive to those 
neutrinos produced directly in B^^'^B^^'^ — > w annihilations. Although the 
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branching ratio for this channel is an order of magnitude smaller than that 
into quark pairs, the neutrinos produced are emitted at the highest available 
energy, E^, = m^d), increasing their probability of being detected. Neutralinos, 
which do not annihilate directly to neutrinos, are more difficult to observe with 
neutrino telescopes. 

In Fig. 1281 the integral flux of gamma-rays is shown for comparison, along 
with the projected sensitivity of future experiments GLAST, MAGIC and HESS. 

6.1.3 The gamma-ray source at the Galactic center 

The EGRET experiment (see Sec. 14. 2. 3|) has reported an excess of gamma-rays in 
the region of the galactic center, in an error circle of 0.2 degree radius including 
the position l=0(deg) and b=0(deg) |377| . The name for this source, in the 
language of high-energy astrophysicists, is 3EG J1746-2851. The radiation is 
well above the gamma-ray emission which would be expected from interactions 
of primary cosmic rays with the interstellar medium (see, e.g. , Refs. |464lll4j] l. 

It is intriguing to imagine that such excess emission could be the product of 
dark matter annihilations near the Galactic center . However, it should be noted 
that some difficulties exist, related to this interpretation. In fact, as shown in 
Refs. 3,02 , ^03 , the EGRET source is not exactly coincident with the galactic 
center. This makes the interpretation of the EGRET signal as dark matter 
annihilation in a density spike (see Sec l2.5.l'|l problematic. 

Furthermore there is some evidence, though weak, that the source could 
be variable. Such a result could rule out completely the interpretation of the 
excess emission as due to annihilation radiation from the galactic center. The 
variability of 3EG J1746-2851 has been recently discussed in Ref. .399,. 

We briefly note here that multiple Atmospheric Cerenkov Telescopes (ACTs) 
have recently reported an excess of gamma-rays from the Galactic center region. 
The VERITAS collaboration, using the Whipple telescope, have reported a flux 
of 1.6±0.5±0.3x 10-**m-2s-i above 2.8 TeV [MSI- The CANGAROO collab- 
oration has reported a flux of approximately 2 x 10~^ m~^ from this region 
in the range of 250 GeV to 1 TeV ,471,. We eagerly await the resuls of HESS, 
which should be the most sensitive to the Galactic center region. 

It is certainly too soon to determine whether the fluxes observed by these 
experiments are the product of dark matter annihilations or are the result of 
another process, most likely astrophysical |307| . Improvements in the mea- 
surement of the gamma-ray spectrum, and improved angular resolution will be 
needed to resolve this issue. 

6.1.4 Upper Limit for the Neutrino Flux from the GC 

Despite the large uncertainties associated with the distribution of dark matter 
in the innermost regions of our galaxy, it is possible to set an upper limit on 
the neutrino flux by requiring that the associated gamma-ray emission does not 
exceed the flux observed by EGRET (see previous section) 
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The maximum neutrino flux is obtained by normalizing the associated flux 
of gamma-rays to the EGRET data. This corresponds to fixing, for each model, 
the product, av N^, with = Ri- Here Ri is the branching ratio of all 

the channels, i, contributing Ni gamma-rays above a given threshold energy. 

Having fixed the particle physics contents of our dark matter candidate, the 
ratio between the number of photons and the number of neutrinos emitted per 
annihilation is known. 

The rescaled flux of muons, 0"°"°(> Eth), will thus be given by 



^NFW(>j^^^) REGRET (jg^ 



C™(> Et.) - ' ^NFw; ^ ^ ' ' . (100) 



where the label NFW indicates that NEW profiles have been used to compute 
profile-independent flux ratios and E^ is the energy at which we decide to nor- 
malize the flux to the gamma-ray data (in our case, E^, = 2 GeV). 

The results are shown in Fig. 123 where shades paler than in the legend 
denote a low value for the muon's magnetic moment (see Ref. 9A for more 
details). The neutrino induced muon flux normalized to the EGRET data rep- 
resents an upper limit, as the observed gamma-ray emission certainly could be 
due to processes other than dark matter annihilation. The comparison with 
the sensitivity of ANTARES shows that only the highest mass neutralinos can 
possibly be detected with neutrino from the galactic center. In this case, conser- 
vatively assuming that the gamma-ray emission observed by EGRET is entirely 
due to neutralino annihilation, the upper limit on the neutrino flux is barely 
above the minimum signal observable by ANTARES in 3 years. 

If neutrinos are nevertheless observed above the given fiuxes, then their 
interpretation as due to neutralino annihilation is problematic and would ac- 
tually require either the adoption of other dark matter candidates annihilating 
dominantly into neutrino pairs or a different explanation, e.g. in terms of as- 
trophysical sources. 



6.2 Synchrotron Radiation from the Galactic Center 

Another interesting means of indirect detection of dark matter is observing the 
synchrotron radiation originating from the propagation of secondary e^'s in the 
galactic magnetic fields. We will focus on what happens at the center of our 
galaxy, where most of the annihilation signal comes from. 

The magnetic field around the Galactic center is thought to be at equiparti- 
tion, i.e. there is equipartition of magnetic, gravitational and kinetic energy of 
the plasma surrounding the central supermassive black hole (see Sec. 12.5.1)1 . 

It is easy to derive the strength of the magnetic field under a few simplifying 
assumptions. Let us consider the existence of a galactic wind of particles with 
velocity, Vg^ . These particles will be captured by the gravitational potential well 
of the black hole at the center of the galaxy within the accretion radius, = 
2GM/Vg^, where M is the mass of the central object. Under the assumption of 
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Figure 29: Neutrino-induced muon flux from tlie Galactic center normalized to 
EGRET. Models are sorted by relic density (left) and leading annihilation channel 
(right). Shades paler than in the legend denote a low 5'^'^'^ value. 



purely radial infall, the radial dependence of the particle velocity is simply 



rs 



1/2 



where rs = GM/c^ is the Schwarzschild radius of the black hole. For a steady 
flow, the number of plasma particles through a sphere of radius, r, around the 
black hole is 

N= — = -ATrr'^n{r)v(r), (102) 
at 

where n(r) is the plasma number density. Solving for n(r) we find 
The accretion rate can be parametrised as follows 



TV /rc\3/2 

"W = = (f) ■ m) 



Nmp w lO^^g s-^WM^, (104) 

where nip is the proton mass, Mq is the mass of the central black hole in units of 
10^ Mq, and W ~ 1 parameterizes the uncertainties of the physical parameters 
around the black hole, namely the velocity and mass-loss rate of the circum- 
nuclear wind (see Ref. |379j for more details). 

We now note that the infalling plasma is expected to be highly ionized, and 
its energy density will reach equipartition with kinetic and gravitational ener- 
gies. After simple calculations, we find that the magnetic field under equiparti- 
tion evaluates to [3791 145U | 

/ \ -5/4 

B = 324^G f j . (105) 
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Far from the galactic center, the equipartition does not apply anymore and 
we assume a flat profile for the magnetic field. More specifically, the profile we 
adopt is 

-5/4 



B{r) — max 



r 



32VG ( — ) , e^iG 



(106) 



which means that the magnetic field is assumed to be in equipartition with the 
plasma out to a galactocentric distance Tc = 0.23 pc and equal to a typical value 
observed throughout the galaxy at larger distances. 

If the actual value of the magnetic field away from the central region was 
smaller than the value we considered, this would imply a shift of the radio 
spectrum toward lower energies and thus, in the range of frequencies we are 
interested in, a higher flux for a given frequency. This would also translate into 
stronger constraints for the mass and annihilation cross section. Nevertheless, 
we prefer to be conservative and consider a quite high value of B. Note that 
magnetic fields stronger than equipartition values are physically unlikely. 

The mechanism of synchrotron emission is well known. We will now esti- 
mate the synchrotron luminosity produced by the propagation of secondary 's 
originating from dark matter annihilation in the galactic magnetic field 

We recall that the critical synchrotron frequency, Vc{E), i.e. the frequency 
around which the synchrotron emission of an electron of energy, E, in a magnetic 
field of strength B, peaks, can be expressed as 

M^, = JLi£(^)\ (107) 

47r irieC \ irieC / 

where me is the electron mass. Inverting this relation, we determine the energy 
of the electrons which give the maximum contribution at that frequency, 

To compute the synchrotron luminosity we also need to know their energy 
distribution, which in our case can be expressed as (see, e.g. , Ref. j262p 

dn ru> E) 

dE " P{E) 
where T is the annihilation rate, 

/ pL4vrr2 dr. (110) 



™2 / t'sp^ 
'"dm "'0 



The function fe{r) is given by 



rsp 



= 2 7 2 ^ (111) 

Jo P;p47rr^ dr 
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and 

-<-^^ 

is the total synchrotron power spectrum. Note that the general expression for 
/e(r) would have to take into account spatial redistribution by diffusion (see, 
e.g. , Ref. |262j ). but this is typically negligible ^T\ . 

The quantity Ye{> E) is the number of e^'s with energy above E produced 
per annihilation, which depends on the annihilation modes, and can be equiv- 
alently expressed as a function of the particle mass, ttidm, and the frequency 
considered, v. Actually, Eq. H1U8|) shows that for the frequencies we are inter- 
ested in, Em{v) « todm, and thus the energy dependence of Ye{> E) can be 
neglected. We estimate Ye{> E) by the number of charged particles produced 
in quark fragmentation (see below for further details). 

For each electron the total power radiated in the frequency interval between 
V and u -\- dv IS given by 

P{v, E) ^ ^B{r)^ r K,/,{y)dy = ^B{r)F f-J— ) , (113) 

where Kn{y) are the modified Bessel functions of order n (for definitions see 
e.g. gni) and 



Integrating this formula over the dark matter distribution, we obtain the total 
synchrotron luminosity 

/•oo i-muM J 

= / dr Anr^ / dE ^P{i^,E), (115) 
which by substitution becomes 

L,^^r /""" ,E^i^F (^) . (116) 

meC^ Jo Jm, P{E) \Vc{E)J 

It is possible to simplify this formula by introducing the following approxi- 
mation (see Rybicki & Lightman j434| '): 

5{u/uc{E)~Q.2Q) . (117) 



ME) 

The evaluation of the integral then gives 

.3^5 \ 1/2 



L.(^) ^ l-l (^^) ^r.(mDM,-) M" m , (118) 
47r 8 \0.297r e / "^dm 

where 

/•oo 

/(V-)-/ ds p'{T{s,^))B-''^r{s,rP)) (119) 
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and s is the coordinate running along the line-of-sight. 

For frequencies around 400 MHz, used below, and for the lowest value of the 
magnetic field, we find that Em (400 MHz) ^ 2 GeV. In reality, for dark matter 
profiles with central cusps, e.g. the NFW, Kravtsov, and Moore profiles, most 
of the annihilation signal comes from the inner region of the galaxy, where the 
magnetic field is probably higher. For v — 400 MHz and r < rc, 

which at the inner edge of the profile, corresponding to the Schwarzschild ra- 
dius of the supermassive black hole at the galactic center, Rs = 1.3 x 10^^ 
pc, takes the value of i?,„(400MHz) = 2.2 x 10~^ GeV. We thus always have 
Ejn (400 MHz) <C TTiDM , which means that most of the secondary electrons are 
produced above this energy and contribute to the radio flux. 

For a particle of mass todm, the average electron multiplicity per annihi- 
lation, Ye(TODM), is evaluated by adding the contribution of each annihilation 
channel with cross section (cru)i, producing Y^(mj)M) electrons: 

avYe{m-DM) = y^^{av)iY^{m-DM) , (121) 

i 

where av is the total annihilation cross section. 

The main channels contributing to this flux are direct production of leptons 
and annihilation into quarks. The calculations are easily performed for Kaluza- 
Klein dark matter. For direct production of leptons, 

r/* (M) = y/"* (M) ~ 2 (122) 

in the relevant range of masses. In the quark channel, to count the number of 
electrons, Y^^{M), we integrate the fragmentation functions for e^'s from tt^^s. 
This results in 

avYe{l TeV) ~ 6 x 10"^ TeV"^ (123) 

and 

ye(lTeV) ~ 4.5, (124) 

for mDM = 1 TeV. The electron multiplicity in the hadronic channel alone would 
be much larger, roughly 20. 

The case of neutralinos is much more complicated, as the dominating anni- 
hilation modes can vary from model to model (a discussion of branching ratios 
in the framework of the mSUGRA models can be found in Ref. ^2). Such 
calculations must, therefore, be conducted on a model- by- model basis. 

One more step is necessary to calculate the observed radiation. We must 
multiply the synchrotron luminosity, L,^, with the synchrotron self-absorption 
coefficient, which we calculate next. 

Synchrotron emission is accompanied by absorption, in which a photon 
loses its energy due to the interaction with a charge in a magnetic field. The 
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synchrotron self-absorption coefficient is by definition (see Rybicki and Light- 
man gni) 

A^^— (I~e-^(''))27r6d6, (125) 



a 



V Jo 



where t(&) is the optical depth as a function of the cylindrical coordinate b, 

1*00 

T{b)=a, fe{b,z)dz, (126) 
Jdib) 



and the coefficient, a^, is given by 



e^ri?(r) r j^2± 



E^P{E) 



F\ — ] dE. (127) 



The final luminosity is obtained by multiplying Eq. I|115() with A^, given by 
Eq. (|125|l . It is evident that in the limit of small optical depth, the coefficient 
Au — > 1, as can be seen by expanding the exponential. 

The lower limit of integration of Eq. (|126|l is 

d{b) = for b^ + > (4i?s)^ (128) 



d{b) = ^(4i?s)2 - 52 elsewhere. 
Using the approximation introduced in Eq. H117|l . we find 



rr 



(129) 



which can in turn be used to evaluate r(&) in Eq. H126|) and A^y in Eq. (|125|l . 

If a density spike exists at the galactic center, the self- absorption coefficient 
cannot be neglected and can lead to a significant reduction of the observed syn- 
chrotron flux by up to several orders of magnitude. The results of synchrotron 
emission in presence of a spike has been discussed in Ref. 'PIJ. If a spike exists 
at the galactic center, and if neutralinos are the dark matter particle, only small 
values of 7 are compatible with radio observations. Kaluza-Klein dark matter 
has annihilation cross section typically larger than neutralinos. Thus Kaluza- 
Klein dark matter is very problematic in a scenario with a density spike. 

Note that Eqs. I125II127I are valid, strictly speaking, only for position in- 
dependent quantities. A rigorous treatment of synchrotron emission and self- 
absorption would require the solution of the radiative transport equation. Re- 
cently, Aloisio et al. 2004 Ref. |23 have derived the equilibrium distribution of 
electrons and positrons from neutralino annihilation at the Galactic center, and 
the resulting radiation considering adiabatic compression in the accretion flow, 
inverse Compton scattering of synchrotron photons (synchrotron self-Compton 
scattering), and synchrotron self-absorption. Such a detailed analysis allows a 
more precise estimate of the radio emission and confirms that neutralino anni- 
hilation in a NEW profile with a spike would exceed the observed radio emission 
from the Galactic center. 
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If there is no spike at the galactic center, the optical depth is negligible and 
the self- absorption coefficient is of the order of unity. In fact, using Eg. 11291 the 
optical depth in Ea. 11261 can be expressed as 



<7V Ye{M) 1 



dsp^{s), (130) 



where ~ 8kpc is the distance of the Sun from the galactic center. Using 
m-DM = ITeV, (TV — 1.6 x 10~'*TeV~^ (for the cross section for annihilation 
into right-handed up quarks) and a NFW halo profile, we find 

T ^ 1.78 X 10-^ ( J^ ) \ (131) 



lOOMHz 

We can thus neglect self-absorption unless the frequency considered is very 
small. The absorption on relativistic electrons from other sources is also negli- 
gible. Using 

n{E) ^ 10"^ GeV^^cm^^s^^sr"^ (132) 

for the locally observed differential electron flux (see Ref. (387j ) in the relevant 
energy range given by Eq. (|120|l . one obtains an absorption coefficient per length 

o.S6xlO-'«pc-'(|.)(g^)". (133) 

Even if the relativistic electron flux due to non-acceleration processes close to the 
Galactic center is orders of magnitude larger, this effect would still be negligible. 
However, for frequencies below a few MHz, free-free absorption is important (see 
e.g. Ref. [nZI). 

To compare with observations, we integrate over the relevant solid angle. 
The comparison between predicted and observed fluxes can constrain the cross 
sections and masses of annihilating dark matter particles for a given halo proflle. 
In particular, this method be used to provide a lower bound on the mass of a 
Kaluza-Klein dark matter particle. In Fig. 1301 we show predicted and observed 
fluxes for Kaluza-Klein particles, for a NFW profile, as a function of the particle 
mass. Three cases are shown (see Ref. |02| for more details). For each case the 
predicted and observed fluxes are plotted, the latter being represented by a 
horizontal line. The three cases are represented by solid, dashed and dotted 
lines. Case 1 is the most constraining, implying a lower bound on the mass of 
about 0.3 TeV (assuming an NFW halo profile). 

The fluxes predicted at high latitudes can also be compared with observa- 
tions (see Ref. |137| '). The strongest constraints result from the lowest frequen- 
cies at which free-free and synchrotron self- absorption are not yet important, 
i.e. ~ 10 MHz. Here, the observed background emission between 0° and 90° 
from the galactic anti-center is ~ 6 x 10^ Jy. Comparing with the predicted 
emission results in the limit 

^ in-22 f rnBM \'^ ye(lTeV) 3 _i 

CTW ^ 10 -— -cm s . 134 

V TeV / Yeimuu) 



88 



10^ - 



102- 




Figure 30: Predicted (curves) and observed (horizontal lines) radio flux from annihi- 
lating Kaluza-Klein dark matter from regions close to the galactic center. An NFW 
profile has been assumed. Three cases for the radio frequency and observed region are 
shown. 

While this is considerably weaker than the constraints above, it is largely inde- 
pendent of the unknown halo profile near the galactic center. 

The argument can also be turned the other way round and interpreted as a 
measure of the galactic magnetic field. One can in fact decide to select the values 
of 7 reproducing the normalization of the observed gamma-ray emission and look 
for the values of the magnetic field that reproduce the correct normalization of 
the observed radio emission. Due to the uncertainties in the particle physics 
models, the precision obtained is unfortunately poor (see Ref. I92j ) . 

A comparison between the prospects for indirect detection at different wave- 
lengths shows that constraints from synchrotron emission are more stringent 
than those obtained from gamma-ray and neutrino fiuxes (from current ex- 
periments), although they are less robust, being derived under the additional 
assumptions made about the strength of the magnetic field at equipartition. In 
the future, however, the strongest constraints will likely result from gamma-ray 
experiments. 

6.3 Annihilation Radiation from External or Dwarf Galax- 
ies 

It also might be possible to observe annihilation radiation from galaxies outside 
of the Milky Way. In this case, although such galaxies are far more distant 
than the galactic center, the observed emitting region is much larger. Dwarf 
galaxies within the Milky Way may also be observable regions of dark matter 
annihilation. 

Baltz et al. studied the expected flux of gamma-rays from M87 and sev- 
eral local dwarf spheroidal galaxies. The predicted fluxes are calculated using 
the same formulae for indirect detection of secondary particles described above. 
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Using a profile with a central core for the sources, they conclude that predicted 
fluxes are below the sensitivities of next-generation experiments, unless the an- 
nihilation signal is boosted by a significant amount of dense clumps. A similar 
analysis was carried on for the prospect of observing M31 with CELESTE (see 
Falvard et al. [223 )• 

Tasitsiomi, Gaskins and Olinto |466) focused on gamma-ray and synchrotron 
emission from the Large Magellanic Cloud (LMC). Fitting the LMC rotation 
curve with different profiles, they determined that although present data do 
not constrain SUSY parameters, future experiments like GLAST (gamma-rays), 
and LOFAR (low frequencies emission) , could probe a significant portion of the 
SUSY parameter space. Similar, though more optimistic, conclusions have been 
obtained by Fieri and Branchini |418| . 

Finally, particularly interesting are the prospects for dark matter observa- 
tions in the Draco and Sagittarius dwarf galaxies |224l I47It] . 

6.4 High-Energy Neutrinos from the Sun or Earth 

In addition to gamma-rays, neutrinos can be produced in dark matter annihi- 
lations. Unlike gamma-rays, however, neutrinos can escape from dense media 
in which such annihilations may take place. For example, WIMPs which are 
captured in deep gravitational wells such as the Sun or Earth, can annihilate at 
great rates. Although gamma-rays cannot escape these objects, neutrinos often 
can, providing an interesting signature to search for with high-energy neutrino 
telescopes 0551 1^ O 15^ 1^ 1^ 

6.4.1 Capture and annihilation in the Sun 

In order to provide an observable flux of neutrinos, dark matter particles must 
be gathered in high concentrations. In the following calculation, we will focus 
on WIMP capture in the Sun, as these prospects are more promising than for 
capture in the Earth. 

The rate at which WIMPs are captured in the Sun depends on the nature 
of the interaction the WIMP undergoes with nucleons in the Sun. For spin- 
dependent interactions, the capture rate is given by (2691 1319] 

~ 3 3^xin20s-i f ^ /^ 270km/s y / an.SD \ /lOOGeVy 

U-3GeV/cm3j 1, ) [lO-e^h)[ muu ) ' 

(135) 

where piocai is the local dark matter density, cth.sd is the spin-dependent, 
WIMP-on-proton (hydrogen) elastic scattering cross section, wiocai is the local 
rms velocity of halo dark matter particles and muM is the dark matter particle's 
mass. 
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The analogous formula for the capture rate from spin-independent (scalar) 
scattering is EHI 



a . 1.24 X 10'".-' I C^)' C-^)' X 

^ 0.3GeV/cmV V "local / V "^dm / 

2.6gH,si + 0.175 Q-Hc, SI 

10-6 pb 



(136) 



Here, ctr, si is the spin- independent, WIMP-on-proton elastic scattering cross 
section and (Tho.si is the spin- independent, WIMP-on-helium, elastic scattering 
cross section. Typically, ctho.si — 16.0 (Th.si- The factors of 2.6 and 0.175 include 
information on the solar abundances of elements, dynamical factors and form 
factor suppression. 

Although these two rates appear to be comparable in magnitude, the spin- 
dependent and spin-independent cross sections can differ radically. For exam- 
ple, for Kaluza-Klein dark matter, the spin-dependent cross section is typically 
three to four orders of magnitude larger than the spin-independent cross section 
| 446[ 1150) and solar accretion by spin-dependent scattering dominates. Spin- 
dependent capture also dominates for most neutralino models. On the other 
hand, for scalar dark matter candidates (such as sneutrinos, or candidates from 
theory space little Higgs models), the spin-independent cross section can domi- 
nate. 

If the capture rates and annihilation cross sections are sufficiently high, equi- 
librium may be reached between these processes. For N WIMPs in the Sun, the 
rate of change of this quantity is given by 

iV = C®-A®iV2, (137) 

where C® is the capture rate and A® is the annihilation cross section times the 
relative WIMP velocity per volume. C® was given in Eq. 11351 while A® is 

A® = M, (138) 

where Vcs is the effective volume of the core of the Sun determined roughly 
by matching the core temperature with the gravitational potential energy of a 
single WIMP at the core radius. This was found in Refs. .280^ .270) to be 



Fe« = 5.7xlO-cm3(l^V'\ (139) 



The present WIMP annihilation rate is given by 

r = ^A^N^ = ^ tanh^ (VCQA® t©) , (140) 

where ~ 4.5 billion years is the age of the solar system. The annihilation rate 
is maximized when it reaches equilibrium with the capture rate. This occurs 

when 

VC®A0iQ > 1 . (141) 
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For many of the particle physics models which are most often considered (most 
supersymmetry or Kaluza-Klein models, for example), the WIMP capture and 
annihilation rates reach, or nearly reach, equilibrium in the Sun. This is of- 
ten not the case for the Earth. This is true for two reasons. First, the Earth 
is less massive than the Sun and, therefore, provides fewer targets for WIMP 
scattering and a less deep gravitational well for capture. Second, the Earth 
accretes WIMPs only by scalar (spin-independent) interactions. For these rea- 
sons, it is unlikely that the Earth will provide any observable neutrino signals 
from WIMP annihilations in any planned experiments (for a recent analysis of 
WIMP capture in the Earth, see Ref. ^U^). 

The flux of neutrinos produced in WIMP annihilations is highly model de- 
pendent as the annihilation fractions to various products can vary a great deal 
from model to model. We will attempt to be as general in our discussion as 
possible while still considering some speciflc cases as well. 

In supersymmetry, there are no tree level diagrams for direct neutralino 
annihilation to neutrinos. Many indirect channels exist, however. These in- 
clude neutrinos from heavy quarks, gauge bosons, tau leptons and Higgs bosons. 
These processes result in a broad spectrum of neutrinos, but with typical en- 
ergies of 1/2 to 1/3 of the neutralino mass. For experimental (muon) energy 
thresholds of 10-100 GeV, lighter WIMPs can be very difficult or impossible to 
detect for this reason. 

For neutralinos lighter than the mass (80.4 GeV), annihilation to bb 
typically dominates, with a small admixture of t~^t~ as well. In these cases, 
neutrinos with less than about 30 GeV energy are produced and detection is 
difficult. For heavier neutralinos, annihilation into gauge bosons, top quarks 
and Higgs bosons are important in addition to bb and r+r" . In particular, 
gauge bosons can undergo two body decays {Z — > or l^v) producing 

neutrinos with an energy of about half of the WIMP mass. Neutralinos with a 
substantial higgsino component often annihilate mostly into gauge bosons. 

For Kaluza-Klein dark matter, the picture is somewhat different. Kaluza- 
Klein dark matter particles annihilate directly to a pair of neutrinos about 3-4% 
of the time |446lll50ll305| . Although this fraction is small, the neutrinos are of 
higher energy and are, therefore, easier to detect. The more frequent annihila- 
tion channels for Kaluza-Klein dark matter are charged leptons (60-70%) and 
up-type quarks (20-30%). Of these, the r+r~ mode contributes the most to 
the neutrino flux. Unlike in supersymmetry, a large fraction of lightest Kaluza- 
Klein particles annihilate into long lived particles, such as up quarks, electrons 
and muons, which lose their energy in the Sun long before decaying. Bottom 
and charm quarks lose some energy before decaying, but not as dramatically. 

Neutrinos which are produced lose energy as they travel through the Sun 
| 204[ I31ijl I168| . The probability of a neutrino escaping the sun without inter- 
acting is given by 

P = e-^'^/^" (142) 

where Ek is ~ 130 GeV for ly^, ~ 160 GeV for i/^, ~ 200 GeV for and 
~ 230 GeV for i7t- Thus we see that neutrinos above a couple hundred GeV are 
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especially depleted. For a useful parameterization of solar effects, see Ref. |2(J4) . 
Note that neutrino oscillations can also play an important role in calculating 
the flux of muon neutrinos in a detector | lti8 |. 

6.4.2 Detection of high-energy neutrinos from the Sun 

Several experiments are potentially able to detect the flux of high energy neutri- 
nos from dark matter annihilations in the solar core. The AMANDA experiment 
is currently the largest operating neutrino telescope. The AMANDA B-10 array, 
due to its "soda can" geometry, was not very sensitive in the direction of the 
Sun (the horizon), although the current version of the experiment, AMANDA- 
II, does not have this problem and can place limits on dark matter annihilations 
from the center of the Sun and Earth. ANTARES, with a lower energy threshold 
(10 GeV) and IceCube, with a much greater effective area, will each function as 
effective dark matter experiments (see section [4.31 for a description of neutrino 
telescopes). 

The background for this class of experiments consists of atmospheric neutri- 
nos |243j and neutrinos generated in cosmic ray interactions in the Sun's corona 
|82l . In the direction of the Sun (up to the angular resolution of a neutrino 
telescope), tens of events above 100 GeV and on the order of 1 event per year 
above 1 TeV, per square kilometer are expected from the atmospheric neutrino 
flux. The rate of events from neutrinos generated by cosmic ray interactions in 
the Sun's corona is predicted to be less than a few events per year per square 
kilometer above 100 GeV. 

The sensitivity of a square kilometer neutrino detector with a moderate muon 
energy threshold (50 GeV) to supersymmetric dark matter is shown in Fig. 1311 
From this figure, it is clear that high-energy neutrinos will be an observable 
signature in only a small fraction of possible supersymmetry models, although 
such experiments are still certainly an important probe. 

For Kaluza-Klein dark matter, the prospects for detection via high-energy 
neutrinos are substantially better. This is largely due to the dominating annihi- 
lation modes. The spectrum of muons in a detector due to LKP annihilations in 
the Sun is shown in Fig.|221for various annihilation channels and for two choices 
of LKP mass. Unlike in the case of supersymmetry, annihilation to neutrinos 
and taus dominates the neutrino spectrum. In supersymmetry, b quarks and 
gauge bosons dominate, producing fewer observable neutrinos. 

In Fig. 1331 the event rates from Kaluza-Klein dark matter annihilation in the 
Sun are shown for a square kilometer detector with a threshold of 50 GeV. Each 
of the three lines correspond to variations in the Kaluza-Klein spectrum. For 
the spectrum predicted in Ref. |149j , a kilometer scale neutrino telescope could 
be sensitive to a LKP with mass up to about 800 GeV. The relic density of the 
LKP varies from low to high values from left to right in the graph. The range 
of masses of the LKP that gives the appropriate relic density was estimated 
from Refs. [1501 1446j and shown in the figure by the solid sections of the lines. 
Combining the expected size of the one-loop radiative corrections with a relic 
density appropriate for dark matter, we see that IceCube should see between a 
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Figure 31: The number of events from neutralino aimihilation in the Sun per year 
in a neutrino telescope with an effective area equal to one square kilometer and a 50 
GeV muon threshold [201 • The Ughtly shaded region represents the general Minimal 
Supersymmetric Standard Model (MSSM), the darker region corresponds to mSUGRA 
models, a subset of the MSSM. For each point shown, the relic density is below the 
maximum value allowed by the WMAP data (Q.-^h^ < 0.129). The sensitivity projected 
for IceCube is shown as a dashed line |205l . 
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Figure 32: The spectrum of muons at the Earth generated in charged-current inter- 
actions of muon neutrinos generated in the annihilation of 600 GeV (left side) and 
1000 GeV (right side) Kaluza-Klein dark matter particles in the Sun |305| . The elastic 
scattering cross section used for capture in the Sun was fixed at 10~^ pb for both 
graphs. The rates are proportional to that cross section. 



few events and tens of events per year. 

For detectors with smaller effective areas one simply has to scale the curves 
down by a factor A/{1 km^) to obtain the event rate. In particular, for the 
first generation neutrino telescopes including AMANDA and ANTARES, with 
effective areas up to 0.1 km^, the event rate could be as high as ten events per 
year for a 500 GeV LKP. The current limits from AMANDA-II (with data up 
to 2001) is roughly 3000 muons per square kilometer per year from the Sun jl2|. 
This sensitivity is expected to improve significantly with the analysis of more 
recent and future data. 

6.5 e"*" and p from Annihilations in the Galactic Halo 

Charged particles, such as positrons and anti-protons, which are generated in 
dark matter annihilations do not travel in straight lines. Therefore, rather than 
observing a single region, such as the Galactic center or the Sun, the entire 
galactic halo can contribute to the flux of such particles. In this section, we 
will discuss the impact on dark matter annihilations in the galactic halo on the 
cosmic positron and anti-proton spectrum. 

6.5.1 The positron excess 

In 1994 and 1995, the High Energy Antimatter Telescope (HEAT) observed a 
flux of cosmic positrons well in excess of the predicted rate, peaking around 
10 GeV and extending to higher energies [SJ. This result was confirmed by 
another HEAT flight in 2000 t TC^ITC^ . Although the source of these positrons 
is not known, it has been suggested in numerous articles that this signal could 
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Figure 33: The number of events per year in a detector with effective area equal to 
one square kilometer and a muon energy threshold of 50 GeV |3U5| . Contours are 
shown, from bottom to top, for r^i — 0.1, 0.2, and 0.3, where r^i is the mass splitting 
of the LKP and KK quarks over the LKP mass. The expected size of the one-loop 



radiative corrections predict 0.1 



< 



<i 0.2, therefore, the r„i = 0.3 contour is 



shown merely for comparison. The relic density of the LKP's lies within the range 
Q,gih? = 0.16 ± 0.04 for the solid sections of each line. The relic density is smaller 
(larger) for smaller (larger) LKP masses. 



be the product of dark matter annihilations, particularly within the context of 
supersymmetry [^21 EH EHl Ei EZl ElHl EH] and Kaluza-Klein dark matter 

pinii^ . 

If the dark matter is evenly distributed in our local region (within a few kpc) , 
the rate of annihilations may be insufhcient to produce the observed excess. It 
has been suggested, however, that if sufficient clumping were present in the 
galactic halo, that the rate at which such particles annihilate could be enhanced 
enough to accommodate the data. 

Positrons can be produced in a variety of dark matter annihilation modes. 
Direct annihilation to e^e^ is suppressed for neutralinos, but occurs frequently 
for Kaluza-Klein dark matter |150| . Also, annihilations to ZZ or W~^W~ can 
produce positrons with energy of half of the WIMP mass |824j . A continuum 
of positrons, extending to much lower energies, will in most cases also be pro- 
duced in the cascades of annihilation products such as heavy leptons, heavy 
quarks, Higgs bosons and gauge bosons. The spectrum of positrons produced 
in dark matter annihilations can vary significantly depending on the mass and 
annihilation modes of the WIMP. 

As positrons propagate, they move under the influence of the tangled galac- 
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Figure 34: The positron spectrum from neutralino annihilations for the most important 
annihilation modes. Solid lines represent the positron spectrum, per annihilation, for 
X°X° bb, for LSPs with masses of 50, 150 and 600 GeV. The dotted lines are the 
same, but from the process X^X^ ~^ t'^t~ . Dashed lines represent positrons from the 
process XX ~^ W^W^ for LSPs with masses of 150 and 600 GeV. The spectrum from 
XX ~* is very similar. 



tic magnetic fields, travelling in a random walk, and losing energy via inverse 
Compton and synchrotron processes. The diffusion-loss equation describing this 
process is given by 



= V 



b{E^+ , x) 



dn„ 



dE^+ 
dn„+ 



dE, 



Q{E^+,x) 



(143) 



where dne+ /dE^+ is the number density of positrons per unit energy, K{E^+ , x) 
is the diffusion constant, b{E^+,x) is the rate of energy loss and Q{E^+,x) is 
the source term. 

The diffusion constant |489| and rate of energy loss can be parameterized by 



K{E, 



3 X 10^^ 
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and 



b{E, 



10 



(144) 



(145) 



respectively. 6(£'g+) is the result of inverse Compton scattering on both starlight 
and the cosmic microwave background |361| . The diffusion parameters are con- 
strained from analyzing stable nuclei in cosmic rays (primarily by fitting the 
boron to carbon ratio) '374', '375' . 

In equations 1 1 441 and 1 1 451 there is no dependence on location. This is due to 
the assumption of a constant diffusion zone. For our galaxy, the diffusion zone is 
best approximated as a slab of thickness 2L, where L is chosen to be 4 kpc, the 
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best fit to observations [4891 13741 1375| . The radius of the slab is unimportant, 
as it is much larger than the distances which positrons can propagate at these 
energies. Outside of the diffusion zone, the positron density is assumed to be 
(nearly) zero (free escape boundary conditions). For detailed descriptions of 
two zone diffusion models, see Refs. ITMl . 

The effect of propagation on the positron spectrum depends strongly on the 
distance from the source. To compare to the data recorded by HEAT, a quantity 
called the "positron fraction" is typically considered. The positron fraction is 
the ratio of the positron flux to the combined positron and electron fluxes. The 
spectra for secondary positrons, secondary electrons and primary electrons can 
be found in Ref. [557] . 

Figure IJH shows the positron fraction, as a function of positron energy, for 
two scenarios with supersymmetric dark matter candidates. The various lines 
represent clumps of dark matter at different distances from Earth. Note the 
substantial variation in the positron spectrum which results. In all cases, the 
normalization was considered a free parameter. The predicted spectrum is com- 
pared to the error bars of the 1995 and 2000 HEAT data. 

These results show that the spectral shape of the observed positron excess 
can be flt well by dark matter annihilation models. This neglects the issue 
of the annihilation rate (normalization), however. To produce the observed 
excess, a very high annihilation rate is required in the local region (within a few 
kpc). For supersymmtric dark matter, this requires very dramatic dark matter 
substructure j309| . For Kaluza-Klein dark matter, with larger cross sections and 
more favorable annihilation models, it may be more natural to accommodate 
the observed positron excess |3()f)) . 

In the future, new experiments, such as AMS-02 (g^, PAMELA |410j and 
Bess Polar j441) . will refine the positron spectrum considerably. See section 
for more details. 

6.5.2 Anti-protons 

Anti-protons travel much greater distances than positrons before losing their 
energy as they propagate through the galactic magnetic fields. Therefore, the 
dark matter distribution throughout much of the galaxy can contribute to the 
observed anti-proton spectrum 46jl I432| . The measurement of the BESS ex- 
periment finds a cosmic anti-proton fiux of 1.27lg'32 x 10~^cm~^s~^sr~^GeV~^ 
in the range of 400 to 560 MeV. This measurement is difficult to interpret in the 
context of dark matter annihilations due to large uncertainties in the size of the 
diffusion zone and other propagation characteristics |192| . Future experiments, 
especially those with sensitivity at greater energies, will be needed to identify 
signatures of dark matter in the cosmic anti-proton spectrum. For more infor- 
mation on anti-protons from dark matter annihilations, see Refs. [1251 11921 184| . 
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Figure 35: The predicted positron fraction, as a function of positron energy (in GeV), 
for a 150 GeV neutralino which annihilates 96% to bb and 4% to r+r~ (left) or 58% to 
W^W' and 42% to ZZ (right). The solid line represents the distance to the source 
(the dark matter clump) at which the predicted spectrum best fit the data (0.42 and 
0.62 kpc for the left and right panels, respectively). Dotted lines represent the spectra 
for a source at a distance less than found for the best fit (0.23 and 0.19 kpc or 0.20 
and 0.19 kpc for the left and right panels, respectively). For these two lines, the 
is larger by 1 and 4, respectively (1 and 2-a). The dashed lines are the same, but 
for distances greater than found for the best fit (0.85 and 1.3 kpc or 1.1 and 1.6 kpc 
for the left and right panels, respectively) . The normalization was considered a free 
parameter. The error bars shown are for the HEAT experiment. Red (lighter) error 
bars are from the 94-95 flight. The three blue (darker) errors bars between 6 and 12 
GeV are from the 2000 flight. 



6.6 The Role of Substructures 

Annihilation radiation could be enhanced by the presence of substructures in 
the galactic halo. The actual effect depends crucially on the prescription of the 
profile and the spatial distribution of substructures. 

Several groups focused on the signal enhancement due to the presence of 
"clumps" in dark matter distribution, a common feature of N-body simulations. 
The effect of the enhancement of the annihilation radiation on the gamma- 
ray flux has been studied by, e.g. , Bergstrom et al. Calcaneo-Roldan and 
Moore |135| . Tasitsiomi and Olinto |465| . Berezinsky, Dokuchaev and Eroshenko 
2003 [HI and Stoehr et al. gP). 

Recently the problem has been carefully investigated by Koushiappas, Zent- 
ner and Walker j344| , by means of a semi-analytic model of structure formation 
calibrated on high-resolution N-body simulations. The authors concluded that 
previous estimates were optimistic, and that it may be possible for the upcom- 
ing experiments GLAST and VERITAS to detect gamma-rays from dark matter 
clumps only if the neutralino is relatively light, i.e. < 100 GeV. 
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Blasi et al. |114| studied the synchrotron emission produced by secondary 
electron-positron pairs, produced by neutrahno annihilations, in the galactic 
magnetic field. If confirmed, their results would imply a microwave emission 
observable over CMB anisotropics, which is potentially identifiable by its spa- 
tial structure or its radio spectrum. As we mentioned before, the presence of 
substructure is also a possible explanation for the positron excess observed by 
HEAT HUEOII- 

Clumps are not the only substructures that can potentially increase the an- 
nihilation flux. Among other structures considered in the literature, are the 
so-called caustics. In fact, continuous infall of dark matter on our galaxy should 
give rise to ring shaped overdensities, called caustics (see e.g. Ref. |453 |). Un- 
fortunately the prospects for the detection of annihilation radiation from these 
substructures do not appear promising (see Bergstrom, Edsjo and Gunnars- 
son !Hg). 

Apart from galactic substructures, one could ask what the annihilation flux 
from all structures and substructures in the Universe is, i.e. what extra-galactic 
background would be produced by dark matter annihilation. This problem has 
been investigated by Bergstrom, Edsjo and UUio |HZ|, Taylor and Silk |467| . and 
UUio et al. |482) . In particular, the authors of the last reference stressed the 
possibility of observing, for some specific regions of the SUSY parameter space, 
and sufficiently dense substructures, a spectacular feature in the gamma-ray 
spectrum produced by cosmological redshift and absorption along the line-of- 
sight of the gamma-ray line from dark matter annihilation. 

6.7 Constraints from Helioseismology 

The seismic diagnostics of the Sun's interior puts important constraints on the 
internal thermodynamic structure of the Sun. Indeed, such research has led 
to significant improvements in our understanding of microphysics such as the 
equation of state and the opacity calculations, and to a better determination of 
specific cross-sections in the pp chain (see e.g. Ref. |363j and references therein). 

It is intriguing to investigate whether the fact that the Sun evolves in a 
halo of WIMPs affects its internal structure and the details of its evolution. 
Modifying an existing numerical code for the Solar structure, Lopes et al. |363j 
estimated the influence of the WIMP halo on the evolution and structure of 
the Sun, and calculated the deviations of the "modified Sun" with respect to 
the Solar Standard Model and to helioseismic data. They then rejected the 
portions of the WIMP parameter space leading to Solar models in conflict with 
Helioseismic observations. 

Although current measurement do not appear to impose strong constraints 
on dark matter particles (see also Lopes, Silk and Hansen |364| and Bottino et 
al. |127p . it is expected that future helioseismic experiments will be sensitive 
to luminosities from WIMP annihilations in the solar core larger than 10~^ 
times the solar core luminosity (see figure Furthermore, if the dark matter 
density increases toward the central region of our galaxy, as is suggested by N- 
body simulations, stars nearer to the Galactic center would evolve in a WIMP 
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Figure 36: Predicted neutrino-induced muon flux produced by neutralino annihilation 
in the Sun. SmaU squares correspond to models within the phenomenological MSSM, 
triangles correspond to selected benchmark points within the mSUGRA framework. 
Big squares are used to highlight models leading to a local variation of luminosity of 
the solar core larger than 10~^ (which could thus be potentially probed by upcoming 
solar seismic observations). The dotted and dashed curves represent the current limit 
sensitivity of MACRO and the expected sensitivity of IceCube. 



halo which is much more dense, where the effects of dark matter on the stellar 
structure could be of enormous importance. 

Preliminary calculations |362) suggest that evolution times of stars evolving 
in dense dark matter halos are significantly shorter. If confirmed, these results 
would change our understanding of stellar evolution and shed new light on the 
stellar population near the center of our galaxy. 

6.8 Constraints on Superheavy Dark Matter 

In our discussion thus far, we have focused on the detection of weakly interacting 
dark matter particles with masses near the electroweak scale. Here, we will 
extend that discussion to include particles which are much more massive. 

Recently, Albuquerque and Baudis 23 have studied the prospects for the 
direct detection of supermassive dark matter particles. They find that if such 
particles are strongly interacting (simpzillas), masses below ~ 10^^ GeV can be 
probed by current experiments. A superheavy, weakly interacting particle is not 
constrained by this method, however. 

The prospects for the observation of supermassive dark matter annihilation 
from the galactic center are not very promising (see e.g. Ref. HU). Nevertheless, 
portions of the relevant parameter space can be efficiently probed by gamma-ray 
experiments. 

The prospects for the observation of high-energy neutrinos from the Sun 
are potentially interesting |24[ 1168) . For simpzillas, kilometer scale neutrino 
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telescopes, such as IceCube can test a broad range of masses. A signature 
unique to this scenario would be a predominance of tau neutrinos with energies 
above '--^TeV. 

The compilation of results in the work of Starkman et al. |15D| considers 
the constraints on superheavy dark matter found from double-/? decays, cosmic- 
ray detectors, galactic-halo stability, cooling of molecular clouds, proton-decay 
detectors and longevity of neutron stars. 

The constraints derived from old neutron stars is particularly interesting. 
The argument goes as follows: if WIMPs exist they would accrete on neutron 
stars, the same way as they do on the Sun (see Sec. 16. 41 and Sec. l6.7|l |257| . For 
certain regions of the WIMP parameter space, the accretion can be so efficient 
that WIMPs become self-gravitating, then collapse to a mini black hole, which 
finally destroys the star. However, a large portion of the parameter space of 
modern superheavy candidates would escape these constraint, since the collapse 
could be prevented by self-annihilations. 

7 Conclusions 

There is compelling evidence for the existence of dark matter. Although our 
understanding of its nature and distribution is still incomplete, many indepen- 
dent observations suggest that about 30% of the total energy density of the 
Universe is made of some sort of non-baryonic dark matter. We have reviewed 
such observations and discussed how they compare with theoretical predictions, 
and in particular with the results of N-body simulations. 

The dark matter problem is not only relevant to astrophysicists but also to 
the particle and high-energy physics community. In fact, some of the best dark 
matter candidates come from possible extensions of the Standard Model of par- 
ticle physics. There is certainly no shortage of particle dark matter candidates 
found in such models. Among those proposed in literature, we have focused 
on the dark matter particles found in models of supersymmetry (the lightest 
neutralino) and models with Universal Extra Dimensions (Kaluza-Klein dark 
matter). Although many simple models of supersymmetry, extra dimensions or 
other scenarios are widely discussed by the particle and astroparticle communi- 
ties, the phenomenology of the actual physical theory could be more rich and 
complex. Collider experiments are probing significant regions of the parameter 
space of these hypothetical particles. Conversely, a positive astrophysical detec- 
tion of dark matter would provide invaluable information regarding the physics 
"beyond the Standard Model" . 

The astroparticle community has started a vigorous and broad program of 
experiments that may be able to shed new light on the physics and astrophysics 
of dark matter. Before discussing the results of direct and indirect searches, we 
have reviewed the present and future experiments on which they are based. 

Among the most promising dark matter searches appears to be direct detec- 
tion. The current situation is complicated by the claim of a positive detection 
by the DAMA experiment, which have been contradicted by several other ex- 
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periments. It is unclear, but more and more improbable, whether it is possible 
to find a theoretical scenario that accommodates all the experimental findings. 
The much higher (several orders of magnitude) sensitivity of future experiments 
should be able to solve this controversy. 

Indirect dark matter detection via annihilations in the Galactic center region 
is also an exciting possibility, although the prospects for the observation of 
gamma-rays, neutrinos and synchrotron radiation from that direction strongly 
depend on astrophysical parameters, such as the profile of dark matter in the 
innermost regions, which unfortunately are poorly known. Nevertheless, the 
development of next-generation gamma-ray and neutrino telescopes will allow 
us to test many scenarios, especially if effects such as the adiabatic accretion 
onto the central black hole significantly enhance the dark matter density and 
corresponding annihilation signal. If the Galactic center turns out to contain 
less dark matter, observations of dwarf galaxies, external galaxies and local dark 
substructure may play an important role for indirect searches. 

Indirect searches for dark matter through the observation of high energy 
neutrinos produced in dark matter annihilations in the Sun are also promising. 
These rates do not depend strongly on the dark matter halo distribution and are 
thus fairly unambiguous probes of particle dark matter models. Measurement 
of the positron and anti-proton spectra, which are soon to improve dramatically, 
can also provide an opportunity to observe products of dark matter annihilations 
in the galactic halo. 

Collectively, the direct, indirect and collider searches for particle dark matter 
have incredible prospects for discovery in the coming years. We hope that this 
review can be a useful tool in guiding members of the scientific commimity closer 
to the goal of dark matter identification which has eluded us for so long. 
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A Neutralino Mass Eigenstates 

In the Minimal Supersymmetric Standard Model (MSSM), the neutral elec- 
troweak gauginos {B, W'^) and higgsinos (i?", H2) have the same quantum num- 
bers and, therefore, mix into four mass eigenstates called neutralinos. The neu- 
tralino mass matrix in the B-W^-H'^-H2 basis is given by 



Mn = 



Ml 




— Mz cos/3sm9yy Mz sin (3 sin 9 -^^^ \ 

M2 Mz cosPcosO-^fl^ — Mz sin/3cos6'^y 



-Mz cos (3 suvO^ Mz cos [3 cosO^, 
\ Mzsin(3su\dy^ —Mz sin (3 cosOy^ 







(146) 

where Mi, M2 and fi are the bino, wino and higgsino mass parameters, respec- 
tively, 0\Y is the Weinberg angle and tan /3 is the ratio of the vacuum expectation 
values of the Higgs bosons. This matrix can be diagonalized by the matrix, N. 

^diag ^ jsf'iM^oN. (147) 
The masses of the four mass eigenstates are then given by |2(J71 1^ 



eiM^o 



e^M^o = +(ia-ic2)i/2- 



—a 
2 



rC2 



(8a-|C2)i/2 
C3 



e^M^o 



(8a -1^2)1/2 
C3 



1/2 



1/2 



-,1/2 



3a-|C2)i/2 
C3 



(8a-fC2)i/2 



1/2 



i(Mi+M2) , 

(148) 
i(Afi+M2) , 

(149) 
i(Afi+M2) , 

(150) 
i(Mi-HM2) , 

(151) 



where is the sign of the ith eigenvalue of the neutralino mass matrix, and 



C2 



(M1M2 - Ml -t?)~-{Mi+ M2f 



^(Mi + M2)^ + ^(Mi + M2)(MiM2 - M| - + (Mi + M2)^l^ 



l{M,+M2f + \ 
^(Mi cos^ + M2 sin^ 6yv)Ml + ^jlMI sin 2(3 , 



(152) 



(153) 



C4 = -(Micos^ 
1 



M2 sin^ ew)Mln sin 2(3 - MiM2M^ 



-H^(Mi M2)[(Mi + M2)n'^ + (Ml cos^ 6^ + M2 sin^ 6'vi.)M| + ^M| sin 2/3] 



.l(MiM2 - M| - m')(Mi + M2f - -|:(Mi + 1/2)" , 
16 256 



(154) 
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1 



rRe 



21/3 

D = -W^ - 27S^ , U 



-S + i{D/27f/'^ 

1 
3 



1/3 



- 4C4 



n2 2 3 



-C2C4 



(155) 
(156) 



The four masses above are not generally in the order M^o < M^a < M^o < M^o , 

^ ^ Ai X2 A3 A4 ' 

although it is conventional to relabel the states, from lightest to heaviest. 
The mixing matrix, N, is then given by | 2(J7[ I62| 



1 Ml - e^M^o 
'tan6lw Afa - e,;M^ 



(157) 



/i[M2 - e,M o][Mi - e,A// 0] - M| sin/3cos/3[(Afi - M^) cos^ + M2 - e,M 0] 



Afz[-^^2 — EiAf^o] sin0jy[— ^cos/3 + eiAf^^o sin/3] 



(158) 



i4 



-EiM [M2 - eiM 0] [Ml - eiM 0] - M| cos^ P[{Mi - A/a) cos^ 0^^, + M2 - e,M 0] 



Nil 



and 



Mz[M2 - eiM^o] sm9^y[-^lC0sP + eiM^o sin/?; 



(159) 



^^1 = 



N3 
Nl 



-1-1/2 



The lighest neutralino (xi) is a mixture of gauginos and higgsinos: 
The gaugino fraction of Xi is defined as 



(160) 



(161) 



(162) 



and its higgisino fraction as 



(163) 
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Figure 37: Tree level diagrams for neutralino annihilation into fermion pairs. From 
Ref. [3T9j . 

B Neutralino Annihilation Cross Sections in the 
Low Velocity Limit 

In this appendix, we give the amphtudes and cross sections for the most impor- 
tant neutralino annihilation channels in the low velocity limit (the first term in 
the expansion av = a + bv^ + ...). This is sufficient for indirect detection but 
generally insufficient for relic density calculations in which velocity dependent 
contributions are important. For a more complete list, with all S and P-wave 
tree level annihilation amplitudes, see Refs. 15^ ISMHTUS) . 



B.l Annihilation Into Fermions 

Neutralinos can annihilate to fermion pairs by three tree level diagrams |195l 
121:^1 12751 127H| . These processes consist of s-channel exchange of pseudoscalar 
Higgs and Z°-bosons and t-channel exchange of sfermions (see Fig. I37|l. 
The amplitude for pseudoscalar Higgs exchange is given by 

Aa = 4V2 , n 4-(^.M,)^V^r.m.Mr ^'''^ 

Here, niA is the pseudoscalar Higgs mass and Ta is the pseudoscalar Higgs 
width. Tau is the ^'^-neutralino-neutralino coupling and is given by 

Taii = -sin/3Q'/i + cos PS'I^, (165) 

where Q'(^ = iV3,i(iV2a - tan^vyTVij) and S'l^ = A^4a(A^2,i - tan^vi/^i,!). 
TV is the matrix which diagonalizes the neutralino mass matrix in the B-W^- 
Hi-H^ basis, W^a""^ = N'^M^oN (see Appendix . 9w is the Weinberg angle 
and tan/3 is the ratio of the Higgs vacuum expectation values. /lA// is the 
^"-fermion-fermion Yukawa coupling. For up-type fermions, this is given by 

gruf cot /3 



For down- type fermions, it is 



gm/tan/3 
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The amplitude for neutralino annihilation via sfermion exchange to a pair 
of fer miens, fifi, is given by 



= ^ E f + (Kif] ^ + ^x'm^v-^i) ' (168) 



3=1 



where Pj — 1 + (mj./m^)^ — {mf./m^)'^ and the sum is over the six sfermion 
states whieh couple to the final state fermion. The fermion-sfermion-neutralino 
couplings, Xj^j-^ and W^j^, are given by 



and 

where 

and 



X-j, = Xi{IlLef)ij + Zi,kAniiQf)k,j (169) 

w-j^ = Yi{iLRef)ij + Zi^k,i(nLQf)k,j, (170) 

X, = -gV2[T^{fi)Nl^ - tanOwinif,) - e{h) ) iV* i ], (171) 



Fi = g^/2 tanOw e{fi)Nl^. (172) 
For final state up-type quarks, 



9_ 

^/2mw± sin/3 



= - ^ . . e,,, n:,. (173) 



For final state down- type quarks. 



Zi,j,i = ^ ^0^,. ^3,1- (174) 

V^myy± COSp 

And for final state leptons, 

Zi,j,i = ^©iJ ^3,1- (175) 

^/ZmY/± cosp 

Here, Tz{fi) and e(/i) are the weak hypercharge and electric charge of the final 
state fermion. iV, again, is the matrix which diagonalizes the neutralino mass 
matrix. 0/'s are the appropriate 6x6 sfermion mass matrices and TiL,R are 
left and right projection operators: 

1 
Hl = I 1 I , (176) 
1 

1 
Hii = I 1 I . (177) 
1 
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Lastly, the amplitude for neutralino annihilation to fcrmions via Z exchange is 
given by 



Az = 2V2^^0'/,fr3(/,)^^, (178) 
cos^ t^w "f^z 

where T3(/i) is the weak hypercharge of the fermion. The couphng 0'{ \ is given 

by i(-^3aA^3*i + ^4,i^4,i)- 

Summing these three contributions to the amphtude, we can calculate the 

cross section for this process: 

ov{xx kn)v-.o = j^i^i^Aixx ^ .kn)+Aj{xx ^ kn)+Azixx ^ /./.)i', 

(179) 

where /3/ = y^l — m^/m^. c/ is a color factor which is equal to three for quark 

final states and one otherwise. 

We emphasize that all tree level (low velocity) neutralino annihilation dia- 
grams to fermion pairs have amplitudes which are proportional to the final state 
fermion mass. For sfermion and Z° exhange, this is because the Z°-fermion- 
fermion and neutralino-fermion-fermion couplings preserve chirality. For psue- 
doscalar Higgs exchange, the amplitude introduces an explicit factor of the 
fermion mass in the Yukawa coupling. We also note that the Yukawa coupling 
which appears in the psuedoscalar Higgs exchange amplitude is proportional to 
tan/3 for down- type quarks and cot /3 for up- type quarks. The net result of these 
observations is that neutralino annihilation into fermions will be dominated by 
heavy final states, bb, t~t^ and, if kinematically allowed, tt. Furthermore, if 
tan/3 is large, bottom-type fermions may dominate over up- type fermions, even 
if less massive. For example, annihilations to bb may dominate over tt, even for 
heavy neutralinos. 



B.2 Annihilation Into Gauge Bosons 

Generally, neutralinos can annihilate into gauge boson pairs via several processes 
(see Fig. OHl) [T^ ITfSl BUI ITOj . In the low velocity limit, however, only t- 
channel processes via chargino or neutralino exchange are non-vanishing. 

In the low velocity limit, the amplitude for neutralino annihilation to W^- 
pairs is given by 

2 1 

Aixx - w+w-).^o = 2V2 Pwg' [(^t)' + J'] jr^ (180) 

n—l ^ 

where f3w — y^l — m^/m^ and P„ = 1 4- ("i^± /m^Y ~ {mw /fn^Y . The sum 
is over chargino states. Of „ and 0{^„ are the neutralino couplings to charginos 
given by ^^^4,1^2% + ^2,1^1% and ^^^3* it^2,n + ^^2* i^^i,"' respectively, where 
iV, again, is the matrix which diagonalizes the neutralino mass matrix. The V^s, 
and f/'s are components of the chargino mass matrix, in the basis 
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Figure 38: Tree level diagrams for neutralino annihilation into gauge boson pairs. 
From Ref. PT^ . 



sm 0_ cos ( 



u=['Tv -fr^ ) (181) 



and 



where 



and 



( cos0+-sm 0, 
' sm 0+ cos (p- ' 

(M| - + 2m\, cos 2/3) ^ ' 

r- (u COS /3 + Mq sin /?) 

tan20, - ^^-^ (M| "-24 cos /3) - ^^^^^ 

The amplitude for annihilations to Z'^-pairs is similar: 

Aixx - ^°^°).^o = 4x/2/3z-|^ E (Oi,^)' (185) 

cos (7^ -Tj^ 

n— 1 

Here, Pz = y^l — jti^ / , and P„ = 1 + (rn^^/m^)^ — {mz /m^)'^ . The sum is 

over neutrahno states. The coupHng O"^ is given by i(— A^ajTV^ „ + A'^4_iiV|„). 
The low velocity annihilation cross section for this mode is then given by 

avixx ^ GGU^o - ^ (186) 

where G indicates which gauge boson is being considered. Sg is a statistical 
factor equal to one for W^W~ and two for Z^Z^. 

It is useful to note that pure-gaugino neutralinos have a no S-wave annihi- 
lation amplitude to gauge bosons. Pure-higgsinos or mixed higgsino-gauginos, 
however, can annihilate efficiently via these channels, even at low velocities. 
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Figure 39: Tree level diagrams for neutralino annihilation into a Z and a Higgs boson. 
From Ref. ETHl- 



B.3 Annihilation Into Higgs Bosons 

There are many tree level diagrams which contribute to neutralino annihilation 
into Higgs boson pairs or a Higgs boson and a gauge boson (see Figs. 15^ BUI and 
mTl |lMII?7^Bnil¥)Kll^l^lT7^ . 

In the low velocity limit, the amplitude for neutralino annihilation to a Z° 
and a light neutral Higgs, /i'^, is given by 



-4(xx - 



^ sin(/3 — a) 



a 



4 



mz cos dw 

2cos(q! - /3)Tai,i 



/i l.n 



(187) 



TO^cosS^ 4 — m^/rri^ _|_ j p^^^/. 



where Ya is the pseudoscalar Higgs width and T^i^i is the — Xo ~ Xn 
Yukawa coupling (sec below). The couplings, O'/^, is given by N3^i{N2.n — 
tan6'H'^i,n)/2+iV3,„(7V2,i-tan6'vy7Vi,i)/2 and P„ = l+{m^^/m^)^-^{mz/m^y- 
^{rrii^/m^)'^ . tan/3 is the ratio of the Higgs vacuum expectation values and the 
mixing angle, a, is related to beta by 



and 



sin 2a 



cos 2a 



sin 2/3 



cos 2/3 





+ ml 














m\ 


-ml 




-ml 



(188) 



(189) 



The three terms of Eg. 11871 correspond to neutralino, and A° exchange, from 
first to last. 

The expression for neutralino annihilations to a Z*^ and a heavy Higgs boson, 
iJ", is the same, but with sin(/3— a) and cos(q; — /3) exchanged, and the couplings 
and masses of /i° replaced by the couplings and masses of These Yukawa 
couplings are given by 



sin a Q'( „ + cos a S'" „ 



(190) 
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Figure 40: Tree level diagrams for neutralino annihilation into a W and a Higgs 
boson or a pair of Higgs bosons. From Ref. |319| . 



and 

Th i,n = - COS a Q'l^ + sin a S'l„. (191) 

Here, S'l„ = iV44(A^2,„ - tan0wA^i,„)/2 + A^4,„(iV2.i - tan0vviVi^i)/2. g'/„ is 
defined above. 

The amplitude for annihilations to a and a charged Higgs boson is given 

by 



2 



2 mw 

n— 1 



Pn 

rriyTAii 



(192) 



2 mw 

n—l 



Pn 



mvi/(4 — m\/mV) 



where P„ = 1 + im^±/m^)'^ — i(TO^±/m^)^ — ^(tow/tojj,)^. 0{*„ and Of^ 
are couplings given earlier in this appendix. and Q'^^ the chargino- 

neutralino-charged Higgs couplings, given by 



Qnm = cos/3 



and 







nm 



sin/3 



(iV2n + tan 6*14/ iVi„)t/2« 



(193) 



(194) 



where each of the quantities used have been defined earlier in this appendix. 
The first and second terms of Eq. 11921 correspond to chargino exchange. The 
third term corresponds to pseudoscalar Higgs exchange (there is no low velocity 
contribution from scalar Higgs exchange). 

Finally, the low velocity amplitude for neutralino annihilation into one neu- 
tral Higgs boson and one pseudoscalar Higgs boson is given by 
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Figure 41: Tree level diagrams for neutralino annihilation into a neutral Higgs boson 
and a pseudoscalar Higgs boson. From Ref. |319| . 



p mz sin(a+/3) COS2/3 Taoo _ cos(a-/3)05^oo m\-ml 
2 cos 0w 4— m-^/m^ 2 cos-^ 9w 

Here, P„ = 1 + {m^^/m^Y ~ ^{'i^a/^x)'^ ~ ^i^h/m--^)^ . The other quantities 
have been defined earlier in this appendix. Again, the amplitude for the anal- 
ogous process with a heavy rather than light Higgs boson in the final state is 
the same, but with sin(a + f3) and cos(q; — /3) exchanged and the light Higgs 
couplings and masses replaced with those for the heavy Higgs boson. 

In the low velocity limit, there is no amplitude for neutralino annihilations 
to H+H-, h°h°, H°H°, A°A° or Z°A°. 

The low velocity cross section for neutralino annihilation via any of these 
modes is 

avixx -> XY),^o = T^^IMXX ^ ^^.-0!', (196) 
where X and Y are labels referring to the final state particles. 

B.4 Annihilation Into Photons 

Although there are no tree level processes for neutralino annihilation into pho- 
tons, loop level processes to 77 and jZ'^ are very interesting, as they may 
provide a spectral line feature observable in indirect detection experiments. 

In Fig. l42[ all of the one- loop diagrams are shown for neutralino annihilation 
to a pair of photons. In Fig. 031 the corresponding diagrams to a photon and a 
Z° are shown. We do not include the corresponding amplitudes or cross sections 
here. For those results, see Ref. [7H| and Ref. [ISD] for 77 and 7^° final states, 
respectively. Also see Ref. |271| . 
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Figure 42: Diagrams contributing, at one loop level, to neutralino annihilation into 
two photons. From Ref. {79 f. 
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Figure 44: Tree level Feynman diagrams for neutralino-quark scalar (spin-independent) 
elastic scattering. From Ref. |319| . 



C Elastic Scattering Processes 
C.l Scalar Interactions 

Consider a WIMP with scalar interactions with quarks given by 

Cscalar = agXXQq, (197) 

where Uq is the WIMP-quark coupUng. Then the scattering cross section for 
the WIMP off of a proton or neutron is given by 

r^^>' dajv ^ 0) _ 4ml .2 ...o^ 

(^scalar — / jl ^12 ~ Jp.nJ 

Jo d\v\^ TT 

where v is the relative velocity of the WIMP, is the reduced mass of the 
nucleon (m^ ~ "ip „ for WIMPs heavier than ~ 10 GeV) and /p „ is the WIMP 
coupling to protons or neutrons, given by 

q—u.d,s q—c.b.t 

where f^^^ = 0.020 ± 0.004, f^} = 0.026 ± 0.005, f^^] = 0.118 ± 0.062, f^'J = 

0.014 ± 0.003,4'^^ = 0.036 ± 0.008 and f^"^ = 0.118 ± 0.062 f^^^ is 

related to these values by 

/(f^") = l- E frV^- (200) 

q—u,d,s 

The term in Ea. ll99l which includes f^^p^ results from the coupling of the WIMP 
to gluons in the target nuclei through a heavy quark loop. The couplings of 
squarks and Higgs bosons to heavy quarks leads to a loop level coupling of the 
WIMP to gluons [2761 161L I323| . Such diagrams are shown in Fig. 03] 

To attain the scalar cross section for a WIMP scattering off of a target 
nucleus, one should sum over the protons and neutrons in the target: 

^""^ ' Z/p + (A - Z)/„ ) , (201) 



116 




Figure 45: Feynman diagrams for neutralino-gluon scalar (spin-independent) elastic 
scattering. Notice that no tree level processes exist. From Ref. |319| . 

where Z and A — Z are the numbers of protons and neutrons in the nucleus, 
respectively. 

The above expression is valid only at zero momentum transfer between the 
WIMP and the nucleon. For finite momentum transfer, the differential cross 
section must be multiplied by a nuclear form factor. The appropriate factor, 
called the Woods-Saxon form factor, is given by 221 

i^(Q)=(^^^)'exph(qs)2], (202) 

where ji is the first spherical bessel function and the momentum transferred 
\b q — y/sruNQ. Ri is given by y/R^ — 5s^, where R and s are approximately 
equal to 1.2 fmA^/^ and 1 fm, respectively. 

Although less accurate than the Woods-Saxon form factor, the following 
simple form factor is sometimes used in its place |17l I24()| : 

F{Q) ^ exp[-Q/2Qo]. (203) 

Here, Q is the energy tranfcrred from the WIMP to the target and Qo = 
1.5/(TOA,i?g) where Rq = 10"" cm [0.3 0.91(mN/GeV)i/3]^ 

In the context of ncutralino scattering, the value of Oq can be calculated 
from the parameters of the MSSM [2iHI HSi IMI 1111 ■ Following Ref. W, a. 
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is in this case given by: 



2(m2^ - ml) 



Re [(X.) iYT] 



2{ml^ - m2 ) 



Re [{W^) {V^y 



gmq 
ArnwB 



Re {6i[gNi2-g'Nii])DC[- 



Re {52[gm2-9'Nn]) 



(204) 



where 



Y, 



gmqN*r, 



* / TV 

Vi2eig ^^11, 



2mwB 

* (Vi 'AT , rj. j.r \ , * gnigNl, 

Vn [J9N,^+9T,,N,2)+rh2^^^ 



5-i 

2mwB 



;^ 5 '?22ei5 ^11, 



2mwB 



2mwi3 

where and denote hypercharge and isospin, and 

Si = NisiNu), 62 - ^i4(-iVi3), 
-B = sin/3(cos/3), C = sina(cosa), 
D = cos «(— sin a) . 



(205) 



(206) 



Here, i = 1 for up-type and 2 for down-type quarks, a is the Higgs mixing 
angle. mii,m2i denote elements of the appropriate 2x2 squark mass matrix. 
Nin are eilements of the matrix which diagonalizes the neutralino mass matrix 
in the B-W^-Hi-H2 basis (see Appendix fX)l . rj is the matrix which diagonalizes 
the appropriate squark mass matrices. 

To crudely estimate what scale we expect for the scalar cross section between 
a neutralino and nucleon, we can carry out a back-of-the-envelope estimate. 
Considering a gaugino-like neutralino, we see that as and S2 both vanish, so 
do most of the terms in Eq. 12041 We are left with a neutralino-quark coupling 
of ttq ~ A/m'^, where A is the product of the various order 1 couplings, mixing 
matrix parameters, etc. which contribute. For a typical case, A might be 
~ 1G~^ or so, although it can vary a great deal. Inserting this coupling into 
Eqs. 11981 and 11991 we estimate a neutralino-nucleon scalar cross section of ~ 
A^mp/m~, which is roughly 10^^ picobarns, for TeV mass squarks. These results 
can vary dramatically, however, depending on the characteristics of the model 
being considered (see Figs. 1211 and I22|l . 

We can contrast this with the much larger neutralino annihilation cross 
sections. Considering again a gaugino-like neutralino, its amplitude for annihi- 
lations into bb via psuedoscaler Higgs exchange (see Eq. I164f) is roughly Aa ~ 
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a a 



Figure 46: Tree level Feynman diagrams for neutralino-quark axial- vector (spin- 
dependent) elastic scattering. From Ref. |319| . 

mbtsin.(3\/Jh/'my^± where fh is the higgsino fraction of the WIMP. The annihi- 
lation cross section fEa. 1179(1 is then roughly a ^ Sm^tan^ /?//j/1287rm^m^± . 
For even a very small higgsino fraction, say 1%, and a 200 GeV neutralino, 
we find a cross section of ~ 10~^ picobarns for small values of tan (3 and a few 
picobarns for tan/3 — 30. 



C.2 Axial- Vector Interactions 

Next, we consider a WIMP with axial-vector interactions with quarks given by 

= dq X7''75X qit^lbq, (207) 

where dq is the generic coupling. 

For such a WIMP, the spin-dependent scattering cross section can be written 
as HHH] 

W ^ ^'''^ 

where v, again, is the relative velocity of the WIMP, and r(w^) is the scattering 
matrix element. This expression can be integrated over the Boltzman veloc- 
ity distribution of halo WIMPs to arrive at an average elastic scattering cross 
section. At zero momentum, the matrix element, T(v'^), is given by 

WW = ^^^^^\iduAP+ddAP+d,AP) <Sp> -fKAC+drfA','-f4A;^) < S„ > \^ 

(209) 

where J is the nuclear spin and the A's are the fraction of the nucleon spin 
carried by a given quark. Their values are measured to be Af^ — A^* — 0.78 ± 
0.02, A^ = AS = -0.48 ± 0.02 and Af = A^' = -0.15 ± 0.02. < Sp > and 
< Sn > are the expectation values of the total spin of protons and neutrons, 
respectively. Notice that for target nuclei with even numbers of protons and 
neutrons, there is zero total spin, and the cross section vanishes. 

The values of < S'p > and < Sn > depend on the nucleus being considered. 
For ''^Ge, the interacting shell model finds < Sp > and < S'n > to be 0.011 and 
0.468, respectively. For ^^Si, they are given by -0.0019 and 0.133. For ^^A, they 
are 0.3430 and 0.269. And for ^^K, they are -0.184 and 0.054 [5^5] . 
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For non-zero momenta, a more complex form of eauation l2U9l is needed. This 
equation is given by 



J 

?0/„ .2 



<S'p-5„ >Fi(z;2)|2, (210) 



where the i^'s are nuclear form factors given by 
and 

2„2 i%2 



F'{v^) ~ exp(-rfi;7?i^ + icv/h) 



(211) 



(212) 

where Tq and ri are parameters which depend on the nucleus being considered, 
with typical values of 1.3 — 2.1 fm^^. 

Again, within the context of neutralino scattering, the value of d2 can be 
calculated from the parameters of the MSSM ,266, ,223. ,2751 EZi EIH 
Following Ref. j209| . ^2 is in this case given by: 



\YA^ + \xf 



4(m2 



4m| cos^ 6w 



|iVi3|'- |iVi4|' 



2i 

Ik 

2 ' 



where the quantities used are defined in Appendix Kill 



1X^:1 



\W,\ 



(213) 



C.3 Vector Interactions 

As a third case, consider a WIMP with vector interactions with quarks, given 

by 

^ec = \ xit^x q^q- (214) 

Here, bq is the WIMP-quark vector coupling. In this case, the contributions of 
each quark in the nucleus add coherently and large cross sections result for large 
nuclei. The WIMP-nucleus cross section in this case is straight forward 



647r(m^ + rn]\j)^ ' 



(215) 



where Bn is simply bpf = 2Zbp + {A~ Z)b„ 



As an example of a WIMP with vector interactions, consider a Dirac neu- 
trino. In this case, bq — GpiT^ — 2eq sin^ 6'w)/\/2, where Gp is the Fermi 
constant, and are the weak isospin and electric change of the quark g, re- 
spectively, and 0w is the Weinberg angle. Summing over the quarks in a proton 
or neutron, we get bp = Gf(1 — 4sin^ ^w)/(2\/2) and bn = —Gf/{2^/2). Since 
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4sin^ — 1, the neutron-neutrino cross section is much larger than the anal- 
ogous proton-neutrino interaction. The Dirac neutrino-neutron cross section is 
then given by u^^n — G|,m^m^/(5127r(m^ -I- m„)^). A cross section of this size 
has been ruled out by direct scattering experiments, except perhaps in the case 
of a very light WIMP. Similar calculations show that other similar particles, 
such as sneutrinos, are also excluded by this method |225j . 

Neutralinos, being Majorana fermions, do not have vector interactions with 
quarks. 
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